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Bewegungen eines Schwingers unter dem Einfluß von 
a. achwarzeWeißsReselungen... 000 008 
NT "Bewegungen eines Schwingers von einem Freiheitsgrad; Regelung mit 
00.0... ‚Stellungszuordnung ohne Schaltverschiebungen. . 
3 Von I. Flügge-Lotz in Paris und K. Klotter in Karlsruhe. !) 


RN Sogenannte Schwarz:Weiß-Regelungen (auch Ja-Nein-Regelungen genannt) werden wegen ihres einfachen er 
> konstruktiven Aufbaus vielfach stetigen Regelungen vorgezogen. Wir beschränken uns in der vorliegenden 
. Arbeit darauf, eine ‚ideale‘ Regelung zu untersuchen, bei der das: Stellglied dem Kommando, das durch 
. die „‚Regelungsfunktion‘“ gegeben wird, sofort gehorcht. Die Schwingungseigenschaften des Stellgliedes und 
die Unvollkommenheiten im Gestänge, Fühler und Mischaggregat werden also zunächst vernachlässigt. Dem- 
‚nach erfolgt die Bewegung des Stellgliedes, solange auch für die Praxis die Abweichnugen von der „idealen“ 
. Regelung vernachlässigt werden können, im Takt der Bewegung des Hauptsystems. Das Ziel unserer Unter- 
" suchung ist, eine Übersicht über den Einfluß der-System- und Regelbeiwerte auf das zu erreichende Abkling- 
verhalten zu gewinnen. ?) 
= 'So-called ““on-off” regulations are often preferred to continuous regulations owing to their simple construc- 
3a ER tion. In the present report, we confine ourselves to investigate am “deal”. regulation, i.e. the regulating 
SE) device instantly obeys the command given by the “function of regulation’’. T'he elastic properties of the regu- \ 
ae > Iating device and the defectiveness of the other parts are conseqwently neglected. Thus the motion of the regu- 
er -, . . lating device synchronizes with ihe motion of the main system provided that ihe deviations from the “ideal” 
A £ . . regulation can be neglected in the practice likewise. T'he aim of our research is to gain information about the 
u : influence of the data of the system and of ihe regulation. on the obtainable effect of the diminishing of the .os- 
RIESE: cillations. £ : | 3 
Les reglages dits «noirs-blancs) sont souvent pr£feres aux reglages continus gräce 4 leur construction 
| simple. Dans la publication presente les auteurs se bornent & un reglage «ideal ), c’est-ü-dire le systeme de 
a reglage obeit im&diatement a la « fonction de reglage). Les propri ts lastiques du sysiöme de röglage ei la 
defectuosite des autres parties sont megligees. Ainsi le mouvement du systeme de röglage est en correspondance 
avec le mouvement du systöme principal, & condition que les deviations du röglage ideal puissent ätre megligees 
dans la pratigue. Le but. de la recherche est de trouver l’influence des donndes du sysiöme et du reglage sur 
m la diminution des oscillati ons. E 


Bi Peryınpyiomne MmexanuaMmbl, pa6oTamınne NO IPMHIUIY ‚„UePHEIH-OeJIBIH‘“ Mm „‚Aa- 

Sn f ‚HET‘, yACTO NPEAHOYHTAMT CH, BCHIEACTBUE UX 6o1ee HPOCTOU KOHCTPYKIHH, HeITPepHBHO 
peryılmpyloımmm mexanusmam. B naHnoü pa6oTe AaBTOPbI OTPAHHYUBAMTCH UCCHEeNOBAHUEHM 
„HNeANBHO‘ PeryAaupylolmero MEXAHU3MA. PeTyMIMPmINaA CHCTEMA KOTOPoTO CPasy HOAUH- 
Haerca, „peryanpysomei dyukunn“ Taxum o6pasoM, sec CHayama He YYHTBIBAIOTCA 
yNpyrue cBOÄCTB yaperyınpmmmeü. CHCTeMbI U HeCOBEPLIeHCTBA APYTUX yacTef MAXaHH3Ma. 
 BenencTBue 9TOTO, ABUSKeHHE peryiımpyloineii CHCTEMBI COOTBETCTBYeT ABHKEHNIO TIIABHON 
CHCTEMbI, EeCJIM HA IIPARTUKE MO3KHO npeHedpe4b. OTKJIOHeHHAMN OT ‚„‚UNeANBHOTO‘‘ Pery- 
- IHpyrolmero MeXaHusMa, Aannoe wecrexoBaHne CTABUT cEebe MEI BEIACHUTB BIIUSIHHE TIO- 
CTOAHHBIX CUCTEMBI U PeryAupyiIMeTo, MEXAHU3MA HA 3ATyXAaHme KONe baHud. 
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1. Die Sehwarz-Weiß-Regelung (Systemgleichung, Reglergleichung, Bewegungsgleichung, 
Regelfunktion, Systeme vom Typ A und B, Schaltstellen vom Typa und b). 


Die geregelten Bewegungen eines Fahrzeugs (das selbst als System von einem Freiheitsgrad 
betrachtet wird) werden durch zwei Gleichungen bestimmt: a) die Bewegungsgleichung des 
Systems (‚‚Systemgleichung“), b) die Reglergleichung. Die Systemgleic hung kann 

“meist mit mehr oder weniger großen Vernachlässigungen auf die Form - 


N De a Eee th) 


2) In dieser Zeitschrift, Bd. 25/27, S. 97...113, hat I. Flügge-Lotz unter dem gleichen Obertitel zusammenfassend über 
Ergebnisse einer Reihe bislang unveröffentlichter Arbeiten verschiedener Verfasser aus den Jahren 1942 und 1943 berichtet. 
Die erste dieser Arbeiten, auf denen der genannte Bericht fußt, wird hier nunin ausführlicher Fassung vorgelegt, damit 
5 neben den vollständigen Ergebnissen auch die verwendeten Methoden im einzelnen gezeigt werden können. Die Verfasser 
! sind sich dabei bewußt, daß die Darstellung noch nicht endgültig sein Kann. 

2») Das gleiche Problem, wenn auch in etwas engerem Rahmen (nur System A) wurde schon von H. Bilharz (Über 
eine gesteuerte eindimensionale Bewegung, Z. angew. Math. u. Mechan. Bd. 22, S. 206 bis 215) behandelt. Die Unter- 
suchungsmethoden weichen voneinander ab. H.Bilharz hat besonderen Wert auf die mathematische Darstellung gelegt, 
während wir uns bemüht haben, die mechanische Seite des Problems in den Vordergrund zu rücken. Nachdem die Ar- 
heiten unabhängig voneinander entstanden waren, zeigtengemeinsame Diskussionen völligeÜbereinstimmung der Ergebnisse. 
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gebracht werden. Darin bedeutet g die die Bewegung kennzeichnenc 
einer Fahrzeugkursregelung also z. B. die Abweichung vom Kı 
"der Kursregelung den Ruderausschlag —; der Faktor N wird gelegentlich als „| 

bezeichnet up a0 00, W EB BERATER > 
Erfolgt die Regelung der Bewegung durch eine stetige Änderung des 


die Reglergleichung allgemein 


x 


an 7 Be mr Zug) Anne en A 
Bi "N S y ö | w=0 „ BR rs 5 DNA , 
7... Die Anzahl der Glieder in den Summen hängt dabei von der Anzahl der gefohlten Ab, I 
0... Jeitungen 9, @, ... und von dem Aufbau des Reglers ab, durch den die Zahl der Ableitungen B ag 
SFR . bestimmt wird. E CH | ur 
N In (1.2) ist z.B. die Gleichung der einfachen Stellungszuordnung en 
2.0. RS RE: ae ee 2 TEAK BE 
0° oder der einfachen Laufgeschwindigkeitszuordnung ERERENF 

RR B = ap +a,9 DEWERT AU ae EEE TR (1.2b) 


(bei Benutzung von zwei gefühlten Größen) enthalten. ' re + | 
Durch Einsetzen der nach ß aufgelösten und so oft als nötig differentiierten Gl. (1.1) in 


Pr 


Br | die linke Seite von (1.2) entsteht für die Koordinate p eine homogene lineare Differentialgleichung 

5 höherer Ordnung mit konstanten Koeffizienten. Sie beschreibt die zu untersuchende, geregelte 
2 Bewegung. Die Wurzeln ihrer Säkulargleichung geben z.B. Auskunft über das Stabilitäts- . 
ih; verhalten. 

B | . In dieser Arbeit soll im Gegensatz zu einer stetig arbeitenden Regelung eine sog. Schwarz- 


-Weiß-Regelung (auch Ja-Nein-Regelung genannt) untersucht werden. Das Kennzeichen einer 
‘solchen Regelung besteht darin, daß die Stellgröße ß sich nicht stetig und proportional einer 
linearen Kombination der gefühlten Größen 9,9, $, .... ändert, sondern nur zwei feste Werte +ß, 
oder —ß, annehmen kann. Die Reglergleichung lautet in diesem Fall, wenn eine lineare Kom- 
bination der beiden Größen p und g, nämlich F= 9,9 + 02, gefühlt wird, 


Ber-ßsen (ie +39 2 ann (1.3a) | 


| BEHBSDIOF ORTEN Ne (1.3b). 

Die Funktion F = 0,9 + 059, deren Vorzeichen sowohl in (1.3a) wie in (1.3b) über das Vor- 

zeichen der Stellgröße ß entscheidet, heiße im folgenden die Regelungsfunktion. 

Wenn die Regelungsfunktion F durch O0 geht, wird die Stellgröße ß in unstetiger Weise 

von ihrem pesitiven Wert ß, auf ihren negativen Wert —ß, (oder umgekehrt) „umgeschaltet“, 

Die Stellen des Nulldurchgangs der Funktion F sind .also die ‚‚Schaltstellen‘“ der Regelung. 

n Solange (wie dies in unserer Betrachtung der Fällist) für das Umschalten der Nulldurchgang 
der Regelungsfunktion maßgebend ist, darf man an Stelle von 


F= 0,9 +09 


oder 


die Funktion 


Fer ERTL SR ee (1.4) 
betrachten (wo 0 = o3;/o, ist). In der Form (1.4) werden wir in dieser Arbeit die Regelungs- 
funktion stets benutzen. 

Ist 0 > 0, so eilt die Regelungsfunktion F dem Ausschlag p um die Zeit 2, vor, ist e< 0, 
so eilt sie um t„ nach (s. Bild 1). 

Wir zählen die Zeit übrigens so, daß von jedem Schaltpunkt aus eine Zählung mit t= 0 
beginnt. Die Werte g(0)=g, und 9 (0) —=9, heißen auch die Schaltwerte. Wegen 
F (0) = 0 besteht zwischen ihnen die Beziehung 

Po 


En a TEEN ee 1.4a). 
[ I (1.4a) 


Die Länge des Intervalls bis zum nächsten Schaltpunkt sei 2}. 

Ein System, dessen Reglergleichung wie (1.3a) lautet, wird im folgenden stets als S ystem 
A, eines, dessen Reglergleichung wie (1.3b) lautet, als System B bezeichnet. Man beachte 
also, daß zwar mit der Anderung des Vorzeichens der Regelungsfunktion F das Vorzeichen der 
Stellgröße ß sich ändert, daß es dabei aber zwei Möglichkeiten von Zuordnungen gibt; in 
System A ist 
sen ß= —sgn , 
in System B ist dagegen 

sgnß= +senF. 


Be re RA 
e gel ten Be wegungdes Systems RN Dabei gilt das Be 
System vom Typ 4, das untere für eines vom Typ B. en -N Er ist 6: 
tegelmomentes : Be. 


wi nn m=-F(snMM r 
auf der Techn Fa der GE (1.5a) steht; M ist also eine wesentlich PERmE Größe. 
N Den Verlauf der Regelfunktion F und des Sn 


 Regelmomentes m über der Zeit zeigt Bild 2a für 


ein System A, Bild 2b für ein System B. Wir 


\ Base weiterhin. ‚eine Schaltstelle, De der das Regel- 


Bild 1. ‘Voreilende (e>0) und nacheilende ve < 0) Bild 2a und 2b. Verlauf der Regelfunktion F und 
Regelungsfunktionen. des Regelmomentes m für System A und System B, 


moment m von seinem negativen Wert auf den positiven springt, eine Schaltstelle a, eine 


solche, bei der das Regelmoment vom positiven auf den negativen Wert springt, eine 


 Schaltstelle b. Nach dem vorher Gesagten ist klar, daß im System A bei einer Schaltstelea 


die Steuerungsfunktion F von positiven nach negativen Werten durch Null geht (F >0), bei 


- Schaltstellen b dagegen von negativen nach positiven Werten (F < 0), während diese Zuordnung 
im System B umgekehrt ist. 

In der Form (1.5a) der Bewegungsgleichung treten vier Parameter, a,b’, c' und M, auf; 
ferner haben alle Glieder die Dimension eines Momentes. (einer verallgemeinerten Kraft). Nach 


Baeion durch den Trägheitsfaktor (das Trägheitsmoment) a’ nimmt die Gl. (1.5a) die Gestalt 


$+2Do9 +o®p= 7 (sgn F)b a a PN (1.5b) 


"an; die Zahl der Parameter ist auf drei, D, , b, zurückgegangen; jetzt haben alle Glieder der 


Gleichung die Dimension einer Winkelbeschleunigung. Es bedeuten dabei (wie in der Schwingungs- 
lehre üblich) 


und b=M/ja'. 


a ae 
a’ ver ” 
Das obere Vorzeichen auf der rechten Seite gilt BAR für Systeme vom Typ A, das untere 
für Systeme vom Typ B. 
Gelegentlich ist statt der Form (1.5b) der Bewegungsgleichung noch eine weitere Form 
AR sie entsteht a aus (1.5b) nach Division durch &® und lautet 


ee 
DE +p=Fisnf) 5 EDER Er SIESCH 


Hier haben nun alle Baer die een der RE p selbst. Schließlich ließe sich mit 
einer ‚neuen unabhängigen Veränderlichen 7= wt, die zur Folge hat, daß die Ableitung 
1d d 
@di dr 


wird, die Differentialgleichung (1.5c) auf die Form 
b 
po’ +2Dp +9 = Feen ER N ee 


bringen, wenn die Ableitungen nach 7 durch Striche beieichnet werden. 

In allen numerischen Rechnungen werden wir weiterhin den Parameter b/o® = M/e' 
gleich Eins setzen. Das bedeutet, daß wir über die Maßstäbe verfügen, in denen wir die Ergeb- 
nisse aufzeichnen. Der Parameter b/o® bedeutet nämlich jenen Ausschlag (Winkel) 9 = 94, den 
das.Regelmoment M an der Rückstellfeder c’ erzeugt. Einheit der Veränderlichen ist dann 

Eike 


wei $ = REEN 


820 Fiügge-Lotz und Klotter, 


Y 


wir ?n 
Dr 


dieser „statische“ Auschlag g4. Oder anders ausgedrückt: 


a 


„ Als z bt änglg 

an Stelle der Größe p die neue Größe $—= p/pa benutzt. So erhält die Differe: 

schließlich die letzte Form ER EI, EN: 
+2 DE +HF= FM ......22 22. lböe) 


Von den vier Parametern a‘, b‘, e' und M, die in (1.5a) auftreten, oder den drei wesentlichen 
Parametern D, ® und 5/o® = M/e', die in (1.5c) auftreten, sind durch die genannten Maßnahmen, 


ge nämlich erstens andere Zählung der Zeit (bzw. Verwendung von Ableitungen p/o, Pl? usw.), 
a "zweitens Verfügung über die Maßstäbe der Veränderlichen g und @/o (bzw. andere Zählung des 
E: '  Ausschlags) zwei Parameter, nämlich & und b/o entfernt worden, so daß als einziger Parameter 
a die Dämpfungszahl D übrig bleibt. Für sie ist, falls nichts anderes bemerkt wird, in allennume- 
ERS, rischen Rechnungen D = 0,1 gesetzt. BEN / 
RR Zur besseren Wahrung der Anschaulichkeit benutzen wir die Bewegungsgleichung jedoch 


nicht in der Form (1.5e), sondern bleiben bei der Form (1.5b) oder (1.5c) stehen. Wir rechnen 
also mit der Zeit t, benutzen aber als Ableitungen nur Quotienten 9/o, $/@® usw.; wir rechnen 
auch mit der Koordinate . selbst, erinnern uns aber, daß sie in einem besonderen Maßstab . 
gemessen wird, und zwar so, daß die Einheit den statischen Ausschlag 9, bezeichnet. 


2. Die endliche Bewegungsgleiehung, die Gleichung für die Regelungsfunktion und ihre 
ET | | Ableitung. i 

Er Die Differentialgleichung (1.5b) beschreibt die Bewegung jeweils nur abschnittsweise. 
Sobald die Regelungsfunktion F ihr Vorzeichen wechselt, tritt eine andere Differentialgleichung 
E- auf, da das Störglied in (1.5b) sein Vorzeichen ändert. Auch das Integral der Bewegungsgleichung, 
die Dauergleichung der Bewegung, gilt deshalb jeweils nur für einen Abschnitt zwischen zwei 
Nulldurchgängen der Regelungsfunktion F; es lautet 


PN, = EN (2.1), 


Be et nn ee 


wobei 
h\=-ödö+Ww=ol-D+ifi—D] 

| ed nreof- DI MID re (2.2) 

ist, mit ö= Doundv=@yJ1-—. D. Daß die Zählung der Zeit in jedem Intervall von Null 


an beginnen soll, wurde schon erwähnt. t= 0 bezeichnet den Beginn eines Intervalls; dort 
befindet sich eine Schaltstelle. Die Ableitung von (2.1) lautet: 


En 9p=A,Acktı2,Bet ee A SEN (2. 3). 
Da zur Erzielung einer reellen Lösung die Integrationskonstanten A und B konjugiert komplex 
sein müssen, lassen sie sich schreiben r 
\ A=(el®, B=Üet®; \ 

damit kommt für (2.1) und (2.3) 
= O fe: Era eAat —ie] X (sgn F) 30 
10) 
BA BR a a N RER DEN EN EZ RR (2.3 a). 
Eine weitere, zweckmäßige Fassung der Gl. (2.3a) erhält man dadurch, daß die konjugiert 


komplexen Faktoren A, und A,, die wegen / 6% +»? = w den Betrag » aufweisen, durch diesen 
Betrag ® und das Argument +0 ausgedrückt werden: 


ER. A ET re Var (2.4); 


| 
| 
Ä 


und 


umgekehrt gilt dann 
0) Ne 


cso=————D; sine= = y1-D8 TREE RAR 
10) 


[497 
So kommt an Stelle von (2. 3a) 


= Ofehttichn Heutmieto], 000 2a). 


[097 
In entsprechender Weise läßt sich auch der Ausdruck für die Regelungsfunktion 
F=9-+09 umschreiben. Unter Benutzung von (2.la) und (2.3b) kommt zunächst 
F-( [et + ie (1 + 0 ei°) + ert—ie (1 +0 0 e-'e)] nm (sgn F) > 
0) 
zustande, Drückt man die in den runden Klammern stehenden, konjugiert komplexen Faktoren 
‚ durch ihren Betrag (’ und ihr Argument + r aus, 


r+Wwe#t=(l' er; DHepnwelierie  Nure (2.7), 


_ 


“ \ 


Z.angew. Math. 


BESTEN ENEt Den Fiügge-Lotz und Klotter, Uber Bewegungen einen Schwingere 321 
- so wird 
F=COfektt ern au mie+n] F (sen PL, (2.63); 
dabei ist der Betrag 
DES UISID DR A ER, (2.82) 
und das Argument 
woy1— D: 
== LE RN ER RE Re ? 
T=Aärctg TI Dev | (2.8b) 


Die Ableitung F der Steuerungsfunktion weist vor den e-Funktionen die Faktoren A, und A, 

auf. Sie lassen sich nach (2.4) wieder durch Betrag und Argument ausdrücken. So kommt 
EL odjertietnte BLU N REIHE a ae (2.9). 
@ 


Schließlich geben wir noch an, wie die Integrationskonstanten C und e mit den Anfangs- 


werten 9, und 2 zusammenhängen. Es ist 
U a D. 1° 90 b 90)” 
SR N Posen) Be DE de u 
er = 
1 (0) 
tge= > 
Rn b a 
Fe + (sen MT 
ö, b RE re lan. (2.11). 
role 
sın € E + Po —E (sgn ) 2 9 6) yı RE? D: b) 
bı1 
20 6= (po (sen H) ala e). 
PA) Fl) er; 
Die Gl. (2.1a), (2.3b), (2.6a) und (2.9) zeigen die Funktionen (ti), ae (t) und Eh Wir 


schreiben sie noch einmal an, indem wir die Koeffizienten A, und A, in ihre reellen und ima- 
ginären Teile aufspalten: 


ot) = Met tpei"t +9. er up Fer + (sgn F) S 


? 


ro En Werte Her 6) 4 BL RIER OH 


) 12.12). 
Fi=(0l’endtfeiriretnteTiette tn] TF (sgn F 
F (t) Er (620. endtfeirtts Ho-+7) POTT AMT EA 9] 
M) | 


Die eckigen Klammern lassen sich nachdem Moivreschen Satz auch in reeller Form schreiben: 


ol) =20e?icos(vt+e) 6a), 

PÜ_2Ge=ttcostwi te BE REN ee (213), 
Fiyr=200e cos, wi He +7) Flsgn m 

Et) 


_— —-20(0 etcoswtte+to-tT) 
{63} 


Um spätere Benutzungen vorzubereiten, geben wir noch explizit an, welche Form die dritte 
der Gl. (2.13) annimmt, wenn die Integrationskönstanten durch die Anfangswerte ausgedrückt 
werden: 


b ‚ 
Fti)=+# (sgn Pzerreosri 
1 . 
Tor m BO} Doo + (+ (sgn N) (D—o o)|e- sin vi (sgn N (2.14). 


(02) 


322 Flügge-Lotz und 


‚An den Schaltstellen = 0 wird demnach 


Er 


“ y FR h in » 
’% Ba: Y BP A) 

BAR Ka Kae re E 
ott I, EN NOIR EEE 


h were N 
a, ET Te 
3 Ba s- : SUÜmEEEN 


a) free nante | a RER 


wer [RI L0) | 
wie man aus #=# +09 unter Benutzung der Differentialgleichung (1.5b) erkennt. _ 


© f 


3. Darstellung des Bewegungsablaufs. Die erzeugende Bewegung. Die Phasenebene. 
| Die Ebene der Schaltwerte. 


'Um den Ablauf einer Bewegung darzustellen, kann man sieh verschiedener Hilfsmittel & 


bedienen. Das erste, übliche, besteht darin, die Veränderliche @ als F unktion der Zeit auf- 
zutragen; so erhält man das Ausschlag-Zeit-Diagramm (g-t-Diagramm). Die Ordinate dieses 
Diagramms, der Ausschlag selber, läßt sich dabei für besondere Klassen von Bewegungen, zu 
denen auch die hier interessierenden rechnen, erhalten als Projektion einer „‚erzeugenden Be- 
wegung. Die erzeugende Bewegung einer harmonischen Schwingung ist z. B. eine Kreisbewegung 


"mit konstanter Winkelgeschwindigkeit &; eine Tatsache, die bekannt ist, und von der häufig 


Gebrauch gemacht wird. | 
Die Funktionen g (f), a und F (t) (vgl. die Gl. (2.13)) weisen alle denselben Abklingfaktor 


auf; sie lassen sich also alle auf die Form 


ER aA ELCH De 2% 3 ee 25) Ausee ER (3.1) 


bringen. Solche Funktionen lassen sich aber darstellen als Projektion der Bewegung eines Punktes, 
der mit konstanter Winkelgeschwindigkeit » sich auf einer (um die Strecke y, aus dem Nullpunkt 
verschobenen) Spirale bewegt. Hier ist die Bewegung auf der Spirale die „‚erzeugende‘“ Be- 
wegung. Die Spiralen, die zu den Funktionen p, @/w, F und F/o® führen, sind dabei kongruent. 
Sie unterscheiden sich nur um den Winkel, an dem der Anfangspunkt der Spirale liegt (bzw. 
um den die Spirale gegenüber einer Anfangslage gedreht ist). 

Die genannten Tatsachen sind leicht einzusehen. Nimmt man zu 


te U: (3.2) 
noch die Funktion 


ben nl na. ET Re A (3.2b) 
hinzu und bildet entweder 


z=y-+tie ee r=V + 2, 


so kommt man auf die komplexe oder reelle Form der Gleichung einer logarithmischen Spirale 
in der Form 


= ldriniHte ,...., (3.33) 


ur ER ar 


oder 


Se 
v 


r=(e REN EEE EN 


y(t) ist also die eine, & (t) die andere rechtwinklige Projektion der Bewegung eines Punktes auf 
der Spirale. 


Wir haben im Verlauf unserer Rechnungen von diesem Hilfsmittel der erzeugenden Be- 
wegung auf der Spirale ausgiebig Gebrauch gemacht. Das Auflösen transzendenter Gleichungen, 
die Ausdrücke von der Bauart (3.1) enthalten, z. B. F (t) = 0, gelingt so auf einfache Weise. 

Ein zweites Mittel, einen Bewegungsablauf zu verfolgen, bildet die Phasenkurve (Kurve 
in der-Phasenebene). Koordinaten der Phasenebene sind die Veränderlichen op und 9 oder Viel- 
fache von ihnen; wir werden stets, um Veränderliche gleicher Dimension zu benutzen, p und 9/w 
wählen. Die Zeit £ spielt für solche „‚Phasenkurven‘“ die Rolle eines Parameters; jedem Punkt 
einer Phasenkurve gehört ein fester Wert tzu®). In Formeln: Aus der Gleichung 9 = (t) ent- 
steht durch Ableitung die zweite Gleichung p/o = ? (t)/o; beide Gleichungen zusammen bilden 
die Parameterdarstellung der Phasenkurve. 


°) Phase heißt in der allgemeinen Physik jeder Komplex der einen Zustand bestimmenden Größen (Druck, Tem- 
peratur ...). Ein Bewegungszustand wird vollständig durch die Koordinate und alle ihre zeitlichen Ableitungen bestimmt; 
PP, 9... bestimmen die Phase eines Bewegungszustandes (sie bestimmen einen Punktim Phasenraum). und 9 sind zwei 
der die Phase einer Bewegung bestimmenden Größen. Siesind die Koordinaten der Phasenebene. Ineinem Problem ‚‚zweiter 
Ordnung‘ wird durch zwei Größen, @ und 9, die Phase aber vollständig bestimmt. 
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Für eine harmonisc % eSchwin & ung lautet die Gleichung Ei Ausschlag- -Zeit- 
„Diagramms 


BE NIEREN LINSE, iR (3.4); 
das Diagramm stellt eine Kosinuslinie dar. : Für die Kurve in der, Phasenebene lautet die Para- 
meterdarstellung 


i | ee (ti ta) Wa RR N NN EN (da) 


Fe 4 sin NEO le RE (3.4b). 
Die Gleichung der Kurve, die durch Elimination von t entsteht, lautet 


g? el = 4 } N BR OH A N Ysdoy 


" sie ist ein Kreis vom Halbmesser A. Zu einer harmonischen Schwingung gehört also sowohl als 


erzeugende Bewegung wie als Phasenkurve ein Kreis, Bild 3, 


Wird eine gesteuerte Bewegung der in Ziff.1 und 2 genannten Arti in der Phasenebene 
beschrieben, so liegen die Schaltpunkte für ein bestimmtes System alle auf einer Geraden, der 
„sSchaltgeraden‘‘. Geschaltet wird, wenn #=®-+og‘=0ist; da.dann aber auch die Zählung 
der Zeit neu einsetzt, so istp=@, undo=9% Bo kommt 


Po 
0 
oder 
a a N EN (3.5) 
Po/® 
Die Schaltpunkte liegen also auf der Geraden, deren Na 
Steigung tga& den Wert —wehat. (Istao >0,50 2... .ee-m ZEN 


gehen die Geraden durch den 2. und 4. Quadranten, 
ist®o< 0, gehen sie durch den 1. und 3. Quadranten.) 


(9a; 


IR wg = const: 


Bild 3. Phasenkurve einer Bild 4. Logarithmische Spirale als Phasenkurve 
harmonischen Schwingung. in einem. schiefwinkligen Koordinatensystem. 


. Die Phasenkurven ähneln logarithmischen Spiralen. Es zeigt sich jedoch, daß sie nicht 
genau mit logarithmischen Spiralen übereinstimmen. Wirkliche logarithmische Spiralen erhält 
man als Phasenkurven jedoch dann, wenn man in die Phasenebene nicht ein rechtwinkliges 
(geradliniges) Koordinatensystem, sondern ein schiefwinkliges legt. 

Der Nachweis läßt sich an der Hand von Bild 4 leicht führen: Der Winkel zwischen der 
negativen Abszissenachse (— p/w-Achse) und der positiven Ordinatenachse (9-Achse) sei der 
stumpfe Winkel o, wie er er die Gl. (2.5) als Funktion von D festgelegt ist. N man den 
Kosinussatz auf das schräffierte Be an, so kommt 


2 .\2 bla f 


[02) 
Drückt man hierin 9 und $/® nach den Gin. (2.13) als Funktionen der Zeit aus, so erhält man 
r?—=40?er?dt[cos® (vb te) +cos?(vt He-+o) 
—2cos(wt te)cos(wt te -+o)coseo]. ..... . .. (8.6b). 
Eine Zerlegung von cos [(vt +&)-+-o] nach dem Additionstheorem liefert dann 


ö 
N N HE 


also die Gleichung einer logarithmischen Spirale. 


EN 
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- Die Schaltgeraden, die in einem rechtwinkligen ‚Koordinatensyst wie oben ruhe 

die Steigung tg« = — wo haben, haben in dem genannten schief inklige ee 

' system die Steigung tg« = —tgr. Man erkennt dieses an der Hand von Bild4: Für den 
dort eingezeichneten Winkel a findet man BETRA y % Si 


.y 


un wosno Se RR 
e 5975 0 —p9C08 0 .. 1+m»oc0oso 


= ® E Unter Benutzung der Gln. (2.5) und (2.8b) kommt daraus | En 
RR | dass PER rail AR Pl a ee > 88) 


% 


also 
Hen—T. 


Aus Gl. (3.7) geht auch hervor, daß der Winkel, der am Konvergenzpunkt einer Spirale .. 
zwischen den Strahlen nach irgend zwei Phasenpunkten aufgespannt ist, den Wert »t hat; er 
ist also ein Maß für die Zeit, die die Bewegung in Anspruch nahm. Daraus folgt dann ohne weiteres, 
daß der Winkel-zwischen den Strahlen nach den Schaltpunkten am Anfang und am Ende des 
Intervalls den Wert »t, hat. 

Da für System A der Konvergenzpunkt des Spiralenbogens stets auf der dem Bogen ab- 
gewandten Seite der Schaltgeraden liegt, für System B aber auf der dem Bogen zugewandten 
Seite, folgt sofort, vgl. Bild 5, daß für System A der Winkel vt, <z, für System B aber rt; zn 
ist (s. auch Ziff. 5). 


Dr 
nn 


Dr 


W272) 


Bild 5. Intervall-Längen für System 4 Bild 6. Phasenkurve einer Bewegung aus Stücken von 
und System B, 4 . logarithmischen Spiralen zusammengesetzt. 


Die Tatsache, daß die Phasenkurve im schiefwinkligen System eine logarithmische Spirale 

ist, erleichtert die Verfolgung des Bewegungsablaufs auch durch viele Intervalle sehr. Man 

ö 

zeichnet sich mehrere Windungen einer Spiraler =e ” "N aut. Außerdem zeichnet man in 
ein mit schiefwinkligen Koordinatenlinien g und ®/o, vgl. Bild 6, versehenes durchsichtiges 
Papier die Punkte + b/w® und die Schaltgerade ein. Die Schaltgerade, für deren Punkte die 
Steuerungsfunktion F (t) = ist, teilt dann die Phasenebene in solche 9-#/o-Werte, zu denen 
ein positives F (t) und in solche, zu denen ein negatives F (t) gehört. Der Anfangspunkt der 
Bewegung sei P,, vgl. Bild 6. Es sollentsprechend System 4 gesteuert werden. Die Steuerungs- 
funktiin# = +.» hat im Punkte P, für das gewählte positive o positives Vorzeichen. 
Demnach müssen wir den Punkt — b/o® auf den Konvergenzpunkt der Spirale heften und die 
Spirale so drehen, daß sie durch P, geht. Der erste Schnittpunkt der Spirale mit der Schalt- 
geraden ist dann der Anfangsschaltwert S,. Im Schaltpunkt S, wechselt F das Vorzeichen. 
Wir müssen deshalb beim Weiterkonstruieren der Phasenkurve den Punkt + b/o® an den Kon- 
vergenzpunkt der Spirale heften und können nun von 8, aus bis zum nächsten Schaltpunkt $, 
die Phasenkurve zeichnen. Auf diese Weise kann man den gesamten Bewegungungsablauf 
schnell konstruieren. 

Beim Konstruieren solcher Phasenkurven entdeckt man leicht zwei Besonderheiten, die 
beim Bewegungsablauf vorkommen können. Erstens kann es in der Nähe des Nullpunktes der 
Phasenebene solche Schaltpunkte geben, von denen ausgehend die Phasenkurve die Schalt- 
gerade nicht wieder trifft, Bild 7. Ein solches ‚‚Verfehlen‘‘ der Schältgeraden kann nur im 
Fall B eintreten, im Fall A dagegen nicht, wie stark auch die „Dämpfung‘ sein mag; denn im 
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ie, werden diesen Pall noch. gehenden. nern wir sprechen dabei von ‚sinom „Ruhen Sach de 
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E Bil 7 Phasenkurve einer Rena ungt mit zur , „Ruhe“ 
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2 Au nages Bild 8. ee einer Bewegung. Bild 9. Phasenkurve einer Bewegung mit nachhinkendem 
i x = . Stellglied (konstante Nachhinkzeit » At = 5°). 


Auch die zweite Besonderheit des Bewegungsablaufs sieht man am Phasenbild leicht ein, 
vgl. Bild 8. Wir haben z.B. ein System A vor uns und kommen von einem Schaltpunkt $,, 
dessen Abstand vom Nullpunkt die Größenordnung b/®? hat, zum weiteren S,, ohne die Ylo- 

. Achse geschnitten zu haben, bleiben also links vom Nullpunkt auf der Schaltgeraden. Wenn . 

- wir jetzt in der üblichen Weise weiter konstruieren, finden wir, daß die Phasenkurye auf der- 

selben Seite der Schaltgeraden weiterlaufen will,d.h. die Regelfunktion hat wohl in 8, den Wert 
Null erreicht, hat aber nicht ihr Vorzeichen gewechselt. Da aber der ‚‚Durchgang“ durch Null 
Bedingung für Schalten ist, ist hier das Verhalten des Stellgliedes nicht festgelegt. Nur wenn 
z.B. das Schalten etwas nach dem Nulldurchgang von F einträte (Nachhinken des Stellgliedes), 
könnten wir weiter konstruieren, vgl. Bild 9. Man sieht, daß die Bewegung dann sehr hoch- 

Be; frequent wird und ihre Ausschläge abnehmen. Die Fälle, die auf ein solches undefiniertes: Ver- 

ä "halten der ‚idealen‘ Regelung führen, werden in Ziff. 4 und 7 ausführlich behandelt. 

Im weiteren Verlauf der Untersuchungen wird neben der Phasenebene eine andere Ebene, 
die mit ihr eine enge Verwandtschaft zeigt, häufig benutzt werden. In der Phasenebene mit den 
Koordinaten » und ®/w läßt sich die Bewegung noch lückenlos beschreiben ;'eine Phasenkurve 
bildet die Aufeinanderfolge aller Zustandspunkte. Wir werden uns bei unseren weiteren Unter- 
suchungen aber oft damit begnügen, statt der Aufeinanderfolge aller Zustandspunkte die Reihen- 
folgeder ,‚Schaltpunkte‘‘ zu betrachten,d.h. derjenigen Werte 9, und $/w=®,/w, die im Augen- 
blick des Nulldurchgangs der Funktion F vorhanden sind. Aus den Punkten der Phasenebene 
werden also die (auf einer Geraden, der Schaltgeraden, liegenden) Schaltpunkte @,, 9,/® heraus- 
genommen und zur Kennzeichnung der Bewegung verwendet. Fügt man aus den zu allen mög- 
lichen Systemen gehörenden Schaltgeraden wieder eine Ebene zusammen, so erhält man die pg- 
o ,/@- De oder Ebene der Schaltwerte. In dieser Ebene der Schaltpunkte wird ‚die Bewegung 


"eines Seins durch Si eis solle BR n. 


einem entweder im Punkte + b/o® oder — b/o* liegenden Kon 


wegung zwischen den Schaltpunkten erhalten wir keine Auskı nftı 
male der Bewegung, vor allem das Abklingen, Kehl ao lasse 
einer solchen Darstellung erkennen. Denn die Phasenkurven sind logz 


vergenzpunkt. Man 

daß die Bewegung zwischen zwei Schaltpunkten keine Besonderheiten Br s Abkl 

der Bewegung. ist gesichert, wenn die Folge der Schaltpunkte nach Null ART 

Während wir.aber in die Phasenebene ein schiefwinkliges Achsenkreuz gelegt ‚werden | 

wir in die Schaltebene stets wieder ein rechtwinkliges legen. Die Schaltgeraden erhalten dann ER ei; 

wieder die Steigungen tg« = —w o. Von dieser Ebene der Schaltpunkte, der Yo-Pol@-Ebene, Per % | 

wird weiterhin viel die Rede sein. | 
Zur Erläuterung zeigen wir ein Beispiel: Von der Bewegung eines Systems, die durch die n 

Phasenkurven, Bild 10 (hier in kartesischen Koordinaten gezeichnet), wiedergegeben wird, = 

erscheinen in der 9,-9,/w-Ebene nur die Schaltpunkte, Bild Hin Für das vorgegebene System 


A=-arc tgl-wg) 


Bild 10. Phasenkurve in kartesischen Bild 11. Schaltpunkte. Bild 12. Schaltpunkte einer. 
‚ Koordinaten (schematisch). periodischen Bewegung. 


mit festen Werten o und wliegen die Schaltpunkte auf der Geraden mit der Steigungiga = 
— oo. Die Punkte auf anders geneigten Geraden gehören zu anderen Systemen. Man muß 
jedoch beachten, daß die Differenz der Entfernungen zweier Schaltpunkte vom Nullpunkt noch | 
kein Maß ist für die Abklinggeschwindigke it. Diese Differenz mißt das Abklingen i 
je Intervall. Um zu einer Abklinggeschwindigkeit zu kommen, muß noch die Intervall-Länge f 
berücksichtigt werden, vgl. Ziff. 6, S. 12. 

Eine periodische Bewegung wird z. B. durch zwei Schaltpunkte wiedergegeben, die nach- 
einander immer wieder durchlaufen werden, Bikl 12. Diese Punkte brauchen aber keineswegs 
den Maximalausschlag zu bezeichnen, der während der Bewegung auftritt; Maximalausschläge 
sind vielmehr nur dann Schaltpunkte, wenn 9,/® = 0 ist; solche Punkte liegen auf der Or- 
dinatenachse. 


4. Aufteilung der Ebene der Schaltpunkte in „Durehgangspunkte“ und „Randpunkte“ 
(„Anfangspunkte‘“ bzw. ‚„‚Endpunkte“) und in Schaltpunkte a bzw. Schaltpunkte b. 


Wir gehen aus von der Betrachtung des Zeitablaufs der Funktionen F (t) und m (t), wie 
erin Bild 2 dargestellt ist. Die Nullstellen der Regelungsfunktion F (£) bezeichnen die Schalt- 
stellen. An diesen ihren Nullstellen weist die Funktion #=g9 +» o(P/o) stets einen Knick 
auf, denn po weist einen Knick auf, weil 9/&® — wie aus der Differentialgleichung, etwa (1.5b), 
folgt — einen Sprung macht. Der Sprung von 9@/w* an einer Schaltstelle ist 


2b 
a (+0) — a (0) 7 für Schaltstellen a, 


R ” A Er (4.1a) 
| (+0) RE (—0)= für Schaltstellen 5. 
Wegen der aus der Definition von F folgenden Beziehung 
Fo $ 
re 


ergibt sich daher für den Sprung von Flo, der den Knick der Funktion # kennzeichnet, an den 
Schaltstellen a der Wert 


F A] b 
HIT Vega 
an den Schaltstellen 5 En a SR N N > A (4.T’b). 
F Be b 
4959) =—- 2005 


> ‚sind. Ein er Wert, der ee ee en SR 
IL größe anal bei fallender Funktion F- dagegen ein 
 „Brechen hin zum Lot‘, ‚Bild 13b, bei steigender Funk- 
tion F ein „Brechen weg vom Lot“. „Brechen hin zum N is Be 
‚SEE vergrößert den Betrag der Neigung, „Brechen RN 
TR weg ‚vom Lot‘ verkleinert ihn. Bi ” 
Ar aka " Beachtet man nun, was in Gl. (4.1b) uhr das een der Sprunggsße Be ist, 
"zusammen mit dem in Ziffer 1 über das Auftreten von’Schaltstellen « und bi in den nteruen A k, 
ER und B Be so. ER sich er, Möglichkeiten: BR DE 1 a 3 A 23 ” 


| Mater a. VER LETN 3 
| OS Ba NT ee 
"o0> 0. |  Brechen’ weg vom Lot | Brocken. hi Tales Ba, VL: 


@< 2 Brechen Bm zum Lot  Brechen NR vom Lot 


Bu 


Sure - Die Zuordnungen sind unabhängig vom „Charakter“ der Schaltstellen, oba BR HAAR 
= la a ‚Beispiel erläutern wir den ersten Fall: Ist ® 3 >: 0, so sind nach (4. 1b) die Bine 7 


weg vom it eintritt. In Weise kann man die aa drei Fälle erörtern. 
Aus.der Tatsache, daß der Knick (Neigungssprung), den die Funktion F an. den Schalt- 
stellen aufweist. einen endlichen Betrag hat, folgt nun, daß es in jenen Fällen (d i. Quadranten), a 
in denen ein Brechen ‚‚weg vom Lot“ auftritt, ‚‚Einfallsrichtungen‘ gibt, denen keine ‚„‚Ausfalls- ‘B 
- riehtungen‘“ auf der anderen Seite der Zeitachse mehr zugehören, so daß also kein Nulldurchgang - PER 
der Funktion F, also auch kein Schalten mehr zustande kommt. Diese Fälle entsprechen der Br 
Totalreflexion des optischen Analogons (schraffierter Winkelraum in. Bild 13 a).  Um- 7 ee 
| . gekehrt gibt es beim Brechen „hin zum Lot“ „Ausfallsrichtungen‘, denen keine „Einfalls- IE 
cK richtungen‘ zugehören, Bild 13b. In beiden Fällen tritt also kein Schalten mehr auf: die Null- 
“ "punkte (Schaltpunkte) der Funktion F sind keine ‚‚Durchgangspunkte‘‘ (der gesteuerten Be- 
wegung) mehr, sondern ‚„‚Randpunkte‘“, bei denen die Bewegung entweder enden muß (,‚End- 
"punkte‘‘) oder beginnen kann (,‚Anfangspunkte‘‘). Sowohl die Einfallsrichtungen ohne Ausfalls- 
richtungen, wie die Ausfallsrichtungen ohne Einfallsrichtungen liegen in den in Bild 13a und b 
schraffierten Winkelräumen. Die Begrenzung des Winkelraums wird jeweils durch die Bedingung 


geliefert, daß für die einfallende oder die ausfallende ruhe F = 0 sei. Setzt man für F den 
Ausdruck nach (2.15) ein, so kommt 


p b | | a 
() +2D Po (san F) a Pu +0, =. GE Lea Le Te (4.2). Y 2 


Wenn D< List, stellt (4.2) sowohl für System A (oberes Vorzeichen) wie auch für System B 
(unteres Vorzeichen) in der Schaltebene zwei Ellipsen dar, da sgn F sowohl +1 als auch — 1 
sein kann, Bild 14. Hat das dritte Glied negatives Vorzeichen, so beschreibt die Gleichung die 
obere Ellipse, hat es positives Vorzeichen, so beschreibt sie die untere Ellipse. Auf dem Rand der 
Ellipsen liegen die ‚‚Grenzfälle der Totalreflexion‘, innerhalb der Ellipsen liegen jeweils Rand- 
punkte, außen Durchgangspunkte. 
Ob die im Innern der Ellipse liegenden Randpunkte nun Endpunkte oder Anfangspunkte 
“sind, kann man etwa in folgender Weise entscheiden: Für Endpunkte tritt stets Brechen weg 
vom Lot ein, für Anfangspunkte Brechen hin zum Lot. Die Tafel 4/1, die die Quadranten in 
dieser Hinsicht einteilt, gibt also Auskunft auch darüber, wo Endpunkte und wo Anfangspunkte 
im Innern der Ellipsen liegen. In Bild 14a und 14b sind die Ergebnisse dieser Betrachtung ein- 
getragen. Die durch Gl. (4.2) dargestellten Ellipsen nennen wir zur Unterscheidung von anderen, 
später noch zu erörternden Kurven gelegentlich die ‚‚kleinen Grenzellipsen“. 


Wie in Ziff.1 auseinandergesetzt wurde, unterscheiden sich die Schaltpunkte a und b 


E (0) 
& 


ist in 


nach dem Vorzeichen, das die Funktion F an der Stelle O0 aufweist. Die Funktion 


(2.15) angegeben. F (0)/o ist also positiv, solange beide Faktoren, @,/® und der Klammer- 


Vorzeichen zeigen. Der Grenzfall P/o = 
schwindet; die Bedeutung dieses Falles haben wir 


stellt die oben erörterten ‚‚kleinen Gren ‚ellipsen 
. ‚der Ellipsen ist X < 0, für Schaltwerte aus dem 


Be RE 


N | a "Bild 15. Gebiete der Schaltpunkte@ undd. 
Bild 14. Kleine Grenzellipsen für System 4 _ Bild 15. Gebiete der a 
und System B mit Anfangs- und Endwerten. 7 (Linke) System A. (Becute) TER R BER 


00 Wir fragen nun zunächst: Welche’ Gebiete der 94-%,/o-Ebene sind für das System A von 

BENY- Schaltpunkten aerfüllt ? Zu System A gehört das obere Vorzeichen des dritten Gliedes in X (2.15). 

x Fürjeden Schaltpunkt a gilt (vgl. Bild 2a) im System A ferner, daß (nach dem Schalten) sgn FERN 
0 —L ist. Damit nimmt das in Rede stehende Glied das negative Vorzeichen an; von den beiden EHE 

„„kleinen Grenzellipsen‘‘ kommt also nur die obere in Betracht. Schaltpunkte a im System A % 2 
erfordern außerdem F < 0, also verschiedenes Vorzeichen der Faktoren Polo und K. In der % nr 

. rechten Halbebene, wo 9,/® > 0 ist, muß K < 0 sein; die Schaltpunkte a liegen also im Innen- 

raum der (oberen) Ellipse. In der linken Halbebene, wo @,/® >0 ist, muß K >0 sein; die 
Schaltpunkte a liegen im Außenraum derselben Ellipse, vgl. Bild 15a. RN. 

Untersucht man Schaltpunkte 5, so wirdsen = + 1, das dritte Glied in K wird positiv, EEE 
es kommt die untere der Grenzellipsen in Betracht. Da jetzt F/» > O sein muß, liegen die Schalt- 
punkte bin der rechten Halbebene im Außenraum, in der linken Halbebene im Innenraum der 

(unteren) Ellipse, Bild 15b. RER N SERaea 
Die Untersuchung für System B verläuft ganz entsprechend. Hier gilt für das dritte Glied 
in K (2.15) zunächst das untere Zeichen. Für Schaltpunkte « ist hier sgn #= +1, so daß das 
Glied das negative Zeichen erhält; damit handelt es sich wieder (wie bei den Schaltpunkten a 
‘ des Systems A) um die obere der beiden Ellipsen. Die Schaltpunkte « erfordern nun F/o >, 
sie liegen also in der rechten Halbebene im Außenraum, in der linken Halbebene im Innenraum 

. der (oberen) Ellipse (Bild 15e). Für Schaltpunkte b kommt wegen sgn #= —. 1 die untere 

Be Ellipse in Betracht; wegen F/o< 0 liegen die Schaltpunkte 5 in der rechten Halbebene im 

0 Innenraum, in der linken Halbebene im Außenraum der (unteren) Ellipse, Bild 1öd. 

Denkt man sich Bild 15a und 15b von System A übereinandergelegt, so sieht man, daß 
innerhalb der Ellipsen links oben und rechts unten überhaupt keine Schaltpunkte a oder 'b 
liegen. Diese Punkte erkannten wir in Ziff.3als Endpunkte der Bewegung. Die Anfangs- 

‚ punkte, die innerhalb der Ellipsen rechts oben und links unten liegen, werden von der Betrach- 
tung dagegen mit erfaßt, da die Definition der Schaltpunkte von dem Vorzeichen der Funktion # 
nach dem Nulldurchgung Gebrauch macht, und dieses Verzeichen auch für Anfangspunkte 
tatsächlich erklärt ist. Ja, es zeigt sich, daß diese zuletzt genannten Gebiete beim Übereinander- - 
legen von Bild 15a und 15b sogar doppelt überdeckt sind, also sowohl Punkte vom Typ a wie 
solche vom Typ benthalten. Die Bewegung kann also von allen Anfangspunkten aus in zweierlei 


Art beginnen. Außerhalb der Ellipsen liegen Durchgangspunkte, die nun eindeutig in Schalt- 
punkte « und 5 aufgeteilt sind. | 2 


5. Kurven konstanter Intervall-Länge »t, in der Ebene der Sehaltwerte 
(Isoehronen) und ihre Sonderfälle. 


Die Länge der Intervalle vi, ist für die einzelnen Sehaltpunkte verschieden. Wir fragen 
hier nun nach den geometrischen Örtern aller der Schaltpunkte, für die das nach dem Schalten 


sich anschließende Intervall gleich groß ist. 'Es läßt sich zeigen, daß diese geometrische Örter 
Ellipsen sind. / 


Aus der Bedingung, daß 


u 
” « 


Fi) 


ng. ner Ellipse i in der Schaltebene mit vtyals RE 35 
zeigt, daß für das System A nur solche Ellipsen auftreten 

En B nur solche, für die vt, Zr ist. Diese Tatsache ist 
‚daß das Regelmoment das Rückstellmoment im System 


System. Bi im wesentlichen schwächt. uch ein ‚Blick auf die 


g heat im Fall Brergiäliert wird. * 
Ziff. 3 hatten wir auch schon erkannt, daß bei RE Fe Tilmann > R 
die P asenkurven ‚die ‚Schaltgerade schließlich nicht "mehr schneiden, so daß kein weiteres 
ten mehr eintritt. Den letzten 'Schaltpunkt nennen wir. „Ruhepunkt‘. Im ee 4 % 
t-Diagramm sieht. , der\.. 5 Een LA z 

enzfall einer solchen B- 

gung so aus, wie Bild 16. 
ngibt.. Dieser Grenzfall ist 


ZU 5.1) hinzutretende Be ; 
Ru, BER HERE 
on gen: 263. 
Reel. (5.3) wählt 
' nun aus der Schar der Ellip- 
sen (9.2) eine einzige aus und 
‚legt den Parameter t, fest. 


=" aa 16. Branztol, des „Ruhens“. _ Bild 17. Kleine und große aa für verschiedene Werte D. 


Aus den Bedingungen (5. 1) und (5. 3) ER: sich unter Benutzung der Gl. (2.14) die trans- 
 zendente: Beziehung 
h 0 A j —.rt, TR 


6. z > 
OL PUT SE NEE TE AEION) Be: 


Br; . 32 Y ? Wr; 


für den Parameter vt,. Sie ist erfüllt 
1. jüurst, = 
2. für einen von Null verschiedenen Wert vi, — »t.. 
Dieser letztere ist der, der dem ‚‚Grenzfall des Ruhens‘“ entspricht. In der Schaltebene wird 
also durch »t, = vt, eine Ellipse festgelegt; sie heißt gelegentlich die ‚große Grenzellipse“. 
Die aus Gl. (5.4) folgende größte Intervall-Länge hat für D = 0,1,den Wert 


vt, = 300°. 


Solche Hehsrnikka können nur im System B auftreten. 


Übrigens ist leicht zu sehen, daß auch die in Ziff. 4 erwähnten ‚‚kleinen Grenzellipsen‘“ 
in de Schar der Ellipsen konstanter Intervall-Länge gehören. Und zwar hat ihr Parameter 
den Wert vt, = 0; denn man erkennt sofort, daß aus Gl. (5.2), wenn man dort den Grenzübergang 
yt,—0 ausführt, die Gl. (4.2) entsteht. 

Wie wir eingangs erwähnt haben, legen wir in allen numerischen Rechnungen dem 
Dämpfungsmaß den Wert D= 0,1 bei. Demgemäß sind auch alle Bilder hergestellt. In Bild 17a 
und 17b zeigen wir nun einmal (für Schaltpunkte a), wie die kleinen und großen Grenzellipsen 
aussehen, wenn D andere Werte annimmt. Gezeichnet sind die Kurven für D=0, D=0,1 
und D= 0,5. In jedem Fall gehen die Kurven durch die Punkte O und +1. Für D= 0 wird 


SON BEIM 5 : ei f 2 ER Dangı 
BR Über PM RENTE EN 
die kleine und auch die große Grenzellipse zu dem symmetrisch zur go-Achs 


‚ durch die Punkte O0 und +1. SE nr aa ; | | 
Legt man— wie wir das taten — in die Ebene der Schaltwerte ein kartesisches Koordina 


hi Ar 
ter, 


statt des rechtwinkligen Koordinatensystems das schiefwinklige verwendet, das in Ziff, 3 für 


Winkel otrittin Erscheinung, wenn man ihn nach (2.5) an Stelle des Dämpfungsmaßes Deinführt. 
Be ü 


6. Energiebetrachtungen und periodische Lösungen. 


Die Schwarz-Weiß-Regelung hat für die vorliegende Aufgabe den Zweck, ein wesentlich 
stärkeres Abklingen der Bewegung zustande zu bringen, als es mit der vorhandenen Dämpfungs- 
kraft allein eintritt. Es sind in einem solchen geregelten System jedoch nicht nur abklingende, so- 
dern auch periodische und aufschaukelnde Bewegungen möglich, da das Regelmoment dem 
System nicht nur Energie entziehen, sondern ihm auch Energie liefern kann. 

Um uns einen Überblick zu verschaffen, wo in der 94-9,/®-Ebene die Schaltwerte für 
abklingende, wo für periodische und aufschaukelnde Bewegungen liegen, fragen wir zunächst 
nach jenen Punkten in der Ebene der Schaltwerte, von denen solche Bewegungen ausgehen, 
denen das Regelmoment Energie entzieht, und solche, denen ‘es Energie liefert. 


Die Energie, die das Regelmoment m dem System liefert, ist gegeben durch den Ausdruck 


U, = / DE N N (6.1) 

Mit N ea en (6.2), 
wobei wieder das obere Zeichen für Systeme A, das untere für Systeme B gilt, wird daraus 
V,.=rFlsgn FM (lo, =). ee (6.3a) 


Nun untersuchen wir, welches Vorzeichen dieser Ausdruck für Schaltpunkte in den einzel- . 


nen Quadranten hat.- Wir betrachten zunächst System A. In Ziff. 3 stellten wir fest, daß für 
System A Schaltpunkte a im 2. und 3. Quadranten, Schaltpunkte 5 im 1. und 4. Quadranten 
liegen. Nach Schaltpunkten aim System A ist san # = — 1, nach Schaltpunkten 5 ist sgan f = 
+1. Unter der weiteren Voraussetzung, daß die Bewegung stets so verläuft, daß p, und p, auf 
verschiedenen Seiten des Nullpunktes liegen ®), wird der Faktor (9, —g,) für Bewegungen, die 

von Schaltpunkten im 1. und 2. Quadranten aus- 


Tafel 6/1. sgnU, gehen, negativ, für solche, ‘die von Schaltpunkten 
‚Bewegung geht { im 3. und 4. Quadranten ausgehen, positiv. 
aus von einem | System A | System B Unter Beachtung dieser Tatsachen findet man 
Sehaltpunkt im die erste Spalte der nebenstehenden Tafel. 
1. Quadranten 1 et, Für das System B dreht sich zunächst das 
2. = N rt erste Vorzeichen in (6.3a) um. Ferner vertauschen 
“ „ TR a sich die Zuordnungen von sgnF zu den Schalt- 


punkten a und b; da aber zugleich die Schaltpunkte « 
und 5 ihre Lage in den Quadranten ändern, bleibt 
die Zuordnung von (sgn F) zu den. Quadranten erhalten. Ebenso bleibt die Zuordnung 
von sgn (%ı — 9.) zu den Quadranten erhalten. Das heißt, das Vorzeichen von U, ist im 
System B für alle Quadranten dem im System A entgegengesetzt (Spalte 2 der Tafel 6/1). 


Wo sgn U,= +1 ist, wird durch das Regelmoment dem System Energie zugeführt, 
wo sgn U,—= — 1 ist, wird Energie entzogen. Wo aber außer durch das Dämpfungsmoment 
Energie auch durch das Regelmoment entzogen wird, wird die geregelte Bewegung noch stärker 
abklingen als die ungeregelte. Wo dagegen Energie zugeführt wird, Klingt sie weniger ab, ja sie 
kann periodisch werden oder sich aufschaukeln. 


Die Bewegungen, die von Schaltpunkten in jenen Quädranten ausgehen, in denen das Regel- | 


moment Energie zuführt, müssen nun noch genauer betrachtet werden, um zu unterscheiden, 
wo die Schaltpunkte liegen, von denen periodische, abklingende oder aufschaukelnde Bewegungen 
ausgehen. 
LET REN 
‘) Für Intervalle, die zu Endpunkten führen oder von Anfang spunkten ausgehen, kann es vorkommen, daß 9, und 9, 


gleiches Vor zeichen haben. Man muß dann eine ge sonderte Betrachtun ans T edoch ebenfalls zu dem in 
Tafel 6/1 ingetr igenen Ergebni 
g tellen; sie führt j 


: | 
system, so sind die Kurven konstanter Intervall-Länge t, Ellipsen. Ohne die Transformation # 
explizit durchzuführen, erwähnen wir nur, daß alle diese Ellipsen zu Kreisen werden, wenn man e A 
die Phasenebenen benutzt wurde (bei dem die negative $/o-Achse mit der positiven g-Achse Be: 
den Winkel o einschließt). Die Gl. (5.2) ist dann die Gleichung der Schar dieser Kreise; der 


ten- 


u u a 


‚as os dat Moe er © u]. f 6.4 . 
Er Drü man nun die Funktionen rückwärts wieder durch [2 ‚9olo.: pı und dılo aus, so erhält ı man a 


“ 


ns en 1 I ern. = D (e& 


6. 8). 


er eu Ds) ae - (sgn un, = ( Kan + (9 De 
b x er % Wir Re nun ee ‚zeigen, daß ” v für solche Schaltpunkte, ‚die zu periodischen. 
BEN Bewegungen ‚gehören, verschwindet. 

Die periodischen Bewegungen selbst können ohne eine- Energiebetrachtung. tunen 
werden: Für periodische Bewegungen bestehen zwischen den Schaltwerten 9, 9/® am Anfang 


a Intervalls und den Schaltwerten 91, Yılw am Ende des Intervalls die Beziehungen ’ 


“ P=—Ppı PEN UNE . (6.7) 
N ea 


er Periode ist dann die dee Intervall- -Länge, FT = 2b: Mit Hilfe der Gin. (2.13a) und (2.13 En 
| ‚lassen sie sich auf die Formen 


a u Pak: 
a FR, 200876), a7 20eos(e Hr) tbsenhZ. Napa (6.8) 
So r Bi cos(e+)=—cos(ie +0 +vi,) FR RER . (6.8b) 
GR bringen. Die Gln. (6. 8a) und (6.8b) stellen zwei Beziehungen zwischen den RER. : 
j ten und e sowie der he co einerseits und der Intervall-Länge »t, andererseits dar. RN; 


t Drückt man die Winkel e und o mit Hilfe von (2.11) und (2.5) durch die Schaltwerte o, und 
0. Polw aus, so kommen (nach längerer Rechnung) die Gleichungen 


RR er. sin (/1— Dot) — &inDwot, 
7 cos (/1— Dot.) +Coj Dot, 


« ; 2 /T)2 
Ro ion F) a De REIT AT DER (6.96) 
® o*/ 1 — D2cos (1 — Dot,) +&0j Dot, 


Ar Er areleltung der. ‚Kurve der periodischen Lösungen“ in der Ebene der Schaltwerte 
zustande. 


Diese Kurve ist eine spiralenähnliche Kurve (vgl. Bild 28a und b), die Intervall-Länge t, 

ist Parameter. Die Kurve beginnt für t, = 0 im Koordinatennullpunkt und läuft für 1,— © 

in den Punkt, =+1,9,/®=0 ein. Für vi, = schneidet die Kurve die Ordinatenachse. 

Die Punkte vt, < x gehören (nach dem früher Gesagten) zum System A, die Punkte »i,> x 

zum System B. Überdies hat die Kurve ein Ende dort, wo vt, = vi, wird. Sie endet also 

auf der „großen Grenzellipse‘ (die das Gebiet der Ruhepunkte begrenzt); und zwar. läßt 
‚sich zeigen, daß a Punkt zugleieh Schnittpunkt mit der kleinen Grenzellipse ist. 


(6.9 a) 
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Setzt man nun in (6.6), nachdem man auch ( noch durch g 
die Schaltwerte nach (6.9) ein und berücksichtigt, daß überdies 9, = 
ist, so findet man ven 


0 


oder 


PATE 
D=-U, N 


Wie stets, so ist auch für die Kurven von Bild 27 und 28 der Wert D = 0,1 des Dämpfungs- EN 
" maßes zugrunde gelegt. Für andere Werte des Dämpfungsmaßes sehen die Kurven der perio- 
dischen Bewegungen so aus wie Bild 18 angibt, das sich auf Schaltpunkte «a bezieht. FürD=0 


artet die Kurve zu einer Geraden aus, die mit der 9,/w-Achse zusammenfällt. u 


Weiterhin interessiert nun das Vor- Re 


| B zeichen, das die Energiedifferenz A U N 
ystem | zu System B halb und innerhalb der Kurve der perio: 
= zu a 6 1. 2 5 4,5 Schen Bewegungen (6.9) annimmt. 
De er a BE Re Eine allgemein gültige Diskussion des 


Vorzeichens an Hand der Gl. (6.6) scheint 
schwierig zu sein. Eine große Zahl von 
numerischen Beispielen zeigt jedoch, daß für 
System A im 1. und 3. Quadranten und für 
System B in dem Sektor zwischen der 99 
Achse und dem Strahl O8Q, Bild 28, die Be- 
wegungen stets gegen die periodischen kon- 
vergieren, und zwar von Schaltpunkten 
außerhalb der Kurve der periodischen Be- 
wegungen mit abnehmenden Ausschlägen, von 
Schaltpunkten innerhalb jener Kurve mit 
zunehmenden Ausschlägen. Im Sektor zwi- 
schen den Strahlen O8SQ und OR gehen die 


Bild 18. Schaltwerte a für periodische Bewegungen bei Bewegungen von „inneren“ Punkten mit _ 


verschiedenen Dämpfungsmaßen D. wachsenden Ausschlägen gegen die periodi- 
schen, Bild 29. 

Der Bogen RS der Kurve der periodischen Bewegungen ist „instabil“, d.h. von Schalt- 
punkten in seiner Nähe geht die Bewegung, falls der Punkt „innen“ liegt, gegen eine periodische 
Bewegung auf’ dem weiter außen liegenden Bogenstück, falls er „außen“ liegt, klingt die Be- 
wegung in der üblichen Weise ab (und zwar führt sie zu einem Ruhepunkt). 


Eine gute Übersicht über den eben besprochenen Verlauf der Bewegung und damit über‘ 


die Frage des Entzuges oder der Zufuhr von Energie erhält man aus Bild 19 bis 24. In Bild 20,21 
und 24 ist über der dimensionslösen Intervall-Länge »t, der zugehörige Schaltwert @, am An- 
fang und 9, am Ende des Intervalls für verschiedene Werte o aufgetragen. Wählt man einen 
ersten Schaltpunkt, z.B. 8, in Bild 21, so sind alle folgenden durch Einzeichnen eines Treppen- 
zuges leicht zu ermitteln. Vom Schaltpunkt S, aus erhält man keinen Schnittpunkt mehr mit der 
Kurve |p,|. Der Punkt S, entspricht einem Punkt im Gebiet der kleinen Grenzellipse. Für 
e= (0 muß man, weil alle, = 0 sind (n=0,1,2...), 9,/® auftragen, vgl. Bild 19. Für 
e—+ @ist vt, stets gleich z und außerdem @, stets gleich Null. Man muß also auch hier eine 
andere Auftragung wählen. Wir haben die von K.Bögel°) benutzte gewählt, denn für 0— © 
gehtja unser Problem in das von Bögel behandelte über. Für diesen Fall oE—+ % läßt sich 
eine einfache Beziehung zwischen den Schaltwerten p, und g, angeben, s. Bild 22 und 23. 


Noch eine weitere Frage tritt im Zusammenhang mit den Energieüberlegungen auf. Schon. 


die freie Schwingung eines gedämpften Systems verläuft mit Energieentzug. Soll die Regelung 
des Systems die Bewegung rascher zum Abklingen bringen, dann muß die Abnahme der Maximal- 


ausschläge stärker sein als bei der freien Schwingung. Für die freie Schwingung ist das Verhältnis 
D 


a JENA 
zweier aufeinanderfolgender Maximalausschläge e et . Für die gesteuerte Schwingung,.ist 
das Verhältnis nicht konstant und auch die Zeit tyr zwischen zwei Maximalausschlägen ist nicht 
konstant. Wie man aber vermutet, liegt Zj, stets zwischen dem zugeordneten Intervall t,1 und 
ku V8l. Bild 25a. Beweisen kann man diese Beziehung mit Hilfe des Phasendiagramms, 
Bild 25b. Zu diesem Zweck ist vom Schaltpunkt $, aus nicht wie üblich weiter konstruiert, 
sondern da wir uns nur für die Absolutwerte der maximalen p-Werte interessieren, S, am Null- 


°) K. Bögel, Das Verhalten gedämpfter und aufschaukelnder freier Schwi a i 
kung, einer konstanten Reibungskraft. Ing.-Arch. Bd.12, S. 247, es eu 
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Wir können nunmehr das Verhältnis zweier aufeinanderfolgender Maximalwerte, 


DER-==uRE 
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Durch die Konstruktion ist eine einfache Zuordnung von iy zu den benachbarten Interva 
S.5 


Längen gegeben. Deshalb ist für System A auch vty <r und für System 


vgl. 


( 


a e. PER UN TE 


Bild 25. Zur Konstruktion von iy und gmaxn. 


berücksichtigen, daß im geregelten. Systeht die Anderung des Maximalausschlags i in ENT 
Zeiten erfolgt. Die Abklinggeschwindigkeit im geregelten System ist (es wurde @maxn der 
ungeregelten Bewegung mit un, Ymaxn der geregelten Bewegung mit Pgn bezeichnet): 
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Bild 26. Verhältnis zweier aufeinanderfolgender Maximalausschläge, aufgetragen über ihrem zeitlichen Abstand. 
Im ungeregelten System, d.h. bei einer freien gedämpften Schwingung ist 
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Nur wenn das Verhältnis. velor grüße als 1 ist, erfüllt die ee | 
deshalb das Verhältnis vg/vy bei gleichem PRSSARREENT Pan = 


Den Zähler 7@r+!__ 1 können wir aus Bild 26 sofort entnehmen ®). Der Nenner 
Pun D St 
"MER, _ > BEE | v vim 

7 

stellt eine Gerade dar, die zum Vergleich eingetragen ist. Sie geht durch den Nullpunkt End durch 
- den Punkt $. Für alle Punkte, die unterhalb der eingezeichneten Geraden $ liegen, ist das Verhältnis 

valun >1,d.h.diegeregelte Bewegung klingt schneller abals die ungeregelte. Dasbedeutetaber, daß _ 

in der Regel das System A vorzuziehen sein wird. Auf diese Frage kommen wir in Ziff.7noch zurück. 


e 


E 7. Zusammenfassende Diskussion der möglichen Bewegungsabläufe 

Ir r (Diskussion Bild 27 und 28). : 

Die in den vorangegangenen Ziffern gewonnenen Erkenntnisse sind in Bild 27 und 28 

3: ' zusammengefaßt. In beiden Bildern deutet die Strichart die Art der Schaltpunkte an: Aus- 

gezogene Kurven gehören zu den Schaltpunkten a, gestrichelte zu Schaltpunkten b. Bild 27a 

er und 27b zeigen jeweils für System A und System B die Kurven konstanter Intervall-Länge ti, 
(Isochronen). Man erkennt dabei die in Ziff. 5 erwähnte Tatsache, daß im System A der Para- 


Zi De a Ei BE Zn 


System A 


E m on\ ya Zndp. 
Kautenen A hier \ giesseisd 


Null-P 
“Von hier } 
zu fngpunkten N 


\venseits des 
Mulpunktes Fi 


Nach zweı Intervallen 
werden Endpunkfe erreicht 
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Bild 28a und b. Aufteilung der Ebene der Schaltwerte. 
°) Zur Grenze des Ruhens in Bild 26 ist folgendes zu bemerken: Auf S.10 ist fü 
r den zeitlichen A i 
aufeinanderfolgender Maximalwerte im Bereiche normalen Wirkens der Regelung betrachtet, z. B. a % a 
? 


Skizze Bild 26 oben rechts. Wenn aber im System Bi 
L »ystem B im Verlauf einer Bewegung Ruhen auftritt, so ist m i i 
der Betrachtung des Abklingens des geregelten Systems als letztes 'Maximalwertepaar Pq, und O7, zu tere N elle 


Abstand (viy)r dieser Maximalwerte kann ab 3 
N B er größer sein als vb, wie man leicht aus der Phasenkurve entnehmen kann: 


‚Bda, char dieser Kurven konstanter il % 
g An N) und = In die in Ziff. 4 und 5 als die „kleinen“ und 
% großen Gre psen bezeichnet worden waren, und deren Eigenschaften dort erörtert 
Ei wurden. 'In Bild 29 sind für das System. Bund Sehaltpunkte a die kleine und die große Grenz- 
De ellipse noch einmal herausgezeichnet, undes ist die sich aus der Betrachtung i in Zi. 4 ergebende, 
Bedeutung der einzelnen Bereiche kenntlich gemacht. , 

In Bild 28a und 28b sind außer den kleinen Grenzen die das Gebiet: der Randpunkte H 
umschließen, noch j jene Gebiete kenntlich gemacht, von denen aus man nach einem oder 2 
. zwei Intervallen zu einem Endpunkt gelangt. RER | 
 DasGebiet,von demauseineinziges Interyallzu _ Anfongopunkte, 
. einem Endpunkt führt, ist noch einmal unter- Fangspunkte. "2 
teilt, und esist angemerkt, obdieser Endpunkt  . 
Auer, diesseits oder jenseits des Nullpunktes liegt. 
Für System B sind neben Endpunkten : 


Puegangsunkte Rt: f 
Ana Dp.) 


\ ‚Ruhepunkte Bat 


a auch Ruhepunkte der Bewegung möglich, RER, 0 
_  Schaltpunkte in den halbmondförmigen Ge- ink) 2 


punkten) 
„bieten führen nach einem oder nach zwei DE 
' Schritten zu solchen Ruhepunkten. 
Er. Nach allem bisher Gesagten ist nun der 
0. sich anirgendeinen Schaltpunkt anschließende BERLIN 
Bewegungsablauf klarunddamitderZeitpunkt un NOYSIENB Re 
B.:- gekommen, sich über die für einen gewünschten BER  Schaltpunkte a ER 
Verlauf bei gegebenen Anfangswerten zweck- "Bild 29. Große und kleine Grenzellipsen. Malen; 
B- mäßigen Parameter klar zu werden. In Ziff. 1 u 
>... „wurde gezeigt, daß die Bewegungsgleichung, wenn die EI SHEEhe Zeit t als unabhängige ‚Ver- 2 = 
—  änderliche benutzt wird, die Parameter D, w, b und .o enthält. Nehmen wir einmal an, daß D 
und w vorgegeben sind, d.h. daß wir ein schwingendes System vor uns haben, dessen Dämpfung 
und Eigenfrequenz gegebenist. DieDämpfung sei unzureichend und es werde verlangt, durch An- 
“ wendung einer intermittierenden Regelung das Abklingen zu verbessern. Man möchte also bei. 
E ' vorgegebenen Anfangswerten 9, und $, in möglichst kurzer Zeit auf vernachlässigbar kleine 
-J Werte und g kommen. Wir können nun b/w®und o wählen und uns außerdem für Regelsystem A 
' oder .B entscheiden. Wählen wir z. B: b/w? groß, so kann es vorkommen, daß die zu unseren Be 
größtmöglichen Anfangswerten gehörende Schaltwerte fast alle in das Gebiet der kleinen Grenz- NA 
ellipsefallen. Das würde bedeuten, daß auf Grund der Ergebnisse dieser Arbeit ( „ideale“ Schwarz- Be 
Weiß-Steuerung ohne Schaltverschiebungen) die Bewegung nach kurzer Zeit undefiniert wird Br 
oder sich gegen die periodischen Bewegungen aufschaukelt (vgl. Bild. 14 und Ziff. 6). Wenn wir Eu 
dagegen b/w® klein wählen, werden die in Frage kommenden Schaltwerte innerhalb eines Ge- EA 
bietes ihrer Ebene liegen, das so groß ist, daß die Grenzellipsen und die beiden Kurven der N 
periodischen Bewegungen ganz in seinem Innern liegen. Würden wir nun z. B. o negativ wählen, 
dann würden wir im System A für einen großen Anfangsschaltwertbereich Aufschaukeln, bei \ 
positivem o dagegen im allgemeinen Abklingen bis in das Gebiet der kleinen Grenzellipse hinein Ba 
haben. Außerdem wäre noch zu untersuchen,, ob die Abklinggeschwindigkeit groß genug ist. u 
Wählen wir System B, so würde für negatives o Abnehmen bis zur Ruheellipse, bei positivem 0 er 
für einen großen Bereich innerhalb der Kurve ‘der periodischen Lösungen Aufschaukeln ein- sun 
treten. Da die Ellipse der Ruhepunkte stets größer als die der Endpunkte ist (kleine Grenz- 
ellipse),; wird System B mit negativem o bestimmt ungünstiger als A mit positivem o. Ent- 
scheiden wir uns also für System A und positives o. Nun kommt die weitere Frage: Wie groß 
soll oe gewählt werden. Dazu müssen wir das Abklingen über der Zeit 2 betrachten, d.h. Bild 19 
bis 22 und Bild 26. Man ist nach diesen zunächst sehr geneigt, ein ziemlich großes o zu wählen. 
Dem steht aber folgende Tatsache entgegen: Wie wir wissen, kommt man recht bald in das 
- Gebiet der kleinen Grenzellipse, in dem die Bewegung unter der Wirkung einer Regelung ohne 
-  Sehältverschiebungen nicht mehr definiert ist. Schon in Ziff.3 wurde aber darauf, hingewiesen, 
daß man in diesem Gebiet wieder zu definierten Bewegungsvorgängen kommt, wenn man die 
Schaltverschiebungen berücksichtigt (vgl. dazu die im Anschluß än unsere Arbeit durchgeführte 
ArbeitvonK.Scholz’”)). Eszeigtsich nun, daß bei Berücksichtigung dieser Tatsachen kleinere 
Werte o vorzuziehen sind, weil sie die Bewegung innerhalb der kleinen Grenzellipse schneller 
zum Abklingen bringen ®). 


?) K. Scholz, Über die Bewegungen eines Systems von einem Freiheitsgrad unter dem Einfluß einer Schwarz- 
Weiß-Regelung mit Schaltverschiebungen, noch unveröffentlicht; Auszugim Bericht I, Flügge-Lotz, Z. angew. Math. 


Mech. Bd. 25/27, S. 97, Ziff. 3. 
8) Herr Fischel machte uns darauf aufmerksam, daß in einer noch zu veröffentlichenden Arbeit von Gol- 


Fin g gezeigt wird, daß die Bewegung für-sehr kleine Nachhinkzeiten in ein Kriechen gegen neh Di CeÜ u) übergeht, 
Eingegangen: 11. 6. 44. 
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 Stoßwellenverzweigung. 
DEREN AS Von Franz Wecken in Ha 
0.0.0 Be wird untersucht, welche Beziehung zwischen Einfallswinke 
Bar reflexion für die Entstehung der Machschen Welle entscheidend ist, und 
ist, Dabei zeigt sich, daß für eine starke einfallende Welle die von Oranz ı chardin, 
dagegen die von Schultz-Grunow angegebene Erklärung u en äg % PR " 
0 Ik is investigated,: which relation, between the angel of i idence and L ; intensity of U ‚shock 
 determining I rallın of the M ach wave in case of the reflection of shock waves, Jurthermore 
' Mach wave is determined in a single way. It is shown that the explanation of the p enomen 
by Cranz and Schardin proves correct in case of ineidence of strong waves, while the 
given by Schultz-Grunow comes Irue in case of a fin me. 
L’auteur cherche quelle relation entre l’angle d’incidence et l’intensite de l’onde de choc est decisive po 
la formation d’une onde de Mach. De plus il examine, si ceite onde est d#erminde univoquement. Pour k 
une onde forte le phenomene est explique par P’hypothösedeCranzetSchardin. Quant d une onde faible, 
on doit se servir de l’explication donnee par Schultz MINE ME RN er 
B sTo4 _pabore uecuenyerTcH, KAKAA BABUCHMOCTb Mey YIIOM TaneHuN MH HHTeHCHB- 
HOCTBI Napa IIPH OTpaskeHHuH TUNPaBIMyecKoro y1apa OÖyCIOBJIHBACT BOSHHKHOBEHHE BONHET 
Maxa u HACKONBKO TOCHENHAA ONHOSHAYHO OnpeneneHua. Ilpu 9TOM BbiAchAerca, YTO NaHHOe 
Kparnom un Mapıunom OoÖBAcHeHne 9TOTO ABAIeHHA IIPHTOAHO B cayyae CHAIBHOH BOAHEL, # 
. a oösacHenne Ill yısu-I'pynoB — B cayyae c1a6oH BoAHE. Eh n 


a3, 


w 
Fi ‘ = 


1. Entstehung der Machschen Welle, x N 


| E.Macht) beobachtete bei der Durchdringung von zwei kugelförmigen Knallwellen ds 
Auftreten einer neuen Welle von erhöhter Stärke, die seitdem als,,Machsche Welle“ bezeichnet 
‘ „wird 1%). Man spricht auch von der ‚V-Ausbreitung‘ (wegen der beiMachs Versuchen aufgetre- 
_ tenen V-förmigen Rußfiguren) oder schlechthin vom M ach-Eifekt. Dieselbe Erscheinung kommt 
zustande bei der Reflexion einer kugelförmigen Stoßwelle an einer ebenen Wand (Bild 1, m 
dem der Reflexionspunkt P, der anfangs an der Wand entlanglief, sich an einer Stelle Avon 
ihr ablöst und unter Bildung einer verzweigten Stoßfront in die freie Strömung hinauswandert. 
‚Dieser Effekt zeigt sich immer dann, wenn der Einfallswinkel « einen bestimmten kritischen 
Wert überschreitet, der von der Stoßintensität abhängt. Als Maß für 
letztere wählen wir das Druckverhältnis er 


_P 


x <22 
NE | ER BER | 
"Fir ..wo ?, ünd 9, die Drücke vor und hinter der Stoßwelle sind. ’ 
2% Da die Erscheinungen (für ein bestimmtes Gas) nur durch a und & bestimmt sind, ver- 
Br: anschaulichen wir sie in einem quadratischen Diagramm: 1> x>0, 0<sina<1 (Bild2; 


die Stärke der Stoßwelle nimmt nach rechts zu). Bei kleinem & kommt 
stets die reguläre Reflexion (rechts in Bild 1) zustande. Es muß eine 
Grenzkurve im Diagramm Bild 2 geben, bei deren Überschreitung von 
unten nach oben die reguläre Reflexion in den Mach-Effekt umschlägt. 
Eine Erklärung für die Entstehung des Mach-Effektes wurde 
bereits von Cranz und Schardin?) angegeben; eine davon ver- 
schiedene gab kürzlich F. Schultz-Grunow). Hier soll eine Ent- 
Ä scheidung zwischen beiden Auffassungen versucht und dabei die genaue 
Bine: Lage der erwähnten Grenzkurve ermittelt werden. 


zz 


” 
SIna—= 


2. Grenze der regulären Reflexion. 


Die reguläre Reflexion läßt sich bekanntlich durch Mitführung des Bezugssystems mit 
dem Reflexionspunkt als stationärer Vorgang beschreiben, und man kann die reflektierte Stoß- 
front leicht mittels der Busemannschen Stoßpolaren konstruieren %). Für ideale Gase ist sie 
nach W.Döring?) auch analytisch zu berechnen. Die folgenden für ideale Gase mit konstantem 
c = A ee ee e E 

= = y gewonnenen Ergebnisse dürften qualitativ weitgehend auch für reale Gase gültig 
v 

. bleiben. Die reflektierte Wellenfront ist, ebenso wie die anliegende Kopfwelle eines Profiles im 
Überschallstrom, durch die Anströmgeschwindigkeit tv, und die Umlenkung ® bestimmt (Bild 3). 
Diese Bestimmung ist bekanntlich zweideutig, solange 9 nicht den von w, abhängigen maximalen 
PARAT re \ 

‘) E. Mach und Mitarbeiter, Wien. Ber. 72 (1875), 78 (1876), 75 (1877), 77 (1878), 78 (1879), 98 (1889). 

'8) (Anm. b. d. Korr.): Daß eine infinitesimale stationäre Stoßwelle, gewöhnlich Machsche Linie genannt, bei W uest 
(Z. ang. Math. Mech. 28 (1948), 73—80) als Machsche Welle bezeichnet wird, kann irreführen. 

») C. Cranz und H. Schardin, ZS.f. Phys. 56 (1929), 163—169. 


°) F.Schultz-Grunow, Z. angew. Math. Mech. 28 (1948), 30—31. 
4) 1.c. Anm.3, Bild 4, B 
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Ke lich, ‚undes lag nahe, zu vermuten, daß X eben die gesuchte Grenze zwi- 
schen den Bereichen der regulären Reflexion und des Mach-Effektes dar-- Sea 3. 
‚ stellt. Diese Annahme ist implizite von Schultz-Grunow °) gemacht er 
‘ worden. Cranz und Schardin gingen dagegen von folgender Überlegung a aus. Am Dale Fe 
"punkt berühren sich zwei Gebiete 'mit den Drücken. ?, und p, (Bild 3). Ist nun die Stoßwellen- 
‚geschwindigkeit D eines Drucksprunges von ?; auf », größer als die Anströmgeschwindigkeit N 
so wird aın Reflexionspunkt der rechts herrschende höhere Druck p, sich nach links Bahn brechen 
und die Machsche Welle erzeugen. Die gesuchte Grenze wäre danach die Kurve 8, auf der 


- Kriterium also nicht anwendbar sein. 


\ gegeben‘ haben. Y.T7 


' der Wand. ?, und «’ ändern sich hierbei stetig; der Weg des Punktes P be- 


Sur # RL eibt: es eine; für > Dan Keine A Ä 
sung. Da w,, w, 9, max durch & und « bestimmt: sind, erscheint ‚die Gleichung le Omas 
R als Kurve K (nach Döring) in Bild 2. “ 
Lösungen, die sich Entecheib der Kulryk K stets geben, SER Be RR: 
: konntlich | in der regulären Reflexion (wie auch in der anliegenden a nur die. mit dem iR 

geringeren Druck (p,, Bild 3) realisiert („schwache Lösung‘‘). Eine physi- Re 


Von den zwei reellen 


kalisch befriedigende Erklärung hierfür dürfte erstmalig INRICHL EN | 


Oberhalb von X (Bild 2) ist sicher boss reguläre Reflexion mög- 


D=w, ist. Diese existiert in der Tat und ist in Bild 2 eingezeichnet. Hiernach ergibt sich 


folgende Vorstellung vom Eintritt des Mac h-Effektes: Unterhalb $ ist. die reguläre Reflexion 


stabil, da D< w, ist; eine Machsche Welle müßte, wenn sie sich gebildet hätte, zurückbleiben 
und wieder abgebaut werden. Auf $ ist bereits eine stabile Machsche Welle möglich, aber sie 
kann nicht wachsen, sondern der Gabelpunkt P schreitet parallel zur Wand fort (Bild 4). 
Wird bei zunehmendem & (Bild 1) S überschritten, so entfernt sich Pvon __ 


ginnt tangential zur Wand. Da 8 unterhalb von K verläuft, würde hiernach 
K nie erreicht werden und das von Schultz-Grunow angegebene 


Jedoch enthält diese Begründung eine Lücke. Das Gebiet unter- 


halb X (Bild 2) ist ja von der Mannigfaltigkeit der regulären Reflexionen -  Bilda 


zweifach überdeckt, und es ist zu fragen, in welchem der beiden 

„Blätter‘“ die Kurve $ verläuft. Nur das Blatt der schwachen Reflexionen ist ja physikalisch 
reell. Nun zeigt sich, daß X und ‚S einander in einem Punkt B (für Luft: x = 0,44; 
& = 41,3°) berühren, und daß $ rechts von B zwar schwache, links von B aber (in Bild 2 
gestrichelt) starke Reflexionen durchläuft. Links, d.h. für schwächere Primärstöße, wird 
also bei Vergrößerung von & bis K hinauf gar nicht $ überschritten, da 8 im anderen Blatt 
liegt; es bleibt hier überall D<< w, und damit die reguläre Reflexion stabil. Unser Ergebnis ist 
also dies, daß die wahre Grenze der regulären Reflexion die in Bild 2 links von B mit X, rechts 


‚mit 8 bezeichnete Kurve ist, oder: für schwache Stoßwellen (x > 0,44) ist das Kriterium von 


Schultz-Grunow, für starke (« < 0,44) das von Cranz-Schardin ‚maßgebend. 

Daß auf Klinks D< w,. rechts D > w, ist, hat auch Schultz-Grunow bemerkt (l. c. 
Bild 9), jedoch ist nach unseren Überlegungen nur Bild 9a (l. c.) reell, während der in 9b ge- 
meinte Zustand gar nicht erreicht wird. 

Während rechts von B auf S der Übergang zum Mach-Effekt stetig (hinsichtlich a’ und 
p,) erfolgt, ist das links auf K nicht möglich; &' und p, nehmen hier sprunghaft zu. 


3. Gegabelte Verdichtungsstöße. 

Es bleibt nun zu untersuchen, ob auf X, links von B, eine Stoßwellengabel mit den vor- 
gegebenen Werten x und « stets existiert, und ob sie eindeutig bestimmt ist. 

Dies führt auf die etwas allgemeinere Fragstellung: ein verzweigter Verdichtungsstoß (Bild 5) 
ist zu bestimmen aus x’und «a. Der Einfallswinkel « zwischen Primärwelle und Wand kann dabei 
ersetzt werden durch den um 90° größeren Winkel zwischen Primäryelle und 
Machscher Welle, da diese anfangs senkrecht zur Wand stehen muß. 

Die Bestimmung einer Stoßgabel durch zund «ist physikalisch nahezu 
zu realisieren durch eine auf eine Kante auftreffende ebene Stoßfront (Bild 6), 
wobei dann auch je nach dem Ort im x-&-Diagramm reguläre oder Mach- 
sche Reflexion eintreten wird. & wird freilich, wegen der Krümmung der 
entstehenden Machschen Welle, nicht genau gleich dem Winkel zwischen rum 
Stoßfront und Wand sein. 2 Bild 5. 


5) W. Döring und G. Burkhardt, Beiträge zur Theorie der Detonation. Luftfahrtforschungsbericht 


Nr. 1939 (1944), 8.120 ff. e 
6) H. Richter, Z. angew. Math. Mech. 28 (1948), 341—345. 


' gleicher Druck herrscht, so ist, wie leicht zu sehen, bei vorgegebenem thermodyn 
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Ein verzweigter oder gegabelter Verdichtungsstoß ist nun einem E 
in dem P ruht, stationär zu behandeln, jedoch nur bei Außeracht . die sich in 


diesem System in der Regel bewegt (Bild 5). Die Anströmung wi Ä 
Macht man die Annahme, daß im Abstrom hinter der reflektierten und der Mach 
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Ausgangszustand die Stoßgabel durch zwei Parameter bestimmt. Die genannte Annahme ist 
ER i ae Falle des Bildes 5 sicher erfüllt; denn beide Wellenfronten | 
verlaufen von P aus stromabwärts, und es‘ist also dahinter kein Ver- 
dünnungsfächer mehr möglich. Somit ist es sinnvoll, aus x und & die 
Gabel zu ermitteln. Aber die numerische Berechnung oder graphische 
Konstruktion erweist sich als äußerst umständlich und praktisch un- 
durchführbar. A 
Nun haben A. Weise”) sowie Eggink®), ausgehend von Beob- 
Bild 6. achtungen in Windkanälen, eine Theorie des stationären gegabelten 
Verdichtungsstoßes entwickelt, und dabei ein elegantes graphisches 
Konstruktionsverfahren angegeben, das jede Gabel aus x und w, zu ‚ermitteln gestattet. 
Es werden dabei jeweils zwei „‚Herzkurven‘‘ miteinander zum Schnitt gebracht; ‚dabei 
ist eine Herzkurve (nach W eise) eine auf die Koordinaten 9 und 1/x oder log 1/x umgezeichnete 
Stoßpolare®?). 
Die Gesamtheit @ der oo? möglichen Gabeln, über die man sich nach der Herzkurven- 
methode einen Überblick verschafft hat, läßt sich dann in das Diagramm Bild 2 einordnen. 
Es wurde so für ein ideales Gas mit y = 1,405 (Luft) folgendes Ergebnis gefunden. Die 
zweidimensionale Mannigfaltigkeit @ der gasdynamisch überhaupt möglichen Stoßgabeln be- 
steht aus) zwei getrennten Teilen G, und G,, zwischen denen es keinen stetigen Übergang gibt. 
Die einfachste Darstellung erhält man, wenn man neben x als zweite Koordinate das Druck- | 


verhältnis y = Pe der reflektierten Stoßfront °) benutzt: jedes der beiden Blätter @, und @, 


überdeckt lückenlos und einfach das Quadrat 0 <z<1, 0 <y<< 1. Man stößt also bei stetiger 
Abänderung der Gabelgrößen (Winkel, Drücke, Geschwindigkeiten) dann und nur dann auf 
eine Grenze, wenn eine der Stoßwellen unendlich stark oder unendlich schwach wird, 

Bei y—1entartet die reflektierte Welle zur Machschen Linie; die primäre bildet mit der 
Machschen eine Gerade, es gilt «— 90°. So müssen in Bild 2 derlinke, rechte und obere Rand 
Bilder der entsprechenden Grenzen = 1, <= 0, y=1 von @, und @, sein. In der Tat über- 
deckt G, sowie auch @, einen Teil der x-&-Ebene lückenlos und einfach (Bild 7); die vierte Grenze 
y = 0 erscheint als je eine Kurve für @, und @,. Beide Kurven liegen 
spiegelbildlich zu «= 0. Bezeichnen wir mit @, das Blatt, das die 
senkrechten Gabeln (Kurve 8) enthält, so reicht @, weit in den Be- 
reich negativer & hinein, während der untere Rand von @, (im Diagramm) 
S dicht links von B schneidet, 

Da also K links von B größtenteils von @, und @, überdeckt. 
wird, ist die hier aus der Reflexion entstehende Gabel durch x und & 
nicht eindeutig bestimmt. 

Da nun rechts von B der Übergang stetig auf 8, also in das 
Blatt @, erfolgt, wird man aus Stetigkeitsgründen vermuten, daß auch 
links nur die Gabeln aus @, realisiert werden, womit die Eindeutig- 
keit hergestellt wäre. Doch ist diese Begründung physikalisch un- 
zureichend, da die Stetigkeit auch sonst verletzt wird. Dagegen kann 
man aus dem Umstand, daß bei festem x und «@, stets kleineres % 
(also größeres 9,) als @, ergibt, schließen, daß @, nicht realisiert wird. 

Bild 7. Denn der Druck p,, der bei Überschreitung von X einen Sprung macht, 

wird nur soweit steigen, als es zur Gabelbildung erforderlich ist; ein 

höherer Druck würde sich nicht halten, da die Abströmung Unterschall ist und nachfolgende 
Verdünnungen ihn abbauen müßten 19), 


. 
Nr, 
fr 


") A. Weise, Theorie des gegabelten Verdichtungsstoßes. Luftfahrtforschung, Techn. Berichte 10 (1943), 
Heft 12, Anlage II A. 

®) H. Eggink, Über Verdichtungsstöße bei abgelöster Strömung. Luftfahrtforschungsbericht Nr. 1850 (1948). 

8) (Anm.b.d. Korr.): Vgl. hierzu: W. Wuest(l.c. Anm. la); F.Wecken, Grenzlagen gegabelter Verdichtungs- 
stöße (Z. ang. Math. Mech., erscheint demnächst). 
25 i ®) Diese ist bei Weise als „Nebenstoß b‘‘ bezeichnet, die primäre als Nebenstoß a, die Mach sche Welle als 

auptstoß. 

1°) Die gleiche Überlegung ist auch auf die schwache und stärke Reflexion anwendbar, leistet aber weniger als die 

von Richter (l.c. Anm. 6): angestellte Stabilitätsbetrachtung. 


h; & Er Wand. trifft, erfolgt En im ee Senioren Gebiet: von Bild 2, mit Mach- Elfekt 
y ER brarieten Gebiet. ‚In diesen Teilgebieten sowie auch rechts von B auf der Grenzkurve hängen 


u, a’ und 7, (Bilder 3, 4, 5) stetig von zund «a ab; links von Berfahren sie einen sprunghafte 
| Zuwachs bei Überschreitung der Grenze von unten nach oben. { 
- Es bleibt. hierbei offen, ob unter anderen Versuchsbedingungen auch ei dene 1% 


-tantekchen: können, die zwar als BE IEE seh möglich, ar im Mach-Eiffekt als nicht reali- 
siert ‚erkannt yurden. E 


Ken Eingegangen: 14. 5. 48. 
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Die Stabilität des ee in einer konkaven Ecke, h 


Von Hans Richter in Haitingen. 


Bei der Berechnung des Verdichtungssioßes in einer konkaven Ecke oder bei Annliihen datlondireh, v 
7. Problemen ist die Lösung im allgemeinen zweideutig (,„‚schwache‘‘ und ‚‚starke‘“ Lösung). Die bisher üblichen 
Begründungen für die Bevorzugung der schwachen Lösung werden abgelehnt. ‚Statt dessen wird eine instatiomäre, 
rein mechanische Stabilitätsbetrachtung durchgeführt, die sich auf die Abhängigkeit der Stromlinienknickung 
von der Stoßintensität stützt. Diese Se zeigt, ne ein an einer starren Wand amsetzender ER 
starker Stoß instabil ist. EN 


Computing the shock of. compression in a concave corner or similar on ne the solution NEE 
is in general two-valued (‘“faint’” and “strong” solution). T’he arguments till now usual for the preferene 
given to the faint solution are rejected. In place of that, an instationary, purely mechanical inquiry of sta- 
biltiy is made based on the dependence of the breaking of the sircam-lines on the a iR “he shock. Th 
inguiry shows that a strong shock coming from a rigid wall is unstable. 


En caleulant le choc de compression dans un angle concave (ou en traitant des an semblables) A, 
on trouve en general deux valeurs (solution « faible» et solution « forte »). Les arguments pour la preference 
de la solution faible, jusqu’a maintenat en usage, sont refusees. L’auteur donne ü sa placerune recherche de 
stabiliteE purement mechanique et non-stationnaire, qui se base sur la dependance du brisement des lignes de 


couramt de Uintensite de,choc. La recherche montre que le choc provenamt d’un mur rigide est instable. Bl 
i x u 
Ilpn packere runpasıugecnoro ynapa o BOTHYTEIH yronN HIM MOXOKUX CTANNOHAPHEIX RR: 
sayay MONYYAIMTCH B OÖINeM CIYy4yae Ba Pa3HEIX pemenHus (,‚caa60e‘‘ u „‚cunsHoe‘“). B aTofi er: 
pa6oTe OTBepramTca IIPHHATHEe MO CHX NOP OÖ60CHOBAHHA IA IpenuoyreHnnn „cuaboro‘‘ RE 
pemeHus. BMmecTo 9TOTO 9TH ABIeHHA PACCMATPUBAIWTCA KAK HECTAIMOHAPHHE YMHCTOMeXa- . 


HuyeckHe CIyYau YCTOÄYHBOCTH, IIPHYeM 9TO UCCTENOBAHHe HCXOAUT OT 3ABHCHMOCTH U3JIOMA 
AIHHMUM TOKA OT UHTEHCHBHOCTH yAapa. ITO HecHeNoBaHne YCTOHYHBOCTH NOKASHBAET, ITO br 
yrap 0 HecMemamımmyioca TBePAyM CTeHy HeyCTqüyuB. 


1. Die Bevorzugung des schwachen Stoßes. 


Liegt ein keilförmiges Profil mit nicht zu großem halben Öffnungswinkel ® in einer zwei- 
dimensionalen parallelen Überschallströmung, deren Richtung mit der Profilachse übereinstimmt, 
so bildet sich bei genügend großer Mach scher Zahl M, der Anströmung eine anliegende Kopf- 
welle aus. Die Richtung und die Stärke dieser Stoßfront an der Profilspitze ist dadurch be- 

stimmt, daß die Stromlinien um den Winkel # abgeknickt werden müssen. Berechnet man — 
ee . etwa unter Benutzung der idealen Gasgleichung — den dieser Forderung entsprechenden Ver- 
dichtungsstoß, so erhält man mathematisch zwei mögliche Lösungen. Die Lösung mit dem 
größeren Verhältnis der Drucke vor und hinter dem Stoß wollen wir die ‚„‚starke‘, die andere 
die „‚schwache‘‘ Lösung nennen. Es ist bekannt, daß im ERpPLImER! stets die schwache Lösung 
realisiert wird. 


Als Begründung für dieses Verhalten wird oft die Tatsache herangezogen, daß unter 
konstant gehaltenen Anströmdaten beim Grenzübergang zu = 0 die schwache Lösung die- 
jenige ist, bei der das Druckverhältnis gegen 1 geht, während die starke Lösung für dieses Grenz- 
fall ein endliches Druckverhältnis < 1 liefert. Dieses Argument ist natürlich nicht befriedigend, 
da es erstens ja senkrecht stehende Verdichtungsstöße, also solche mit # = 0, unter geeigneten 
Versuchsbedingungen tatsächlich gibt und zweitens bei größer werdendem # ja trotzdem ein 
unstetiges Umklappen der einen Lösung in die andere aus sonst unbekannten Gründen denk- 
# bar wäre, so wie man dies in der Theorie des Mach-Effektes kennt, wo die gewöhnliche Reflexion 
=... einer Stoßfront an einer starren Wand bei Veränderung des Einfallswinkels plötzlich in eine 
= -Mach- Reflexion umschlagen kann, obgleich bei den betreffenden Einfallswinkeln auch die ge-" 
wöhnliche Reflexion mathematisch noch möglich ist-!). 


1) vgl. F.Wecken, Z.ang, Math. Mech. 28 (1948), 3383—341. 


‚Sehen w 
N: 


n wir von der vollkommen unzure 
on der schwachen und der starken Lösung einfa 
die mitunter angezogene Möglichkeit, einen Grund für die Bevor 
ET darin zu suchen, daß eine bestimmte Größe wie etwa die Entrop 
aber erstens gar nicht klar, wie überhaupt im vorliegenden Fal \ 
Thermodynamik. streng zu einer Entscheidung zwischen den beiden Lösungen hı 
‚werden soll, da das System nicht abgeschlossen ist. Zweitens muß beachtet werd 


gerade beim starken Stoß ein größerer Entropiesprung stattfindet als beim schx 
Stoß, daß also gerade die Lösung mit der geringeren Entropievermehrung 
‚siert wird. Drittens aber kann man leicht durch die folgende Über egung einsehen, 
"man mit derartigen stationären Betrachtungsweisen prinzipiell zu keiner Lösung des Prob E 
_ kommen kann: GE NE TEE SR. a 
In der gegebenen Parallelströmung möge ein so stumpfes Profil liegen, daß eine abgelöste BA 
Kopfwelle entsteht. Die längs der nach vorn verlängerten Profilachse ankommende Stromlinie _ Br 
_ erfährt an der Kopfwelle keine Abknickung; mit wachsender Entfernung der Stromlinie von 
der Profilachse wird ihre Abkniekung zunächst monoton wachsen, um dann nach Erreichen 
eines Maximalbetrages wieder monoton bis zu Null abzunehmen. Zu jeder überhaupt n dr 
. ‚abgelösten Kopfwelle vorhandenen Abknickung ® sind somit die schwache und die starke Lö- 
‘sung beide an der Kopfwelle je einmal realisiert, wobei die Realisierung der starken Lösung 
näher an der Symmetrielinie des Bildes liegt. Wir ersetzen jetzt eine Stromlinie, die etwas ober- 
- halb der Symmetrielinie liegt, durch ein starres Profil, für das die betrachtete Kopfwelle nun- 
mehr als anliegende Kopfwelle erscheint. Für die Strömung oberhalb des Profils wird 
durch diese Ersetzung nichts geändert. Wir müßten daher eine Strömung erhalten, beideran 
der Profilspitze die Kopfwelle mit der starken Lösung ansetzt. Für die durchgeführte ein 
stationäre Betrachtungsweise würde dies Strömung dann ebenso stabil sein, wie dies für die urs 
 abgelöste Kopfwelle gilt, von der wir ausgegangen waren. Die Entscheidung zwischen starker Fa 
und schwacher Lösung kann daher stationär nicht getroffen werden. Im Gegenteil erscheint 
‘ die Bevorzugung der schwachen oder der starken Lösung an der Profilspitze.hier als eine Frage 
. . der Profilgestaltung stromabwärts; m. a. W. man könnte durch geeignete Profilgestaltung die 
starke Lösung erzwingen. i - 

Dies ist nun, wie wir im folgenden zeigen wollen, deshalb falsch, weil eine der- 
artige Ersetzung einer Stromlinie durch eine starre Begrenzung zwar sehr wohl für eine 
stationäre Betrachtungsweise ohne Einfluß ist, während für eine instationäre Stabilitäts- Kile 
betrachtung ein wesentlicher Unterschied zwischen Stromlinie und starrer Begrenzung besteht. 
Dabei braucht in die Betrachtung nicht die Tatsache der Reibung des Gases an der starren = 
Wand einbezogen zu werden. f | 


2. Die Herzkurve. 
Re Zur Auffindung der schiefen Verdichtungsstöße, die bei gegebenen Werten von Druck ?,, 
PR spezifischem Volumen Y, und Geschwindigkeit W, der Anströmung die Umlenkung der Strom- 
linien um einen Winkel ® bewirken, ist es zweckmäßig, die Beziehung zwischen Druck 2, hinter 
2 dem Stoß und ® bei konstant gehaltenen Anströmdaten aufzustellen. Die entsprechende Kurve 
Re in einem d—p,-Diagramm wird oft als ‚„‚Herzkurve 2)‘ bezeichnet, ; 
Für den geraden Verdichtungsstoß erhält man zunächst aus den Erhaltungssätzen 
von Maße und Impuls die bekannten Formeln: - 


REIT ET : 
‚= YpVı- RE, den SE EN ER TR U (2.1). 
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uw=w—w=7pmh Ya—1) i—u) NE 


Hierbei sind w, und w, die Geschwindigkeiten vor und hinter dem Stoß, « ist der Geschwin- 


HR N Wasl : 
digkeitssprung an der- Stoßfront, u 7 ist die reziproke Volumkompression. a als Funktion 


1 
von betrachtet ist die Hugoniot- Gleichung. 
Für einen schiefen Verdichtungsstoß mit den Geschwindigkeiten W, und W, vor 


‚und hinter dem Stoß bleibt die Geschwindigkeitskomponente parallel zur Stoßfront beim 
Durchgang durch dieselbe ungeändert, während die zur Stoßfront senkrechte Komponente 


®) Vgl. A. Weise, Theorie d. gegabelten Verdicht n \ 
H.12, Anl. ITA. ’ ungsstoßes. Luftfahrtforschung, Techn. Ber. 10 (1943), 
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Ss als s Funktion v von aufzufassen. Wegen ie = ai ist dabei nach 0. RR 
Ko Bun FE N n-1 
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erden. nun m v, ons m, festgehalten, so hängt 0) an (2. 3) von # EN! ee die: Re 
Herzkurve definiert ist. Wegen (2.4) ist (9)a,1 = 0. Bei wachsendem x wächst zunächst auch 
9. Da für Verdichtungsstöße u <1'ist, kann der Nenner von (2.3) im Realitätsbereich des 
RE Zählers nicht verschwinden. Es ist also tg ® endlich. Schließlich wird BEE 0 0, ‚wenn für 
D nn Be RT j Ba £ 
Bee, u re 
E:: RE Ya a1, m 
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= | werden ist. Tier zugehörige Wert max entspricht dann a Raketen Stoßfront. Wie | 
er: - ‚man leicht aus (2.3) nachrechnet, beginnt die Herzkurve bei m = 1 mit der endlichen Nr 


) 2 94: 2 en En 
er m Rt t; während bei mas gilt: 1: a 


PZRmax. 
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2 In Bild 1 ist eine Herzkurve schematisch dargestellt, wobei angenommen wurde, daß BR 

auf der Herzkurve als Funktion von » betrachtet nur ein Extremum besitzt. Wie man sich aus Ne 
4.(2- 3) leicht ableitet, ist dies z. B. sicher dann, aber' nicht a 
‘nur dann, der Fall, wenn die Hugoniot-Kurve zur - 
Achse konvex verläuft, so wie dies bei idealen Gasen der 
- Fall ist. 


Die Herzkurve charakterisiert die Gesamtheit aller 
schiefen Verdichtungsstöße, die bei festgehaltenen Anström- 
daten möglich sind. Stöße, die die entgegengesetzte Neigung 
zur Anströmung haben, entsprechen dabei natürlich dem 
Spiegelbild der gezeichneten Herzkurve an der z-Achse.. Der ‘ 
Winkel zwischen Stoßfront und Anströmrichtung ist bei 
z=1 durch den Machschen Winkel gegeben. Gehen wir 

5. ‘ nun längs der Herzkurve zu größeren z-Werten über, so wird 
Er die Stoßfront immer steiler, bis sie bei % = max senkrecht 
zur Anströmüung geworden ist; gleichzeitig nimmt dabei‘ 
W, ab. 

Normalerweise besitzt die Herzkurve die folgenden charakteristischen Eigenschaften, ie 
man insbesondere bei Verwendung der Abelschen Zustandsgleichung mit konstanter spezi- 
fischer Wärme nachweisen, kann: 


1. Durchläuft man die Herzkurve von z=1 bis 7 = max; SO 


a) wächst 9 monoton bis zu einem d,, mit zugehörigem 7%, um dann monoton wieder zu 
fallen, 
b) fällt die Machsche Zahl 7 der Abströmung monoton von M;= M,>1 bis zu einem 
„Werte: M, <1. 


Y 08 VrE j j | u { =; 

St 2. Der Punkt m,, wo M,;=list, ht m <am 0 a 

Beet. . Diese Eigenschaften der Herzkurve, resp. des betreffenden Gases, se 1 
vorausgesetzt. | | 


der Kopfwelle einem Punkte auf der Herzkurve entsprechen. 

» Beieiner abgelösten Kopfwelle gehört dann der Punkt auf der nach vorn verlängerten Profil- 
achse zum Punkte = %nax auf der Herzkurve. Durchlaufen wir nun die Kopfwelle nach außen 
bis ins Unendliche, so wandert der entsprechende Punkt auf der Herzkurve bis zu deren Spitze 
r—=1. Bei einer anliegenden Kopfwelle ist dagegen nur der Teil der Herzkurve unterhalb 
eines solchen Punktes realisiert, bei dem die Abknickung ® gerade gleich dem halben Profil- 
öffnungswinkel ist. Wie Bild 1 zeigt, gibt es bei einer normalen Herzkurve genau zwei 
solche Punkte, wobei der ‚schwache‘ 7, < 7, und der starke 7, > m hat. 

Die beiden Lösungen x, und , unterscheiden sich nun wesentlich’durch den Sinn der 


Änderung von ® bei kleinen Änderungen von #. Dies hat zunächst die folgende bemerkens- 


werte Konsequenz: 

Es liege in der Parallelströmung ein dreieckiges Profil (Bild 2) mit dem halben Öffnungs- 
winkel 9, das gegen die Anströmrichtung um den kleinen Winkel e gedreht ist und an der Spitze © 
drehbar festgehalten wird. Entspricht nun die Kopfwelle der schwachen Lösung, so herrscht 
auf der Oberseite des Dreiecks. entsprechend der Ablenkung # + ein höherer Druck als auf 
der Unterseite, wo die Ablenkung ® — e beträgt. Das Profil wird daher in die symmetrische 
Lage zurückgedreht. Entspricht dagegen. die Kopfwelle der starken Lösung, so gehört zu 
der Ablenkung ® -+e auf der Oberseite gemäß Bild 1 ein kleinerer Druck als auf der Unterseite, 
so daß das Dreieck die Tendenz hat, sich weiter von der Symmetrielage zu entfernen. Allerdings 
wird dieser Effekt durch die Strömungsverhältnisse an der Basis des Dreiecks beeinflußt werden. 
Immerhin weist dieses eigenartige Verhalten der starken Lösung wieder darauf hin, daß sie in 
Wirklichkeit nicht auftritt. 


Bild 2, Blid 3. 


Um nun hierüber eine Aussage zu erhalten, betrachten wir zunächst einen unendlich 
langen stationär stehenden geraden Verdichtungsstoß. Im Abstromgebiet liege ein punktförmiges 
Hindernis, das also die Strömung nicht beeinflußt. Jetzt stellen wir uns vor, daß dieses Hinder- 
nis sehr langsam aufgeweitet wird. Die mit der Stauung an dem Hindernis verbundene Druck- 
erhöhung wirkt auf die Stoßfront, da hinter derselben ja Unterschallströmung herrscht. Der 
Erfolg ist der, daß die Stoßfront an der dem Hindernis zunächst liegenden Stelle etwas nach vorn 
gedrückt wird, während sie im'weiteren Verlauf sich wieder asymptotisch der alten Lage nähert. 
Ein solches Hindernis hat also die Wirkung, die gesamte Stoßfront nach vorn zu drücken, wo- 
bei seine Wirksamkeit mit wachsender Entfernung abnimmt. 


Denken wir uns jetzt eine konkave Ecke (vgl. Bild 3) mit dem Ablenkungswinkel 9, 
in der der Stoß x steht. Wir werden die Strömung stabil nennen, wenn sie bei Anbringung 
eines kleinen Hindernisses im Abstromgebiet in eine von ihr wenig verschiedene stationäre 
Strömung übergeht. | 

Nehmen wir zunächst an, daß wir es mit der schwachen Lösung zu tun haben: z = I. 
Es kann dann sein, daß im Abstrom Überschallströmung herrscht; dann kann durch ein kleines 
Hindernis hinter der Stoßfront dieselbe nicht beeinflußt werden. Herrscht dagegen Unterschall- 
abströmung, so würde die Anbringung eines allmählich wachsenden Hindernisses die Stoß- 
front nach vorn drücken. Wenn die Lösung stabil sein soll und das Hindernis als sehr klein an- 
genommen wird, so daß es die Stoßfront nicht in die senkrechte Lage bringen kann, dann muß 
der Ansatzpunkt der Stoßfront bei O bleiben; denn ein Vorrücken der Front entgegen der An- 
strömung als senkrecht zur Wand stehende Front würde bereits eine bestimmte endliche 
Drucksteigerung erforderlich machen, die durch ein sehr kleines Hindernis bei einer stabilen 


3. Die Stabilitätsbetrachtung. $ N 
Da für alle Punkte einer Kopfwelle die Anströmdaten die gleichen sind, muß jeder Punkt % 
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' in die schwache Lösung übergeht, für den dann Stabilität herrscht. Es braucht dies aber nicht 
unbedingt so zu sein, da vorher durch Stauungen im PRsenmeehiet der Prozeß de en ee 


rt eine Verdünnung eintreten würde. De Stärke der neuen a R 
e sein, da ja nach wie vor bei O um den Winkel ®: abgelenkt werden muß. 
er ebd dafür, daß. die zunächst erfolgende Verstärkung der Stoßfront sofort wieder rück- 
 gängig gemacht wird, ist nun darin zu sehen, daß gemäß Bild 1 zu einer Stoßfront größerer 
Intensität auch eine größere. Ablenkung der Stromlinien gehört, so daß bei der Bildung der 
‚stärkeren Front in ‘O sofort ein freier Winkelraum zwischen Wand und Abströmung auftri 
. dessen "Ausfüllung - zu einer Verdünnungswelle mit Abbau der. "Stoßintensität rührt. Durch 
E, Mechanismus gt ‚das stabils Yorhalin der Stoßfront gegen kleine un gewähr- 
tet. ER, ‚* 
Nehmen wir jetzt unlesn a daß wir die Each Tlösungs vor uns haben: MIT . Dann 
aerscht nach unserer Voraussetzung über die Herzkurve jedenfalls Unterschallabströmung, 
so daß ein Hindernis im Abstromgebiet wirksam wird. Aber jetzt funktioniert der oben erwähnte BR. 
. Ausgleichmechanismus nicht mehr, Eine Verstärkung der Stoßfront in O bedingt janach Bld1 . 
eine kleinere Ablenkung der Stromlinien, so daß jetzt hinter O eine zusätzliche Stauung auftritt, 
‚ die zu einer weiteren Verstärkung führt. Es kann daher keine Stabilisierung eintreten, sondern 
' die Anbringung eines noch so kleinen Hindernisses im Abstromgebiet führt zwangsläufig de | 
‘daß sich die Stoßfront beiO vollkommen aufsteilt und nach vorn läuft. Die entsprechende en. 
Dh, ist für den Fall durchführbar, daß ein negatives, Hindernis etwa dadurch geschaffen 
' wird, daß. man die Wand im Abstromgebiet etwas zurückzieht. In diesem Falle würde analog 
. eine Abschwächung der Stoßfront in O stattfinden, deren Betrag sich automatisch vergrößert. 
Dies kann theoretisch so weit gehen, daß die Intensität x unter rm sinkt und der Stoßeinfach 


‚macht werden kann. - 


Eine Bestätigung dieser Überlegung. wird durch den folgenden experimentellen Befund 
gegeben, dessen Kenntnis ich einer mündlichen Mitteilung von Herrn Oswatitsch verdanke. 
In der Strömung liege ein Profil, das aus einer dreieckigen Spitze mit anschließendem Rechteck 
besteht (vgl. Bild 4). Hierbei ist die der UN Na Seite des Rechteckes so 
groß gewählt, daß die Dreieckspitze 
die zum Rechteck gehörige abgelöste 
Kopfwelle nur wenig durchstechen 
würde. Auf diese Weise müßte man 
erzwingen können, daß sich an der 
Dreieckspitze eine anliegende Kopf- 
welle bildet, die der starken Lösung stationäre Starke | 
entspricht. Das Experiment zeigte, Aopfweile 
daß keine stationär stehende Kopf- 
welle zustande kommt, sondern daß / Aray 

. die Kopfwelle sich periodisch mit —- ---------- -2--- -- 
hoher Frequenz von der Spitze ent- 
fernt und wieder an dieselbe anlegt. 
Diese Tatsache entspricht nun voll- 
kommen der obigen Überlegung. 
Eine kleine Druckerhöhung bei Vor- 
handensein der starken Lösung sorgt zunächst für eine Aufsteilung und Ablösung der Kopf- 
welle. Da die gewählte Profilform aber stationär eine anliegende Kopfwelle fordert, muß die 
Kopfwelle wieder zurückkehren. Der mit der Rückkehrbewegung verbundene Unterdruck 
gegenüber dem Druck der starken Lösung sorgt jetzt dafür, daß die Intensität der wieder an- 
liegenden Kopfwelle weiter absinkt. Es kann sich aber nicht einfach der’ schwache Stoß an der 
"Spitze einstellen, da dies das angeschlossene Rechteck nicht zuläßt. Es treten also so hohe 
Stauungen auf, daß der Abbau der Stoßintensität zur Umkehr gebracht wird. Auf diese Weise 
kann die beobachtete periodische Bewegung der Stoßfront wenigstens qualitativ verstanden 
werden. 
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0 Es wurde das Temperaturfeld, hinter einer schwachen punkt- oder linienf "mige 
* Boden einer Reibungsschicht ausgemessen und die Versuchsergebmisse dieser turbulenten Wär 
in ein einfaches Formelsystem zusammengefaßt. N ae Ki a 


RER, He IRRE 
| The distribution of temperature in an air flow due to a faint point or line source of heat on the b 
of a frietional ER layer was measured and the results of this turbulent heat transfer are cond 
into simple formulas. Be ehe 


LD’auteur mesure la distribution de temperature derriere une source de chaleur (en forme de point ou ge}. 
de ligne) au fond d’une couche de fricion. Ensuite il nähe formules simples pour exprimer les ve 

resultats de cette recherche traitant une emission turbulente de c, 5 2 RN N IE 
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 TemNa Ha AHe IOTPAHHYHOTO CIOA ÖbLIM U3MepeHs U PesyabTaTk 9T0r0 TYPÖyAIeHTHOTO Por 50 
HPOCTPaHeHus Te ILIA HBIOMEHE B POPME IIPOCTEIX POPMy1. N ET De EAN FE $% ER 
1. Einleitung. ER 


. Im folgenden wird ein Bericht [1] teilweise wiedergegeben, der zwar bereits 1944 abge- Ina 

schlossen, aber damals nur in wenigen Exemplaren vom Auftraggeber (Marineobservatorium) 
verteilt wurde. Die nachträgliche, allgemein zugängliche Veröffentlichung erscheint dadurch 24 
‚gerechtfertigt, daß das Thema: stationäre turbulente Ausbreitung von Stoffbeimischungen 
‘ oder Wärme am Boden einer ebenen Reibungsschicht allgemein genug ist, und die Ergebnisse hr 
in verschiedenen Fällen anwendbar sein dürften. Zum Teil wird dieses Problem auch von 
B. Frost [2] in einer inzwischen erschienenen Arbeit behandelt, in der der Feuchtigkeitsgehalt 
des Windes, der über einen Ozean streicht, abgeschätzt wird. Frost untersucht jedoch nur den 
zweidimensionalen Fall und diesen nur theoretisch, wobei er hierfür zwar formelmäßig dieselben 
Ergebnisse erhält, jedoch mit anderen Zahlenwerten und Exponenten, als sich durch die Aus- 
wertung unserer Versuche ergeben hatte. % TREE 

Die turbulente Diffusion chemischer oder mechanischer Beimengungen und die turbulente | 
Ausbreitung von Wärme in Luft können gleichzeitig betrachtet werden, da man wohl den ent- 
sprechenden Austausch ey (für den Massenstrom) und &,; (für den Wärmestrom), gleichsetzen 
darf [3]: eu = &s. Der Impulsaustausch e, ist dagegen je nacb den kinematischen Einzelheiten 
der Strömung verschieden von ey und &5; so ist z.B. bei Wandturbulenz es/e, = 1,4 [4] und 
bei freier Turbulenz sogar &s/& = 2 [5] gemessen worden. Da es nun meßtechnisch einfacher ist, 
die Temperatur in einem Luftstrom zu messen als etwa die Konzentration einer Chemikalie, 
wurde folgende Modellströmung untersucht (vgl. Bild 1). In die Bodenplatte eines Windkanals 
mit rechteckigem Querschnitt wurde eine elektrisch geheizte Drahtspirale eingelegt, und mit 
einem Thermoelement die Temperaturverteilung hinter dieser Wärmequelle in der turbulenten 
Reibungsschicht längs der ebenen Kanalbodenplatte ausgemessen, und zwar sowohl im drei- 
dimensionalen Falleiner nahezu punktförmigen Wärmequelle wie im ebenen Falleiner Linienquelle. 

Zur Erhaltung der Modellähnlichkeit (ey = &5) muß dabei vorausgesetzt werden können,. 
daß die Unterschiede im spezifischen Gewicht einerseits der Beimengung und andererseits der 
erwärmten Luft gegenüber der kalten Luft der Hauptströmung so klein sind, daß die Eigen- 
geschwindigkeit v, der warmen Luftteilchen nach oben durch den Auftrieb gegen die horizontale 
 Windgeschwindigkeit U‘ vernachlässigt werden kann. Zur groben Abschätzung des Verhält- 
nisses v,/U im Versuchsbereich berechnen wir die Aufsteiggeschwindigkeit v4 eines kugelförmigen 
Turbulenzballens ‘mit dem Durchmesser , der gegen die Umgebung (20°C) um 3°C wärmer 
: Pen 0) Anl 
ist. Der Auftrieb ist dann A= I 5 793 (4 
Luft, g = 9,81 m/s?. In ruhiger Luft würde dieser Turbulenzballen mit einer solchen Geschwindig- 


3 N 
) ; e= Dichte der kalten umgebenden 


2 
keit v, hochsteigen, daß der Widerstand W = c 2 var (5) gleich A wird; der Widerstands: 


w 92 « 
koeffizient c,, (für die Kugel) hängt dabei selber noch von v4 1/v ab (v = kinematische Zähig- 
keit der Luft). Nimmt man z. B. 9 = 200° C und Z= 5 mm an, so wird v, = 0,12 m/s und C 
1,8 bei v‚ 1/v #40. Bei den Versuchen war U = 5,4 bis 30,9 m/s, so daß v,/U = 2% bis 
0,4% war. Diese Zahl gilt nur unmittelbar hinter der Wärmequelle; weiter hinten sind die 
Temperaturerhöhungen wesentlich kleiner. Der Wärmeauftrieb war also besonders bei den 
höheren Geschwindigkeiten vernachlässigbar klein. Dementsprechend ergaben sich auch die 
höchsten Temperaturen am Kanalboden; die Temperatur unmittelbar an der nahezu wärme- 
undurchlässigen Wand wurde zwar nicht direkt gemessen, doch war die Extrapolation des 
Temperaturverlaufs an die Wand vor allem bei größeren Entfernungen von der Wärmequelle 
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boden liegen ı müssen, ‚wenn are Se ae eine nero Aufsteiggeschwindigkeit erfahren. Mi 
R Bailey ‚wie z.B. nr ‚einem Kartoffelfeuer bei schwachem Wind. 
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2. Versuchseinrichtung. 


Das ne ist in Bild 1 skizziert. Die untere Kanalwand wurde dirk: glatte, Bee.. 


‚ebene, wärmeundurchlässige Sperrholzplatten gebildet; sie war 1,40 m breit und insgesamt6m 
lang. Die obere, gelenkig ausgeführte Wand wurde so eingestellt, daß der statische Druck im 
$ Kanal ulanp? der en: such reihe gleich dem Außendruck war. Die Kanalhöhe war so 
groß (40—60 em), daß zwischen den Reibungsschich- BE. 

‘ten an der unteren und der oberen Kanalwand noch 


‚ungestörte Strömung mit der Geschwindigkeit U war. 


: bungsschicht mit einer ähnlichen Geschwindigkeits- 

verteilung wie über einer Ebene (oder z. B. über einer 

Meeresoberfläche) bei gleichmäßigem Wind. 
Beieiner Rücklage von X = 2,93m (bzw. 2,70 m 
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Bild 1. Versuchsschema. 


denplatte war in einer kleinen, mit Asbest ausgelegten Vertiefung eine Heizspirale aus 0,8 mm 
starkem Chromnickelstahldraht eingelassen, die elektrisch geheizt werden konnte. Als Punkt- 
quelle wurde eine Spirale mit 3mm Außendurchmesser und 32 mm Länge verwandt; diese wurde 
der Länge nach in Strömungsrichtung angebracht, um dem Fall einer ‚„‚Punkt‘‘- Quelle vor allem 


hinsichtlich der Ausbreitung der erwärmten Luft nach der Seite möglichst nahe zu kommen. 


Als „Linien“-Quelle wurde eine Heizspirale mit demselben Durchmesser quer durch den gan- 
zen Kanalin die Bodenplätte eingelassen. Der Draht wurde bis zu schwacher Rotglut erhitzt. 


Die Temperaturverteilung hinter der Wärmequelle wurde mit einem Kupfer-Konstanten- 
Thermoelement bestimmt, und zwar wurde der Unterschied gegen die Temperatur an einer 
Stelle im Kanal vor der Heizspirale auf etwa 1/10°C genau gemessen. Die Heizung wurde so 
reguliert, daß .der Temperaturunterschied eines kurz hinter der Spirale aus der Bodenplatte 
ragenden Thermoelements und demjenigen vor der Wärmequelle konstant blieb. Die von der 
Heizspirale verbrauchte elektrische Leistung konnte nicht gleich der von der Quelle gelieferten 


Wärme gesetzt werden, weil vor allem bei kleinen Windgeschwindiskeiten ein Teil der Wärme- 


menge durch die Vertiefung im Kanalboden, in der die Heizspirale angebracht war, infolge 
Wärmeleitung nach außen hindurchging. Kurz hinter der Spirale fielen die Lufttemperaturen 
schnell ab, weshalb der Wärmedurchgang dort durch die 25 mm starke Sperrholzplatte ver- 
nachlässigbar war. Die an die Kanalströmung abgegebene Wärmemenge wurde aus dem ge- 
messenen Geschwindigkeits- und Temperaturfeld durch Integration ermittelt. 

Am. Ort der Wärmequelle war die Reibungsschicht je nach der Anströmungsgeschwindig- 
keit 50 bis 70 mm dick; sie stieg mit wachsender Entfernung von der Quelle in Strömungs- 
richtung langsam an entsprechend der bekannten Formel ö = 0,37 X/(UX/v)'®, X = Rück- 
lage von der Plattenvorderkante = x +2,93 m, wenn x den Abstand von ‚der Wärmequelle 
bezeichnet.. Die Temperaturen wurden bei der Linienquelle bis x = 75 cm, bei der Punktquelle 


‘ Auf diese Weise entstand an der Bodenplatte eine Rei- 3 y 


bei der Linienquelle) In, der Vorderkante der Bo- 


w 
et: 


nr 348 Wieghar dt, Über Ausbreitungs org 


a R a ee a 2 
"bis 62,5'cm gemessen. In diesem Meßbereich fand die Wärmeausbreitung noch 
der Reibungsschicht (etwa bis zu y/ö< 0,5) statt. Bi Bay 
Die Genauigkeit der Temperaturmessung war im Fall der Linienquelle bei 

strömungsgeschwindigkeiten sehr gut. Bei der Punktquelle, wo nur ge Wär 
I Verfügung standen, da sonst die Übertemperaturen in der Nähe der pirale zu g' 
er wären, ergaben sich naturgemäß größere Streuungen. Vor allem machte sich hier 


srößere 


u‘ 2 


ker eine Wärmeschichtung im Kanal von wenigen 1/10° C bei den kleinen zu messenden Tem N 
Bes differenzen störend bemerkbar. Diese Wärmeschichtung rührte davon her, daß die Luft im = 
Pro: Versuchsraum, die durch das Gebläse in den Kanal gesaugt wurde, verschieden warm wat; 


Br 


‘durch Abstellen der Heizkörper im Versuchsraum in der Nähe des Gebläseeinlaufs konnten 
die Temperaturschwankungen im Kanal verringert werden. . 


‚3. Versuchsergebnisse. 
a) Ansatz für die Auswertung. 


& 2, Zunächst führen wir folgende Bezeichnungen ein (vgl. Bild 1): 


x = Koordinate in Strömungsrichtung N ; > ER 
52 ande bie, Hohe De a A 
z = Querkoordinate - die’ Janienquslie ie A Aa 
vw = Geschwindigkeitskomponente in x-Richtung 
® TE ” »» Y- ” 
WORT) | ER) „» ar > =(), 
w = Geschwindigkeitsvektor. Der statische Druck p sei überall = 0. 9 (z, y, z) = Tempe- _ 
ratur an der Stelle &, y, z minus Temperatur der anströmenden Luft. » = kinematische 
Zähigkeit der Luft. 
X = Rücklage der Wärmequelle von der Vorderkante der Kanalbodenplatte. ö (x) = Dicke 
der Reibungsschicht. oe = Dichte der Luft, c, = spezifische Wärme der Luft bei kon- 
stantem Druck. . | 
red las?., ; 
o& = Impulsaustausch ’ 


085 Wärmeaustausch. | 25 
Dem Wärmeinhalt der Raumeinheit 0 9 c, 9 in einem Punkt hinter der Wärmequelle im 
Versuch entspräche bei einer Beimengung die Konzentration e in g*/M® (g* = g Masse). 
. Zunächst gilt für zweidimensionale Strömung die Kontinuitätsgleichung 


; ou 
; re ER NE a — 
div wm REN EEE ER EE LRNERS ES (1). 
u und v müssen die Bewegungsgleichungen erfüllen: 
w- grad w = div (e,grad tw), also genähert 
au u 0) =) 
— +9 — (* Fr NET EN EN we (2); 


en 0) au ; ou 
das Glied le e kann dabei wegen Er < weggelassen werden (Grenzschichtvernach- 


dy 
lässigung). Die entsprechende zweite Gleichung gibt wegen > v nur Aussagen über vernach- 
lässigbar kleine Größen. j 

Auf diese Gleichungen gehen wir aber nicht näher ein, sondern entnehmen « und » direkt 


den bekannten Versuchen über turbulente Reibungsschichten an ebenen Platten. Danach 
kann. man setzen: 


| 1 dö 
urU (y/ö)YR, VS FUSS EEE era 
(y/8) Ta 8) 
Für die Kanalströmung kann nach früheren Versuchen n = 7 gesetzt werden und Ö aus 
6 =.0,37 ZU UXyR EN ee (4) 


berechnet werden. Bild 2 zeigt ein solches Geschwindigkeitsprofil und die Annäherung durch das 
1/7-Potenzgesetz. Für andere Entfernungen von der Quelle & und Geschwindigkeiten U wurde ö 
nach der obigen Formel berechnet. Damit ist die Geschwindigkeitsverteilung im Versuch bekannt. 

. Für den Wärmeinhalt gilt zunächst ähnlich der Kontinuitätsgleichung für die Geschwindig- 
keiten ein Wärmeerhaltungsgesetz, wonach durch jeden Querschnitt & = const. hinter der 
Quelle (wegen » < u) dieselbe Wärmemenge pro Zeiteinheit hindurchfließt: 


oo oo © + © 
in. [ f 09,uddydz=2 f 0 yuddydz=const. .... (5). 
v=0 2= yv=0 2=0 eh 


ER e von irale er Zefleinheita an die SEORE abgegebene Wärme 

ee Br 'erner gilt für die Wärmeausbreitung entsprechend den obigen Bewegungs- 
- gleichungen die folgende: Differentialgleichung, da man gegenüber dem konvektiven Wärme» 
' transport die Ne und die reine  Wärmediffusion Ben kann: Re een Kar: = 3% 
div (eo ar 9 also ; N 


N ee 29 5 0% "(98 7 “ EN 
ER N wagte [77 % = ,,(& a 26 un BR 


ei 


har) o 
wobei sioder — (e Ei vernachlässigt wurde. | 


RR 


de\. dx g. 
&% Im zweidimensionalen Fall der Linienquelle ist also die Gleichung für 9 (6) der Bewegungs 
> gleichung für ı (2) sehr ähnlich. Jedoch sind die Randbedingungen ändere: 
Ber 3 
SER ‚y=0:u=v—=0 = =0 (da die Mann ‚ wärmeundurchlässig 2 
ut .$ 7 e Rn Re R 
ee. - DENE: ER, 0. Nor RAR IE 


Im manchen Fällen kann man nun näherungsweise & proportional & Setzen, one &, aus a 
experimentellen «, v-Verlauf zu berechnen ist, und kann # rechnerisch ermitteln. Hier würde 
N das aber z. B. im Fall der EUDEESHEUE abgesehen von den verschiedenen Randbedingungen, für 
u und Öschon auf Schwierigkeiten führen, da “ 
&, nicht von z abhängt, es am seitlichen Rand 
der Wärmeausbreitung jedoch verschwindet. 
Da also ein ausreichender Ansatz für es nicht 
zur. Verfügung steht, beschränken wir uns im 
folgenden auf die reine Versuchsauswertung. Br 
Bei anderen Versuchen über turbulente A 
We bei denen aus einem Spalt ER 


EEE ENDEN de as am . 
UN \ . + N NR A in,‘ 
i \ “ \ 1 ; z 
5 - Ur y 


UN WEN PER Ze Er 


Rücklage X = 2503m 


Geschwindigkeit U= 180 m/s 
UXx/v = 290:10® 


ENTE ERLITT LA 


© ohne Turbulenzgitter 
e mit nd 
4 u/U= (y/s)'l7 


002 
10 20% 4:30 40 50 02 04 066 08 ı ‚2 nizilem] # 
Ylmm] Bild 3. Temperaturprofile hinter der Punktquelle 
Bild 2. Geschwindigkeitsprofile in der Reibungsschicht. bei U = 17,8m/s und x = 25cm. 


3 warme Luft in die Reibungsschicht einer umströmten ebenen Platte ausgeblasen wurde [6], 
hatte sich folgendes ergeben: Die Temperaturprofile 9 (x = const, y) waren in verschiedenen 
Abständen x zueinander affin und ließen sich durch 9 © exp (—y*) annähern. Es ist daher 
naheliegend, auch das hier vorliegende Versuchsmaterial auszuwerten nach dem Ansatz: 


»—= 9, (z)exp {—-(yYlfı (e))*—(lz|/fa(z))P} für die Punktquelle...... (7): 


im Fall der Linienquelle wird f,—. 
23 
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b) Versuchsergebnisse. a it, A 
Mit diesem Ansatz wurden die Messungen folgendermaßen ausgewertet. In den Abständen 


v 7,5, 12,5, 25, 37,5, 50 und 62,5cm wurden die Mittelprofile 9, = 9 (y) für 2 = const 


i i a n bei eiı 
— 0 gemessen und durch Extrapolation an die Wand 9, gefunden. Ferner wurde einigen 
werten auch Querprofile in verschiedenen Wandabständen y gemessen: ® (2) für x = const, 


y= const. (Der Warmluftstrahl pendelte manchmal seitlich etwas, so daß das Temperatur- 


maximum nicht genau bei z = 0 lag; in diesen Fällen wurde bei der Auswertung der Paar au 
der z-Koordinate entsprechend um wenige mm verschoben.) Dann wurde —log Om/do = R % 
sowie für die Querprofile —log ®/Omax (2 = 0) über 2] doppeltlogarithmisch Fa ae ar 
gibt als Beispiel eine solche Auftragung wieder. Man sieht zunächst, daß die Querpro 
sächlich ähnlich sind, da für drei verschiedene Wandabstände die Kurven Ddm rear 
fallen. Die Richtigkeit des Ansatzes (7) wird weiter dadurch bestätigt, daß sich in dieser er 
tragung Geraden ergeben, wenn man von den äußersten Randpunkten absieht, wo aber eo 
die prozentuale Meßgenauigkeit nicht mehr so groß ist wie im Inneren des Warmluftge ” . 
Die Neigung der Geraden gibt die Exponenten « und ß; aus den Ordinaten a und bbeiy=1 cm 
bzw. 2= 1 cm können f} = (0,434/a)V« und f, = (0,434 /b)P berechnet werden (log e = 0,434). 
& und ß zeigten bei der Punktquelle keine systematische Abhängigkeit von der Entfernung ee 
der Quelle « und der Windgeschwindigkeit U; es ergab sich stets « = ß = 1,8. Nur bei der 
höchsten Geschwindigkeit U = 30,9 m/s erhielt man bessere Annäherung durch «= ß=2,0. Bei 
der Linienquelle war infolge der höheren Meßgen auigkeit eine schwache Abhängigkeit desavon U 
zu erkennen ; «lag zwischen 1,6 bis 2,0, siehe Bild 4. Hier genügte es, Mittenprofile zu messen, nach- 
dem durch Stichproben festgestellt war, dab , 

nicht nur die Geschwindigkeits-, sondern auch 
die Temperaturverteilungunabhängigvonzwar. + 


o ohne 
oe mit 


Turbulenzgitter 


” 6 8 


15 20 30 103 2 2 Urlv 
r [mis] Bild5. Höhe f d Breite j, des T turfeldes hinter 
N 2 = öhe f, un rei ‚ des Temperatur e 
Abb.4, Linienquelle: & abhängig von U. Br Punktquelle :fı = 3,90 Y/(U%/r) 1% und h—0,15em = 
4,75 z/(Ux/»)?’® = 1,64 5. 


& und ß bestimmen nur die Art der AbklingungungSfunktionen. Ein Maß für die turbu- 
lente Ausbreitung bilden erst die Längen f, (x) und fs (x). Sie geben an, in welcher Höhe bzw. 
in welchem seitlichen Abstand in einer gewissen Entfernung x hinter der Quelle die Temperatur 
auf den e-ten Teil des örtlichen Temperaturmaximums (9/9, = 1/9,72 = 0,368) abgefallen ist. 
(Bei & bzw. ß = 1,8 wird 9/9, = 0,1 für y= 1,59 f, bzw. z= 1,59 f,.) f} und f. hängen nicht 
nur von x,sondern auch von Uab. Um diese Abhängigkeiten zu erfassen wurden die dimensions- 
losen Reynolds-Zahlen Uf,/v und Uf,/v über Ux/v doppeltlogarithmisch aufgetragen. Es zeigte 
sich, daß dabei die Punkte U f,/v beikleinen x-Werten systematisch von den übrigen abwichen. 
Das erklärt sich daraus, daß dort f, noch nicht groß gegen die seitliche Ausdehnung der Heiz- 
spirale (3 mm) ist. Um daher auf den Fall einer wirklichen „Punkt‘-Quelle näherungsweise zu 
extrapolieren, wurde als Korrekturgröße vonder gesamten seitlichen Ausbreitung 2 f, die Dicke 
der Spirale 0,3 cm abgezogen; dann fielen für alle Entfernungen und Geschwindigkeiten die 
Punkte U (f,—0,15 cm)/v auf eine Straße. Die Ausdehnung der Spirale in Strömungsrichtung 
durch eine Korrektur Ax zu berücksichtigen, war nicht nötig; der Nullpunkt 2=0 konnte 
einfach in der Mitte der Spirale angenommen werden. Bild 5 zeigt das Ergebnis für die Punkt- 
quelle. Die Meßpunkte für alle zund U streuen im ganzen Versuchsbereich um zwei Geraden mit 
gleicher Neigung; nur die beiden Punkte Uf,/v bei der kleinsten Geschwindigkeit und den kleinsten 
Entfernungen x weichen davon wesentlich ab. Für Ux/v > 5-10* gilt also für eine „„Punkt“- Quelle: 
Uf,/v = 0,625-10% (10-5 Ux/v)? und Uf,/v = 1,022-10% (10-5 Ux/v)?’, oder anders geschrieben: 


fi = 30 x Uxfv)!® und 5, —=4,752/Uc) BELA een. (8). 


| 
{ 
; 


EN Luft: ve: om ER 


Wasser: v = 10" En 


Bild. 6 Nomogramm für h und hr hinter einer Pinktaneile So für EN Hinter, einer Linienauelle, ee 


" merkenswert ist die Proportionalität Kon rund fsoodie een in einem Onkrachäller 
© x=const hinter einer Punktquelle sind also einander ähnlich für jedes 2 und Ine Gesehwin: BR 
7 BERN. Ma Ihre Gleichung kantet 22.7 7,7 x N IEEER 


Y 


TER NENEEE 3 (ylfa)® + (lei, 64fıP=in DR RR ENDE En Kol: 
ER Eon für 9%, = 0,2, 0,4, 0,6.und 0, Sbeia=ß= 1,8 sind in Bild 7 see N 
Die Ausbreitung geht danach vor allem nach der Seite von statten. Daß der seitliche turbulente 
Austausch stärker ist als der senkrecht zum Boden, ergab sich auch beieiner von To wnen d [7] 
untersuchten Strömung in. einem quadra- 
tischen Rohr; die in verschiedenen Wand- 
Bi . abständen gemessenen turbulenten Schwan- 
0... kungen parallel zur Wand und senkrecht zur 
©... Strömungsrichtung waren größer als die 
. senkrecht zur Wand. 
IR Dimensionsanalytisch ist anzuneh- 
=. men, daß die Formeln (8) auch für andere 
| Strömungsmedien als Luft gelten. Deshalb 
ist im Nomogramm’Bild 6 die Geschwindig- 
 keitsskala auch für Wasser angegeben, doch 
steht die experimentelle Bestaligun. hier- 
für noch aus. 


Baer, 0 Par die Linienquelle ergibt sich nach Bild 8 ein anderes Gesetz für /, als bei der Puhkt- BR 

quelle und zwar wird hier: Sie 

# RE dire H=0552/(02be5 ud h=®....2. 2.0.0). BER 
' Diese Formel ist auch in dem Nomogramm enthalten. Ein Vergleich mit der oben bereits an- 
gegebenen Gleichung (4) für 6, die sich nur durch den Zahlenfaktor von (10) unterscheidet, zeigt, 
‚daß die Dicke der Ausbreitungszone innerhalb der Reibungsschicht bei ebener Grundströmung 


23* 


5“ 2 ! 2m £ 


Bild 7. Isothermen in einem. Querschnitt y = const Er: 
hinter einer Punktauelle. 2 En 


R 105 10° 
Bild 8. Höhe f, des Temperaturfeldes hinter einer Linienquelle: j, = 0,55 #/(Ux/v):78,. 


Nachdem nun a, ß, fı und f, bestimmt sind, kann die maximale Temperaturdifferenz in 


. Abhängigkeit von. der Entfernung hinter der Quelle, 9, (x), aus der : Wärmeerhaltungs- 
gleichung (5) berechnet werden. Mit dem ®-Ansatz (7) und dem Potenzprofil für die Ge- 


© +9 
schwindigkeitsverteilung (3) ergibt sich für die Punktquelle aus i 0gc,dudydle] 
oo, oo ; : y=0 z=—o 2 

=2 f Sf 0o9g%,dudyde—= W = Ergiebigkeit der Quelle in cal/s: 09 6,0, («) 

y=0, 2-0 j 

n n+1, D., 

= ———tZ1U 2 [OR 

nn Bonn le 


Bei der Integration nach & ist dabei vorausgesetzt, daß’ die Temperaturunterschiede bereits 
innerhalb der Reibungsschicht völlig abklingen; denn das benutzte Potenzgesetz für « gilt nur 
solange y < ö, es muß also für y> du <U 9, gelten. Mit den Formeln (8) für f, und f, und 
n=717 und &—=1,3 wird 


09% do (x) = 0,0446 (Up) 7 (Ur) BUWr 2.2... (il) 


Ebenso wird für die Linienquelle wegen ji 090%, dudy= WB: 
0 


09,9, (x) = en ————— , W = Ergiebigkeit in cal/s längs der Breite B; und 
—. om ! ) 1/n . 
rn OR 
0990, = 2,92 (Uöjv)17 (Va MB BE BA 


Um noch den Einfluß der jeweiligen Turbulenzstruktur der Reibungsschicht im Experi- 
ment abzuschätzen, wurden einige Versuche mit künstlicher Turbulenz durchgeführt. Zu diesem 
Zweck wurde in die Druckkammer des Windkanals vor der Düse ein grobes Gitter angebracht, 
das aus 20 mm breiten waagerechten und senkrechten Blechstreifen bestand; die Maschen- 


BR N A 30020 0 10 zimm] 
Querprofile bei “= em RE NEN 'Querprofile bei © -25cm vi 
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Bild 9. Punktauelle: Temperaturprofile bei Bild: 10. Punktquelle: Temperaturprofile bei 
U = 13m/s, 0 =P=L8: ’ ' U =-30,9m/s; a =ß = 2,0. b 
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Bild 11. Linienquelle: Temperaturprofile bei Bild 12. Linienquelle: Temperaturprofile bei 
U=5b4mj/s; a =1,64,. ı U =18,0mjs, & = 1,77. 
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sein, daß f, kleiner und f, etwas größer wird, so daß die Ausbreitungszone noch flacher wird 
(vgl. Bild 5). Re * Kae: 


CN 
nd 


Für einige Beispiele sind die direkten Meßergebnisse in den Bildern 9 bis 12 aufgetragen; 


die gezeichneten Kurven sind nach den Gleichungen (7) bis (12) berechnet, Die Ab 1 


zwischen den Meßpunkten und den Kurven sind im allgemeinen erträglich, wenn man die ver- L 


hältnismäßig geringe Meßgenauigkeit bei der-Punktquelle berücksichtigt. Der Ansatz (7) sowie 
die aus den Messungen abgeleiteten Formeln für f} und f, stellen demnach diese turbulenten 
Ausbreitungsvorgänge genügend genau dar. N RR 

'4. Übertragung der Versuchsergebnisse auf die meteorologischen Verhältnisse. 

Dieses im Windkanal erhaltene Formelsystem scheint zunächst sofort übertragbar auf die 
Verhältnisse in der freien Atmosphäre, wenn dies auch eine Extrapolation auf wesentlich größere 
ReZahlen bedeutet. Leider ist das aber nicht ohne weiteres möglich, wie in [1] ausführlicher 
gezeigt wird. Ein Vergleich mit Messungen anderer Autoren über Ausbreitungsvorgänge im 
Bodenwind lehrte nämlich, daß die Turbulenz in der freien Atmosphäre offenbar wesentlich 
grobballiger ist, und zwar wird der turbulente Austausch hier rd. neunmal stärker als im Wind- 
kanal. Erst wenn man dies durch eine Änderung der Zahlenkonstanten berücksichtigt, kann 
man das obige Formelsystem auf meteorologische Anwendungsfälle übertragen. (So wird 


1 
Firm. 78-7" Fıran, ebenso für fs, und dementsprechend bei der Punktquelle e,,,, = &77 Kan‘) 


Zur Abschätzung des Geschwindigkeitsprofils in Bodennähe aus der Windgeschwindigkeit in 
z.B. 2m Höhe und der Bodenbeschaffenheit (Wasserfläche, Ackerland usw.) sind in [1] Dia- 
gramme angegeben. 


5. Vergleich mit der Theorie von B. Frost. 

Der Teil der bereits erwähnten Arbeit vonFrost, der sich auf unser Thema bezieht, läßt 
sich in der obigen Schreibweise folgendermaßen beschreiben: Für das Geschwindigkeitsprofil 
im Bodenwind wird ein Potenzgesetz angenommen: u —= U (y/ö)*!", wo ö groß gegen die sonst 
betrachteten Höhen ist und im Gegensatz zu den Windkanalversuchen als konstant angesehen 
werden kann; aus der Kontinuitätsgleichung folgt damit » = 0. Ferner wird der Austausch von 
Beimengungen &y oder von Wärme &, gleichgesetzt dem Impulsaustausch &,, der hier unab- 
hängig von # und nur proportional y1-t'* ist. Dann vereinfacht sich Gl. (6) für die Linienquelle 

E 1m 2% — ER 
zu :yYN—M Fp 
y— © und dem Erhaltungssatz f wddy=const. (Wenn man ey nicht gleich sondern nur 

ö 


0 
(van =) mit den Randbedingungen d—P0 für —0 y #0, 9 —0 für 


proportional e, setzte, würde sich nur die unwesentliche Konstante a, zahlenmäßig ändern.) 
Mit dem Separationsansatz = d, (x) F(H), wo H= y/f(x) ergibt sich speziell bei n — 7: 


RE 1+1/n 
14+2/n -alr2in irn 
I» x e Foe und dmx 
20,78 I e-dl:20 — 20,89 ° 


statt nach den Windkanalversuchen 


8 
ER, F=e-#" und a 
Der Vergleich der einzelnen Exponenten zeigt, daß zwar die Potenzen von & — also der 
Koordinate in Strömungsrichtung — theoretisch genügend genau vorausgesagt werden, das 


Abklingen der Temperatur bzw. der Konzentration nach der Höhe dagegen nicht; 
(&rneor. = 1,29, &versucn = 1,8). Dieser Unterschied zwischen der theoretischen Modellströmung 
und der im Kanal könnte vielleicht damit erklärt werden, daß bei den Versuchen ö, + const 
war, und zwar stieg die Reibungsschichtdicke innerhalb des Meßbereichs um rd. 20% an. Es 
müßte also die Theorie hierfür folgendermaßen verfeinert werden: u = U (y/ö)”, 

m BE m 2 n 
I U ö(yjö)”t! und = m+n)Omtn U ölylöoy ml —(yjö)!+2m], 
Wie die weitere Rechnung zeigt, wird dann aber. der Separationsansatz 9—= 9, (x) F(y/f (x)) 
nicht mehr erfüllt, da ö nicht als Potenz von x (Entfernung von der Wärmequelle), sondern nur 
von X (Rücklage vom Ursprungsort der Reibungsschicht) angenähert werden kann. Ferner 
versagt im Fall der Punktquelle ein entsprechender Separationsansatz auch schon bei » = 0 
ö = const und e (y), wenn nicht statt dessen e3 — const (also ey + e,) gesetzt wird. Da somit 


® 


N rerseits könnte statt wie, Brose mit. einer he über ‚den Wärme- a 2: 
Stofaustausch & die Ausbreitung zu berechnen, jetzt umgedreht aus den Versuchsergebnissen | 
'&, nach GI. (6) berechnen, allerdings wieder nur unter der Vereinfachung ö= const, da sonst 
ee ö Ri nn DE Unahhängigen En a Yy auch ‚noch X in ‚dep, en auftritt. ‘Man ] 
SEHR RO N 2: Re 
en Es ur EHER &9 = an et a yrn- Due 0, 15 v( (VOpyo2 ( Uxfoyoa ‚d4 (yio) ‚34 a EEE Ne J 


statt 89 = &; = «0 y%86 hei Frost. Diss Ergebnis ist aber auch aeenigende a sta 
am Rande der Ausbreitungszone zu verschwinden wächst e3 monoton mit der Höhe an. Es ist? 
Ei danach überhaupt fraglich, ob — abgesehen von der experimentell nicht erfüllten Annahme Ra 
 d= const — die festgestellten Ähnlichkeitsbeziehungen (= 9, (x) F(H )) ausreichend genau 
gelten, um tieferliegende Folgerungen bezüglich der Turbulenzstruktur zuzulassen. Dazu 
"wären weitere, möglichst genauere Messungen nötig, z.B. bei verschiedenen Rücklagen der ’ 
Wärmequelle (X) um den Einfluß von: ö (X) zu erfassen, oder auch Messungen des Wärme- 
f  machlaufs, wenn die Quelle sich nicht am Boden sondern in verschiedenen Höhen innerhalb 
der Reibungsschicht befindet. Vorläufig kann also das aus dem Versuchsmaterial gewonnene - 
. Formelsystem nur als praktisch EeRUgeuNE Näherung ohne thoeretische F Bee ange- 
© ee En x 
Be: ® Schrifttum. 
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Über Rayleigh-Wellen an Ei; Grenzfläche zweier Medien. ® 
Von H. Koppe, Max-Planck-Institut für Physik in Göttingen. or 


ER Die Ra yleighsche Theorie der Oberflächenwellen wird auf die Verhältnisse an der Grenzfläche zweier BE. 
Körper resp. eines Körpers gegen eine Flüssigkeit verallgemeinert. Im ersten Fall ergibt sich, ‘daß nur unier SLR 
bestimmten Bedingungen Grenzflächenwellen entstehen können, die dann immer schneller laufen als die Ray- > 
leigh-Wellen im Körper größerer akustischer Dichte. Im zweiten Fall gibt es immer Grenzflächen, die stets g 
langsamer laufen als die Rayleigh-Wellen ın der fesien Komponente. 
Rayleigh’s theory of surface waves is generalized to the conditions at the interface between iwo ide 
respeclively of a solid and a liquid. In the, first case, it follows that interfacial waves can arise only on 
certain conditions. T'he velocity of the waves is then always larger than that of ihe Rayleigh waves in ihe \ 
body with the larger acustical- density. In the second case, there are always interfacial waves Ihat are 
propagating more slowly than the Rayleigh waves in the solid component. 
D’auteur donne une generalisation de la theorie de Rayleigh sur les ondes de surface pour les cas solide- 33, 
solide.et solide-ligquide. Dans le premier cas les ondes de surface limite ne peuvent se former qu’& certaines: 
conditions. La vitesse de ces ondes est plus grande que celles des ondes de Rayleigh dans le solide de plus 
grande densite acoustique. Dans l6 deuxieme cas ül y a des ondes de surface limite qui se propagent plus len-' 
iement que les ondes de Rayleigh dans le solide. 
Teopua Penea HOBePXHOCTHEX BONH ODoÖMaeTcH Ha cHyyai TPaHuyHOH TOBePXHOCTU 
MEKAy ABYMA TBePAEIMH TEIaMH MIN MEIKILy TBEPBIM TEIIOM H 3KUTKOCTBIO. B mepsom cıyyuae 
- B TPaHHYHOÄ TOBePXHOCTH BO3HUKAIT BONHEI TOAIBKO HPH OCOÖBIX YCHOBHAX. ÜKOPOCTB BOAHEL 
B3TOM CiIyyae BCerxa 6onbITe cKOPpocTH Bonn Peren B Tee c 6onpmeii arycTuygeckof IMOT- 
HOCTbIP. BO BTOPOM cAy4ae TOBePXHOCTHBIE BOJIHBI BO3SHHKAMT IIPH MOÖBIX YCHIOBUAX, BPaSSM 
OHH BCETA& N MeAIeHHEee, YeM BOJIHBI Benen B TBEPAOM, TeıTe. \ 


1. Rayleigh-Wellen. 


= Von Rayleigh!) stammt die Entdeckung, daß sich längs der Oberfläche fester Körper 
B Oberflächenwellen fortpflanzen können, die nur exponentiell abfallend in die Tiefe dringen. 
Br In der Folge ist dieses Problem hauptsächlich nach geophysikalischen Gesichtspunkten behan- 


ENTE NA IOTIN ONONRTE N TERREN ARTS TE  NORTTELST U DNER. N - Vt 


1) Rayleigh, London Math. Soc.17 (1887), vgl. auch A.E.H. Love, Den 2. Aufl. (Cambridge 1906), 
S. 295 und. die folgenden Fußnoten. 
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delt worden). Über Rayleigh hinausgehend, wurden dabei auch g 5 ete Medien 
betrachtet ®), bei denen man eine für die Geophysik wichtige Dispersion u  e 


. 


Hier wird ein anderer Fall untersucht. Es sollte offenbar auch an der Grenzfläche zweier 5 j 
Medien Oberflächenwellen geben, die von der Oberfläche aus exponentiell nach beiden ‚Seiten | 
abfallen. Man erhält sie durch eine einfache Verallgemeinerung der ursprünglichen Rayleighöcken 


Ableitung, die zu diesem Zweck noch einmal kurz wiederholt sei. 


Es sei zunächst ein als isotrop vorausgesetztes Medium gegeben, welches den Halbraum n ei ue 
z2>0 erfüllt und durch die Dichte o und die elastischen Konstanten A und u gekennzeichnet ist. DL: 


Dann gilt die Bewegungsgleichung für den Deformationsvektor u: = 
i u 


err7 — (A + 2u)graddivu—prot ot u, ....2..n. M. 
Der Spannungstensor ist gegeben durch die Komponenten 
, ou 
7, =Adivu +72 \ 
0) (] 
nl tn) ER Ne elle (2) 


und 3 analog gebaute Komponenten, die im folgenden nicht gebraucht werden. 
Wir untersuchen eine Welle von der Kreisfrequenz p und machen demzufolge den Ansatz 


vu = U, e-ipt TE ERTL (3) 


Der Index „ kann unterdrückt werden, da keine Verwechslungen zu befürchten sind. Dann 
wird aus (1) 


(A+2u) graddivu—urotret utoP@u=0 ...ı..... (1’) 
und diese Gleichung ist zu lösen unter der Randbedingung 
Ta FREE RETTET (4). 
Wir setzen für u: : 
Peg WERTE BER (5,1). 
dv. me er a a ED 
Dann genügt p der Gleichung 
29 — BE ENBL DBIS VERS u Ri 6 
Ay+%Rpoy=0, Äh rau? (6) 
und w außer der Nebenbedingung (5,2) der Gleichung 
An+kw=0, Be RE BER (7). 


Die Lösungen sind ebene Wellen. Wir legen die x-Achse in die Ausbreitungsrichtung und können 
dann setzen / 
= Aexp(—rz+ima); R=m?—R (8) 
w=(sB, C,imB) exp (—sz+ima); 2=m?— 2 
Dabei sind bei ıw die drei Komponenten in der Reihenfolge x, y, z in der Klammer zusammen- 
gefaßt und gleich so durch zwei Integrationskonstanten B und © ausgedrückt, daß die Be- 
dingung (5,1) erfüllt ist. Da die Frequenz vorgegeben ist, ist auch der Ausbreitungsvektor m 
oder, was auf dasselbe hinauskommt, die Schallgeschwindigkeit 


N ee BERLIN €; 


als Unbekannte zu betrachten. Um die vier Größen A, B, CO, m zu bestimmen, berechnen wir 
aus (8) und (5) die Komponenten des Spannungstensors nach (2) und setzen in die Randbedingung 
(4) ein. Das liefert 3 lineare Gleichungen für die vier Größen A, B, ©: 
ul2m?— kB) A—2imsB]=0 | 
ull2imr A + (2m? — k2) B]= 0 J 
usC=0 
") Galitzin, Seismometrie (1914). Uller, Ann. d. Phys. 56, 464 (1918); Gerl. Beitr. Geophys. 20, 410 (1928). 


Hardtwig, Z. Geophysik 18, 164 (1943). 
’) K.Sezawa, Bull. Earthquake Res. Inst. Tokyo 3, 1 (1927); 5, 85 (1928). 


... (10). 


> y , A y N 
J rhıchtete 1.4 


2 = K 


| =) 


nen 1; 


ES dh. a 


R | ee 
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Die in Frage kommende Wurzel. ist 


BER EHE = 0,8452933; Hr — 0,919402 . 


he = as Bednnke Ergebnis, daß die Geschwindigkeit der Ra er i en- Wellen um des X 
Faktor v, 92 kleiner ist, als die der transversalen Men 


Kok 


. 4 
iR: ki 


%; Oberflächenwellen an der Grenzfläche zweier Tester Körper. 


E 2 7 "Um das Analogon der Rayleigh-Wellen für die Grenze zweier Medien zu erhalten,er- 
< - füllen wir den Halbraum z > 0 mit dem Medium 1, und den Halbraum z<{ 0 mit dem Medium 2 
2,5% aa unterscheiden die entsprechenden Größen durch Indizes. Diesesind überflüssigbeipundm, 
die für beide Halbräume den gleichen Wert haben müssen. C, und C', sind wieder gleich Null. 
Für den Halbraum 1 gelten die Lösungen (8). Für den Halbraum 2< N) muß man beachten, daß 
‚jetzt die Wellen für > — » abklingen müssen. Man muß also enweder für r, und s, negative 
Werte wählen, oder, "was begemer scheint, die beiden Größen ebenfalls als positiv ansehen, 
‘. dann aber überall vor r, und s, das Vorzeichen umkehren, so daß (8) z.B. lautet 


EL 
m = (— 82B;, 0, im B,) exp (+38, z +29 2). SH 
ö Als Nebenbedingungen sind jetzt zu fordern: “ 
| 1. die Gleichheit der Spannungskomponenten für z= 0: | FEB 
ne Tan Ta ee Fre rer So TEO 
2. die Stetigkeit der Verschiebungskomponenten ” 

u == un F u =: u” A A 3 


Dabei ist stillschweigend die (wohl praktisch immer zutreffende) Voraussetzung gemacht worden, 


ns 
DE EN 


F::: daß die beiden Körper längs der Ebene z = 0 miteinander ‚‚verschweißt‘ sind, und keine Glei- Be 
E- tung möglich ist. 

E: Er» Die beiden Bedingungen (13) kann man nach (10) sofort anschreiben: 

E ne ; vll? m’ —k,)) AR —2ims, Bl =u,[(2 m’ — k,) A, +2ims,B,) | i (15) : 

E ER ei ‚uPimr, 4, + (2m — ki) B)] = url—2imr,A,+ (2m? —k,) BP | 
# Entsprechend erhält man aus (8) für (14) 

B.. im A, +8, B,=imA,— s,B,[= X] | (16) 


—1r]A,+imB,= 1,4, +imB, [= Z] ] 


Das sind wieder 4 lineare Gleichungen für A, bis B,; damit sie eine Lösung haben, muß die 
Determinante verschwinden. Das liefert eine Gleichung für m oder x. Bei x muß natürlich 
definiert werden, ob es mit k} oder k) gemeint ist. Wir setzen fest, daß vr, < vr, und definieren 


E | 2 Ve 


Dabei bedeutet v4—= vr, Yx die gesuchte Geschwindigkeit der Grenzflächenwelle. - 


Ko ppe, Über Rayleigh-Wellen an der Gronztläche zweier Med 
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Die Ermittelung von x ist nur durch numerische 


homogen sind, bekommen wir wieder die entsprechenden Gleichungen in x durch die Sub- 


stitution m = 1, 3 


S | 
Bu; B=a; Mags Magy nr...) 


op’ 


Wir können zunächst einmal (16) nach den dort bereits eingeführten Hilfsvariablen X und Z. 


Verschiebungsamplituden in der x- und z-Richtung für z= 0) auflösen und erhalten 


D 8 D ER 
DR RE Far ee a ie | 
; ai a: N a in (19) 
r Kar Tg 
Ber ETgR 


mit 


De Meg, 


..s=lfl—al; A, =1—1r,8 


| BR REBNN 
| „= 1-5] Fe a Ir 


Setzt man das in (15) ein, so erhält man für X und Z das homogene Gleichungssystem 


ira -ur(2)|x- [se +8 (Z)|2=0 


Ne 


ar Tan (21). 
r(&)|x ira ur (2)] 2-0 
R(&) +4 os Een v5 | 
Tafell. | Dabei sind F, R und $ die Funktionen: 
j i i 22.1—3) für q = 1/3. t 
Hilfsfunktionen nach Gleichung ( ) für q / F=2 = —_ (a), 
® Fe) sa | R@ Yya—ı) (1—g%) 
a s)=[2 —Fi)]yl—t . .! (22,2), 
0.00 0.50000 . 1.50 000 1, h 
0.05 0.50 658 1.45 561 1.48 092 - 2ü=[2 — Fijyl —gt. . (22,3). 
Br u en, ei; 0 u in Sie hängen also an sich noch von g ab. (Wenn 
020 | 0.523844 1.31.621 1.42 166 9 #g,, dann ist q, resp. q, zu nehmen, je 
0.25 0.53 665 1.26 731 1.40 106 nachdem ob das Argument den Wert x oder 
0.30 0.54 567 1 21678 1.37 969 x|vy hat!) Weiterhin wollen wir, um einiger- 
> m Den an > 2ER 2 maßen übersichtliche Verhältnisse zu bekom- 
0.45 0.57 713 1.05 523 131 181 men, wieder annehmen, daß beide Körper die 
0.50 0.58 957 ‚99 733 1.28 754 Cauchysche Bedingung erfüllen und demnach 
Den N er pi a, 1,40 a 41 = 9, = 1/3ist. Für diesen Wert können die 
0.65 0.63 558 30720 120750 Hilfsfunktionen 2 S, R aus der Tabelle 1 
0.70 0.654938 |  .73672 1.17 774 entnommen werden. ; 
0.75 0.67 721 .66 140 1.14 558 Damit das Gleichungssystem (21) eine 
0.80 0.70 347 ‚57 982 1.11 028 Lösung hat, muß die Determinante von (21) 
BR Kun = 1: Dr 1el ER verschwinden; das liefert für x die Gleichung 
0.92 0.79 659 .34 037 1.00 204 Eh 
0.9 | 0.82.061 .28 889 .97 731 Fa—ur)| = 
0.96 0.85 041 .22 992 ‚94 797 7 (23). 
0.98 0.89 136 .13 678 .90 972 \ % x 
1.00 1.00 000 00.000 „81650 Er Is (2) +usS (3) E (2) Hu R 3) 


Setzt man hierin zunächst u=0, so muß man wieder zu den in Abschnitt 1 behandelten Ra y- 
leigh- Wellen an einer freien Oberfläche kommen. Tatsächlich läßt sich durch eine einfache 
Umformung zeigen, daß die dann aus (23) entstehende ‚Gleichung F?— RS=0 mit (12,2) 


wo. 
R 
Ar ER 
ru = 
ur: 


a 


% 


ua... FE Ze 


Sr " 


Ar, 


 . merkwürdig, 


len zu tun, sondern die transversale 


- Komponente im akustisch dichteren Mit- 
tel klingt nicht mehr ab. Nähert man 
sich nun auf der Kurve = -1 dem Punkte 
vn =1, so geht einmalauch $, nach Null, 
. die transversale Komponente klingt also 


BRSEN allg en De ck her, die e Werte von win BI 
a: 3) gibt Bild 10 ‚Dabei fällt zunächst nl auf, daß es durg! aus nicht unter ‚ganz, be, R 
er ‚liebigen ‚Bedingungen ‚Grenzflächenwel- RAN 

len geben kann: Bei allen Wertepaaren 
und v%, die in das schraffierte, von. 
der Kurve =1 umschlossene Gebiet 
; ‚fallen, sind keine Grenzflächenwellen 
"mehr möglich. "Dabei scheint es zunächst 
daß gerade der Punkt 
Se 1,0 — 1,\ bei dem also die „Grenz: 
fläche“ nur noch gedacht ist, im erlaub- 
ten Gebiet liegt. Dieses Paradoxon klärt 
 sichso auf: Bei der Annäherung an die » 
— Grenzlinie &—= 1 wird zunächst einmal 
8, =1. Man hat es also schon hier nicht 
„mehr mit eigentlichen Grenzflächenwel- 


ei von uw und ı & 


N OR GE = RL 


Bild 1. Geschwindigkeit vg der Ürenztlächenwells an der I s R 


zweier elastischer. Medien, die den Caue Kiyschen Bedingungen 
genügen. un Ha = Torsionsmoduln, vm, 2 vom, = Geschwindigkeit 


‚der transversalen Schallwellen, (Der Index 1 ist, so zuzuordnen, 4 


Va \2 
daß Yn A) 7, ’% Aufgetragen sind die Kurven 2 a 


an des schraffierten Gebietes sind keine re > 


Um. 2 


möglich, (Die Grenzkurve <= 1 hatim Punkte Ki or es a 1. BR N. 
KH, "% at Aal 
eine sehr flache Spitze, ae sich i in der FL nicht, darztellen, TR 


1ä8t.) 


=const. Ki 


auch im zweiten Medium immer schwä- 


cher ab. Außerdem gehen aber, wie man mittels der Gleichungen (21) ad: (19) nochrenfiuent Air 
kann, die Amplituden A, und 4, der longitudinalen Anteile nach Null. en artet im Grenz- 


falle u —1, vp—>1 die Gen ztachenweile in eine ebene Transversalwelle aus, die in Richtung 


Ki gedachten ‚‚Grenzfläche‘“ verläuft. 


3. Grenzilächenwellen an der Grenze zwischen festem Körper und einer Flüssigkeit. 


Man kann auch die Grenzflächenwellen an der Grenze zwischen einem festen Körper und 
einer Flüssigkeit untersuchen. Dabei sind zwei Vereinfachungen von vornherein naheliegend: 


1. wird man die Flüssigkeit als reibungslos ansehen, denn die Berücksichtigung der inneren 


Reibung würde zu einer Dämpfung führen, die bei den bisherigen Rechnungen ebenfalls (wenn 


auch vielleicht mit mehr Recht !) vernachlässigt worden ist; 2. kann man sich auf kleine Schwin- 


gungen beschränken und deshalb in der Bewegungsgleichung der Flüssigkeit die nichtlinearen 
Trägheitsglieder weglassen. Damit werden aber die Bewegungsgleichungen einer Flüssigkeit 
identisch mit denen eines festen ne für den u, = Oist, während A jetzt mit der Kompressibili- 
tät x eg 


EN RN WR Na CHOR) E 


4 
zusammenhängt. 


Geben wir dem festen Körper den Index 1 und der Flüssigkeit den Index 2, so haben wir . 


für u, wieder die Lösung ( ) und (8), während u, gegeben ist (Vorzeichenumkehr wie in Ab- 
schnitt 2). VER 
N — gradp, = grad [Azexp(rgz2+im®@)] ».- » - ....2..0.(28). 
Lediglich die Definition von r, ist etwas abzuändern, da es jetzt keine Transversalwellen gibt 
und damit auch der Koeffizient g, sinnlos wird. Wir gehen auf die durch (6) und (8) gegebene 
ursprüngliche Bedeutung zurück und setzen 


nm NEL A eh TORE, 


Anders werden jetzt die Randbedingungen. Zunächst bleibt bestehen 


2 i 
NT RN N BE ar ie: - 


Bas ES en a EEE RSG). 


Da aber eine reibungsfreie Flüssigkeit an dem festen Körper nicht haftet, fällt die andere der 


va Di ’ u d 
zu (27) un 


gerechnet gibt di 


2imnAytam—MB-0, 


x R3 EN { en ur 
Für die Determinante ergibt ich BER. © Far REIT 


Graphische Diskussion der Gi 


.(31). gehörigen Wert 0.845. Man erkennt erstens, dßesandr 


-  wellen gibt, und zweitens, daß diese, im Gegensatz zur Grenzfläche fest-fest, immer lang- RER 
'samer laufen als die Rayleigh-Wellen. RE 
Eingegangen: 29. 5.48. 
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Pfeil-Diagramme für Zylinderfunktionen.. 
Von Fritz Emde in Pretzfeld (Oberfranken). | | e 


d (28 


3A —2ime, Bı +72 MA2=0, 


RL 


rd, Hi mB, BL te eh N 


' En 


RES MEET N A N 
stehen. Für 9,0, d.h. n=0 geht (32) wieder in 2,1) OR 
“über. RED, RE ER 2 ne 
Zur Diskussion von (31) ist in Bild 2 die linke undde 


rechte Seite von (31) als Funktion von = aufgetragen, nd 
. zwar die letztere einmal für „=.0 und einmal für „>0. 4 
"Die Wurzel x ist dann durch die entsprechenden Schnitt- 
punkte gegeben. xzentspricht dem zu den Rayleigh-Welen 


Grenzfläche fest-flüssig unter allen Umständen Grenzflähen- 


at 


Diagramme für verschiedene Zulinderfunktionen desselben Arguments mit gleichen oder entgegengesetzt 
gleichen Parametern. Polardiagramm für eine Besselsche oder Hankelsche Funktion festen rein 


imaginären Arguments mit stetig veränderlichem Parameter. Umlaufs-Diagramme, 


% 


Diagrams are stated for the different Bes sel functions of the same argument and of equal or ne- 


gatively equal parameters, and a 
_ purely imaginary argument am 


olar diagram is given fra Bessel ora Hankel function ofafixed - 
of a continuously variable parameter, further circulation diagrams. x 


L’auteur donne des diagrammes de differentes fonctions de Bessel du mime argument dont les 
parametres sont egaux ou negativement Egaux. Em suite il donne un diagramme polaire pour une fonction de 
Bessel.ou de Hankel d’un argument purement imaginaire et d’un paramötre continüment variable, 

B aToii pa6ore gamrca nuarpamMH I PASIHMYHBIX IMIMHAPHYSCKHX PYHRUMÜ OAHOTO 
MH TOTO-K® APTyMEeHTA, HOPAAEH (TAPAMETPH) KOTOPBIX PABHEI IIO Be/IMYNHe HM PABHBI IH 
oÖparust 10 3Haky. Kpome Toro namrca nNoNapHse Auarpammsuı Beccetesux byHurumk u 
hyarıml Take sta IOCTOAHHOTO KOMIMERCHOTO APTYMeHTA C HEIPEPEIBHO-TEePeMEHHEIM 


TapaMmeTpoM. 


Von den Zylinderfunktionen, die man eingeführt hat, sind nur ‚zwei (beliebig wählbare) 
voneinander linear unabhängig, weil die Zylinderfunktionen Lösungen einer linearen Differential- 
gleichung zweiter Ordnung sind. Die Zylinderfunktionen sind im allgemeinen komplexe Größen. 
Man kann sie daher nach Wesselund Argand durch Pfeile in der Zeichenebene darstellen. 
Wenn man das tut, werden die Beziehungen zwischen verschiedenen Zylinderfunktionen sehr 
übersichtlich. Seltsamerweise scheint das bisher nicht geschehen zu sein. Die Formel sagt uns, 
wie wir in jedem Fall zu verfahren haben, und läßt uns das Ergebnis nur roh vermuten. Das 
Diagramm zeigt uns das Ergebnis, freilich'nur für einen einzelnen Fall, und läßt uns auch er- 
kennen, wie wir zu verfahren haben. Die Diagramme werden beim Studium und beim Gebrauch 
der Zylinderfunktionen nützlich sein. 

Es sollen hier einige solche Pfeil-Diagramme vorgeführt werden. Dabei werden wir den 
Parameter als reell voraussetzen, was fast immer genügt. Wir werden vor allem zwei Fälle ins 
Auge fassen: 1. Das Argument ist reell und positiv. 2. Das Argument ist rein imaginär und 
positiv. Die Buchstaben 7, r, ®, y sollen im folgenden reelle positive Zahlen bedeuten. 


er 


—— 


© Ord i Jate ı der, Spitze des Pfeil 
Ds r a Ir e 


Für 2 Er 68 ist en nt => SU 58 und N... us = 0; 20 (Bild N). Der Beil u 1/3 
hat die Länge = - 0,61 und den. ‚Richtungswinkel n=aretg Nualdus) = 0, 212 > N 
° Wenn wir den Pfeil um 180°:3,— - 60° vordrehen, geht, er über in den Pfeil Hs mit: de 
0, & und EN le ae 8 Ra Ren vn: Seine ee hat, ae Absziss 


“ SON R 
Zu und klein Ne Ru wa und. ER Wahlen aa Be cisp 3 


RE Pe 


m (0.20) 


BE NENE ee Bild 2. 
Pur p=1/3 (mod 2). Reelles Parameter 9 =1/3 (mod 2). Reelles Parameter p gleich dem Argument 
Argument © =1,68. Gegeben J,,, Argument » = 0,20. Gegeben J,,3 (reell und positiv). \ 

und N und Nyjs- 


13° 


in Lange 5.2 1,60 und’ den ‚Richtunaswinkel. 7 = arete 5 2a. 7a 


Drehen wir ihn um 60° vor, so entsteht der Pfeil AHV 3; mit der Länge = 1,60 und dem Rich- 
tungswinkel — 71,2° + 60° = — 11,2°, Die Pfeilspitze hat die Abszisse J_, 7, = + 1,57 und 


die Ordinate N_,/; = — 0,31. je 
5 Man kann ohne weiteres den Pfeil Br hinzuzeichnen, indem man den Pfeil HN an der. 
= reellen Achse spiegelt. Auch dazu ein Beispiel. a ist YpJ,(p = +0,44 bei kleinem 
 pbis p=3 und = 0,45 bei großem p, und Yp N,(p) = — 0,78 bei rel p bis p = 15 und 
— — 0,77 bei großem p (siehe die nen bei er x h son, Bessel functions, auf-S. 746) und 
daher /pHY) (p) = i2"° 0,89 and Vp By) (p) = i*”'%0,89. Für En —8 erhalten wir J, (8) 
FEN, = und N, (8) = N_, (8) = — 0,39, daher H{) (8) = H\ (8) = i72!? 0,45. (Bild 3.) 
Wenn J_, und N, ER sind, so wird den Pfeil 2) um den Winkel p2 
rechts herum drehen und erhält dann den Pfeil Hy) und darauf J, und N,, indem man ihn 
nach der Drehung auf die Achsen projiziert. 

ER Statt den Pfeil zu drehen, kann man die Achsen im entgegengesetzten Sinne um einen 
gleich großen Winkel drehen. So wird man verfahren, wenn von den beiden gegebenen Zylinder- 
funktionen die eine einen- positiven Parameter hat, die andere einen ebenso großen negativen 
Parameter. Für = 0,2 ist J,,a5 = + 0,62 und J.,,35= + 1,43. Das Achsenkreuz wird um 
0,5 = 45° rechts herum gedreht. J,,.; wird waagerecht nach rechts aufgetragen, und J_y,a5 


Eee unter 0,5 Ah Pches unten (Bild A). In den ee werden A h 
- punkt der beiden Lote hat vom Ursprung die Entfernung h = 1,54. ‚Der P 
nach dem Schnittpunkt ist gegen die ursprüngliche positive reelle Achse Be N 
000 == — 66° geneigt, also nach rechts unten. Die Ordinate des Schnittpunkts ist im alten Achsen Ba: 
Be. kreüz REN 4 41 und im neuen Achsenkreuz None = 0,56. Auch im neuen Achsenkreuz ERS 
Be ist der Pfeil H/ nach rechts unten = 2 
geneigt, aber nur noch um den Win- 
kel — 0,74 +0,50” = — 0,24. oder 
68° 445° = 21° N ae 0 
Für 20,4 ist Jon = +0,32 
und J_0, =:+978 (Bild ö). Das 
Achsenkreuz, wird um -1,5° = 135° 
rechts herum gedreht. J,, wird 
waagerecht nach rechts aufgetragen 
und J_9,,, unter 1,5" nachlinks unten. 
Der Schnittpunkt der in den End- 
punkten errichteten Lote hat vom 
| u Ursprung den Abstand h= 1,45. 
4721020) # (040) Der Pfeil Hi}; vom Ursprung zum 
Ren u es Bei Be die reelle 
“r 3 AN chse des alten Achsenkreuzes um 
mente = 0,80. Gegatan Tops“ Mesa Arsument #040. 0,86 — 77” Eeneige Soc | 
und J_9,25. Gegeben Jg,75 und J_0,75. rechts unten. Im neuen Achsenkreuz “TB 
hat er die Richtung — 0,86” +1,50” —= 
+ 0,647 oder — 77° + 135° = +58°, Die Ordinate des Schnittpunktes ist im alten Achsen- 
.. kreuz N,,; = — 1,42 und im neuen Achsenkreuz N_,,,;, = +1,23. | 
- Für zwei Parameter p, und p,, die sich um eine gerade Zahl unterscheiden, p,= p, (mod 2), { 
bekommt man gleiche Drehwinkel, z.B. für p = 1/3, 7/3, —5/3, nämlich 60°, | 
Wir haben den Diagrammen ganz bestimmte Zahlen zugrunde gelegt. Bei einem andern | 
Maßstab gelten diese Diagramme für andere Zahlen, wenn diese nur in demselben Verhältnis 
zueinander stehen wie die von uns angenommenen Zahlen. Bei reellem festem Parameter +p 
oder — p und wachsendem reellem x schwingen nun die Zylinderfunktionen J und N dauernd 
um Null. Also wird jedes beliebige reelle Verhältnis J : N unendlich oft angenommen werden. 
Daher wird auch jedes beliebige Diagramm, wenn es nur richtig konstruiert ist, unendlich oft 
gültig werden. 
Die bisher vorgeführten Bilder gelten für reelles positives Argument x. Von den Zylinder- 
funktionen mit reellem negativem Argument — x werden wir später (unter III) sprechen. 
In unsern Figuren haben besonders die Beziehungen zwischen den Abszissen J, und J_, 
und der Ordinate N, praktische Bedeutung. Die Ordinate N_, pflegt man nicht zu benutzen. 


obgleich dem nichts im Wege stände, Die Pfeile H\) und ya dienten uns nur als bequemes 
Konstruktionsmittel. Denn diesen komplexen Werten kommt keine SCHRERUSIES Bedeutung zu... 


II. Rein imaginäres Argument. 
Setzt man 
’ Jliy)=W@%,(y) und H\) (iy) = ir- HAN: Er a a 
so sind bei reellem positivem p und y auch \} und $ reell und positiv. !) 
Wenn y wächst, wird \} immer größer und $ immer kleiner. Bei veränderlichem 
y haben daher die Pfeile J„und H\) eine feste Richtung. Bei wach- 
sendem p dreht sich der Pfeil J,links herum und der Pfeil HP rechts herum. Wenn p=—0,5 
(mod 2) ist, haben die ‘beiden Pfeile dieselbe Richtung, und wenn p = —+0,5 (mod 2) ist, ent- 


gegengesetzte Richtung. Bei 9 = 0 weist der Pfeil J, nach rechts, der Pfeil H} nach unten. 


Ferner: ist 
RER | 
H=WrH, = ng. ER re ee 
a a Ru ihn a WE ec ! 
% Jy=V"PJ,+iH, sinp 2 =ır7 N, HERD FLIRT ER), 
2) wie bei den Hyperbelfunktionen benutzen wir Fraktur-Buchstaben und bezeichnen die Beträge der 


‚Zylinderfunktionen rein imaginären Arguments mit = und 9. In England pflegt man diese Funktionen „modified Bessel- 
functions“ go nennen und I statt 9 und K statt = 9 zu schreiben. Hier wäre es sehr störend gewesen, statt 9 
beständig EN K schreiben zu müssen. 


esch jes 
Pr sich unterscheiden, läßt sich leicht nach Größe und Richtung konstruieren 


Pe ‚(Bild 2 Man lasse den Pfeil 3 vom Ursprung ausgehen und füge den 3 
Bi. Pieiyi2, Dan. Das Lot vom N auf ne ist die er Strecke Ba ro 
E sinn. ? u 
Bild 6. 


2 0 an x. RE re 


a r : h 'tion von Hi 
vi ind TE i BR DR HN 4 N 
u ES Ir % RR op ir ER an ne EP: 


‚Dazu werden wir den Parameter p weder als ganze Zahl, noch als Hälfte einer ungeraden Zahl 
"annehmen, denn dann entarten die Zylinderfunktionen. Um günstige Lagenbeziehungen zu er- _ 
halten, empfiehlt es sich, den Parameter klein zu wählen. Wir entscheiden uns für v=18 


° positiven gemeinsamen ‚Faktor sei Js =4 +i8 mit dem Betrage 40,8 und HN = = —4 
IM ı mit dem Betrage 20,4. (Genauer ist dann arc tg 0,2 = 0,125 666° = 11 ‚3099° und... 


' wärts ein wenig nach links (Bild 7). Die Differenz ergibt HN — Js =i Ns = — M—i38, 


“tionen mit dem Parameter —p= — 1/8. 


-+77—i1,5. Das ist zu dem Beer 


teil J_ enig länger als das em von Er an der 
GE chtung wie dieses Spiegelbild. Die kurze Strecke, um ‚die 


Aus I und aD foißen schließlich Parameter p=1/8 


| | DL ‚— a0) N 
"Wir wollen jetzt die Konstruktion des Diagramms für rein imaginäres Argument erlänbekn: A 


=.0,125. Der Drehwinkel beträgt dann 0,125” = 11,25°. Abgesehen von einem unbekannten 


p=(, 125 666.) Der Pfeil J|,;s weist nach rechts ein None, aufwärts und der Pfeil Tan 


also N,),, = — 28 5 44 mit dem Bene 51,2. Der Pfeil N,,, weist nach links oben (Bild 7. 
Ferner ergibt sich HN = a His —= 84 +:36 mit dem Betrage 91, 4, weist also nach _ Ai, 


_ rechts aufwärts. Die Pfeilspitze von Js liegt auf der Mitte der Verbindungslinie der Pfeil- $ = 


spitzen von Hs @) und Hy 


"Wir kommen jetzt zu den Funk- 


Zunächst liegt der Pfeil A; an der 
imaginären Achse spiegelbildlich zu dem 
Pfeil HU. Der Pfeil i2? Js 40—i8 
ist das Spiegelbild von J,/san der reellen 
Achse. Wir brauchen jetzt die Werte von 
sin 0,25” = 0,383 und cos 0,25" = 0,924. 


Hiermit ergibt sich i H')k sin 0,25” 


x . un, N, er 
u nd Ve re 
N a ee irn Br 
f Lem ee a h 
— Ze — 7 


bild YoR Ju 8 zuzuschlagen. So erhält Bild 7. Parameter p=1/8 (mod4). Rein imaginäres Argument. 


. man J_/e= +44,7—i9,5 mit dem 


Betrage 48,6. Die Differenz He —J ls liefert i N_1rs = — 337 — 110,5, folglich 
N_jg = — 10,5 +:43,7 mit dem Betrage 44,9. Der Pfeil N_,,s weist as nach links oben 


als der Pfeil N,;s- Schließlich ergibt sich noch Hs = =2 JR Cs = 91,4 +i1,0 mit 
dem Betrage 91,4. Dieser Pfeil weist also fast waagerecht nach rechts. | 


N_, können wir auch nach der Formel 


| BUN er An JOH 00 DENE EN NE SA 
erhalten. Hierin ist — i%25 Jule = —40 +i8 und Hu 05.0,25. = By et 
Summe ist # N_,,g = — 43.7 — i 10,5. Dieser Wert stimmt mit dem sol früher gefundenen 


überein. Siehe die Konstruktion auf der linken Seite von Bild 7. 

Wie ändert sich unser Bild 7, wenn der Parameter p wächst? Gehen 
wir aus von p=0. Dieser Fall wird durch Bild 8 UBEERSEOLE: für das 
y = 0,65 angenommen ist. Dabei wird J,= +1,11 und Hy = —i 0,46. 
N, hat den Betrag 1,20. Wenn jetzt p Er drehen sich die Pfeile Jp 
und H“) links herum, die Pfeile J_,und H\ rechts herum. Bei p=0, 5 H® 
fällt A“ auf J_, und EB und J, erlangen entgegengesetzte Richtun- Bila8. Parameter p=0 
gen. Beip=1 "sind die Funktionen J, N, H\) mit dem Parameter Bi ee 
— 1 entgegengesetzt gleich denen mit dem Parameter +1, was ER 


Bild 9 zeigt. Für dieses ist y—= 1,2 angenommen. Dann wird J, = +10,i 71 und HV’ = —.0,38. 
Der Betrag von N, ist 0,77. Bei p= 2 fallen die Funktionen mit dem Parameter. — 2 mit 


N, 


(2) 
fl o 


J,(i065) 


BEE? °  Emde, Pfeil-Diagramme für Zylinderfunktionen 
(den Funktionen mit dem Parameter +2 zusammen. Für Bild 10, das dies zeigt, ist 18 
angenommen. Dafür wird J,=—0,53 und H= +0,22. Der Betrag vom Nzist 057 
Bei p=3 unterscheiden sich die Funktionen mit dem Parameter —3 N es i Be 
Parameter +3 wieder nur durch das Vorzeichen. Für das zugehörige Bild 11 ist y= RR: 
genommen. Dafür wird J;,= — 0,47 und Hy) = +0,17. Der Betrag von ee Di RER 


m 
" 


Bild 9. Parameter »=1 (mod 4). 


Rein imaginäres Argument -3 
"ty =i1,20, 
H°, (180) 
4 
J,(250) | 
Hr 
\ 
- N; # h 
Bild 10. Parameter p =2 (mod 4). Bild 11. Parameter» =3 Bild 12. Jp(5 i) und J-p (5i) bei veräpderlichem 
Rein imaginäres Argument (mod4). Rein imaginäres Parameter p. Je (di) = 27,2. 
iy=i1,80. Argument iy =i2,50. 2 


Wir zeigen noch die beiden Spiralen, die die Spitzen der Pfeile J, (5) = 5, (9) und 
Hy) (5) = i-v-1 9, (5) beschreiben, wenn sich p ändert (Bild 12 und 13). Die eine Hälfte jeder 
der beiden Spiralen ist das Spiegelbild der andern und zwar bei J, (5) an derreellen Achse, und 
bei H{ an der imaginären Achse. 
Die Spirale für Jstrebt nach dem 
Nullpunkt, die Spirale für HW 
nach dem Horizont ins Unendliche. 
Wenn p eine gerade Zahl wird, R 
schneiden sich die beiden Hälften 
jeder der beiden Spiralen, also 
nach jedem halben Umlauf. 

Die Spirale des Bildes 12 ist 
zwar richtig, gibt aber wegen des 
großen Arguments y=5 eine zu 
einfache Vorstellung von dem Ver- 
halten der Funktion J_,(öy) bei 
festem y. Nach der Formel (III) » 
können J, und J_, nicht zuein- 
ander symmetrisch sein. Des- 
wegen zeigen wir in Bild 14 
noch die Spirale für J,(% 1,6) und 
J_, (1,6). Wenn p eine gerade 
Zahl ist, schneiden sich die bei- 
den Teile der Spirale allen Un- 

 symmetrien zum Trotz. Und 
wenn p eine ungerade Zahl ist, 
sind die beiden entgegengerich- 

se ss) teten Pfeile gleich lang. Für 

Bild13. A, (5°) und H_,(5 i) bei veränderlichem Parameter p. eine ganze Zahl p ist ja sin 22 


#1) (5 i) = i 0,0514. gleich 0. 


ERRERN BE Rz, y enden. ie ae bei ER 
einem andern Maßstab für andere Zah- 
len gelten, die in demselben Verhältnis. 
zueinander stehen wie unsere angenom- SE 
 menen Zahlen. Aber bei Be 
"Argument iy wächst die Besselsche 
Funktion J beständig und nimmt die 
-  Hankelsche Funktion HU beständig 

ab. Daher wächst das Verhältnis 38 
beständig. und nimmt jeden Wert nur. 
einmal an.. Jedes für rein imaginäres 

Argument . ‚gezeichnete Diagramm wird 
„also stets einmal Gültigkeit erlangen, 
aber eben nur einmal. EN ar 

DEAN Ip: merken noch an, daß die bei- Bild 14: 991,6.) und J_ 2 ‚s» bei veränderlichem Parameter. 
den Ausdrücke E en 


RE h Je tiy) HIN (® 2 +iy)= Ä 7) 
| | Ieiy—iNß@iy=HN@—ig). 


einander konjugiert sind. 


Verfolgen wir jetzt die Änderung unseres allgemeinen Diagramm Bild 7 bei Wachen 
dem %, dann sehen wir, daß H{ immer mehr gegen J zurücktritt. Das Diagramm nimmt immer 


7 mehr die Gestalt eines langgestreckten Streckenzuges an. Einige der Strecken werden einander 
Bl immer mehr gleich. Wir lesen so aus unserm Diagramm ab, daß bei wachsendem großem % 
Be die folgenden Funktionen einander immer mehr gleich werden: 
E 1 a Ba 
E Y Nez a HIN, AN=z HN: ER 
= Konjugiert dazu sind bei großem % die Funktionen: 
e : f ; 1 ; : RR RR i 
, J,ty)=-— a = ER:CT Eh) HNiN=,;B in: 
er Irgend zwei von ‘den Zylinderfunktionen sind zwar bei gebrochenem Parameter grund- 
sätzlich voneinander unabhängige partikulare Lösungen der Besselschen Differentialglei- 


chung. Aber bei rein imaginärem und allgemein bei komplexem Argument empfiehlt es sich 


durchaus, für eine Zahlentafel die beiden Funktionen J, und HN zu wählen. Als Gegenbeispiel 
B- betrachte man etwa die Tafel für Y,3(y) und Sys (Y) in unsern Funktionentafeln (S. 285. der 
| 2. Aufl. oder S. 235 der 3. Aufl.). Bei großem y werden beide fast gleich: 1, (8) = 424,39 und 
S115. (8) = 424,40. Von einer linearen Unabhängigkeit kann praktisch keine Rede mehr sein. 
Aus $ kann man übrigens auf folgende Weise einen Näherungswert für 9 berechnen: 


Say 


EDEN 


a a eh et 


| IT 08 
Das im Zähler stehende Produkt der beiden anfangs- en Reihen ist annähernd: 
S2y):8(—2y)= 1 


und weicht gewöhnlich noch merklich von 1 ab. 
| | 24 


5 Bei reellem festem Parameter und beliebiger 


wählen es so: r— 0,9; dann ist J,= 0,99 — i 0,20 und Hi) = 0,34 — 50,23. Daraus fo 
N = — 0,03 +0,65 und AD =1,64—50,17 (Bild 15, 00 Be: 


+0,35 +0,28 und HY’ =— 0,21 — 0,53. Daraus folgt He) 
= 0,91 +i1,10 und N, = — 0,82 +30,56. k 


Be 


10€e1 J e2 # Fe % 
mit veränderlichem Betrag r, aber festem Richtungswinkel o hat keine der Ayllt 
einen festen Richtungswinkel, sondern die sie darstellenden ee rk 
ihre Richtung. Man betrachte etwa die räumlichen Spiralen in de Figuren 1: a | 
Funktionentafeln auf S.260 der 3. Aufl. von 1938. Aber im übrigen g‘ Ye 0 SUR 
Konstruktionsregeln unseres Bildes 7. Wir zeigen das an einem Beispie BORN 


Für dasselbe Argument z = 0,9% (Bild 16) it , = 


Diese Funktionen treten beim Wechselstrom auf, 
"Wechselstrom, der durch einen Runddraht "fließt, verteilt -, 


HR $ 
| #” Jl03i®) 4 re 
Bild 15. Parameterp = 0. Komplexes Bild 16. Parameter p = 1. KRRES 7 
Argument 2 = 0,9 5, Argument u 0,9 05. 


d, ! 5 | 
sich ungleichmäßig über den Drahtquerschnitt (sog. Hautwirkung), er wird nach der Ober- 
fläche hin gedrängt. Die Stromdichte ist nach Größe und Phase proportional zu J,(r #05) | 
und die magnetische Feldstärke zu J,(r i%5). Wenn s den Abstand bedeutet, den ein | 
Punkt im Innern des Drahtes von der Drahtachse hat, und k die komplexe Kreis-Wellendichte, Re: 
soistz— ks. Ferner sei die Stromdichte der reelle Teil von @ (s) Y2 e-i®tund die magnetische 
Feldstärke der reelle Teil von H(s) J2e-ie!. (wo = Kreisfrequenz und t= Zeit). @,ı sei die 
Stromdichte, die bei gleichmäßiger Stromverteilung (also bei Gleichstrom) vorhanden wäre, 
und 25 sei die Drahtdicke. Dann ist 


Gi) _KkbIokk) na Ha _ 5 Jılka) 


G 2 Jı(kb) Ga. 2 Jı(kb)' 
Das Verhältnis des komplexen inneren Scheinwiderstandes zum Widerstand bei Gleichstrom R, 
— auf diese Größe kommt es besonders an — ist 


R—iol; 2J,() Bi 
—- — mit, 
R, 2. Je) 


Dabei ist z = r i%5 für. die Drahtoberfläche s= b zu nehmen. Für r = 0,9 ist m = 1,0084 und 
u = 0,0639” = 7,1°. Bei dem kleinen Wert r = 0,9 ist also die Hautwirkung noch geriug, sie 
beträgt nur 0,84/100 und äußert sich mehr durch eine Phasenverschiebung der Stromdichte an 
der Drahtoberfläche um 7,1°, Die Kurve, die die Spitze des Pfeils m i”“ beschreibt, wenn sich r 

ändeft, zeigt Bild17. Das Wirkwiderstandsverhältnis ist 
R/|R,= mcosu und das Blindwiderstandsverhältnis ist 
o Li/R,=msin u. Wenn ferner Z;. den Grenzwert bedeutet, 
dem die innere Induktivität des Drahtes bei verschwindender 
Frequenz zustrebt, so ist 


R, 8i 8 
und 
ei Re 8 ; 
Lio a tn 5 m (sin u + %c0s u) 
also das Induktivitätsverhältnis 
er RA 2. L; RR, 


m und „ als Funktionen von r. Li La 


etwa ine! Frequenz, von 10 


{5 K 


Bl ee udeR at eine ptäpe! zurü ück, di 
lie auf die Frage nach den Zylinderfunktionen 1 für Be reell 
Hier treffen wir auf einen Umstand, den wir bisher mit Stillschweigen über- 
Nullpunkt Bent sich die ee im allgemeinen un sind also 


” 


radeswegs. auf die negative übergehen, so RERSERR: wir ins Unbewisse: geraten. Er 
E Wer müssen den Nullpunkt auf der Riemann schen Fläche umlaufen, die sich um den Null- 
punkt. ‚herum windet, um zu bestimmten ‚Werten zu gelangen. Im Nullpunkt werden überdies 
ie Besse schen. Funktionen SR, mit einem negativen. ‚Parameter und alle übrigen Zylinder- 
2 eg unendlich groß. ‚Nur In ist Saweruß, wenn ‚Jen Parameter, n aus positive ganze Y 
Zahl ist, i 
3 "Wie sich die. Funkuionswekte, für las Argument ver je El (m ganz) aus Bi Funktions- 
St werten berechnen: lassen, die für z= rie gelten, das lehren die Umlaufsformeln der ee a 
funktionen. (Siehe etwa Weyric "hs Zylinderfunktionen, Leipzig 1937, S. 34.) Wir wollen 
 uns.hier auf reelles Argument « i2% (k ganz) beschränken. Wenn x in x it? übergeht, ‚wenn also e 
das Argument e einen ‚halben Umlauf links herum um den Nullpunkt macht, so wird Ye 


18 SLR RR Parc N (vi = HD (a). E N i 


"Greilen ‚wir ut das Beispiel des Bildes 1 zurück: P=1/3 u 21, 68. Ele ist, 


HN) 1 (1,68) = 0, 614-9212 (Bild 18). Diesen Pfeil haben wir um 4/3” = 1,333” oder um 120° vorzu- 
drehen. Dadurch erhalten wir 0,61 i1-121 mit dem Richtungswinkel — 101°. ‚Das ist also einer 


2 Es ; der Werte, die Hi) für das Argument — 1,68 annimmt. 


7 


3 Wird dagegen .der Nullpunkt rechts herum um- Pal 

E: laufen, so nimmt H\', folgenden Wert an: 
En | HY), (1,68 52) — i+2!3 HÜ), (1,68) + :2/3 7,5, (1,68). 

= Wir drehen also H() um 60° vor und den waagerecht #448) 

mach rechts weisenden Pfeil :2J um 60° zurück. Die Haar) 


Addition der beiden so erhaltenen Pfeile ergibt einen 
Pfeil von der Länge 0,80 und dem Richtungswinkel — 30° 
' (Bild 18) oder die komplexe Zahl 0,69 — i 0,40. Das ist Hrtı88r® 
also ein zweiter Wert fürH für das Argument —1,68. #ısi?) 
Nach einem vollen Umlauf gelangt man zum Aus- H 
gangspunkt zurück, auf. der Riemannschen Fläche 


j aber zu einem anderen Punkt. Man erhält so zwei weitere € 
E Werte für ZU mit dem Argument +1,68, wenn der ä RER 
Nullpunkt erst links herum, dann rechts herum umlaufen a i78Jz IR 
u Nach einem vollen Umlauf links herum erhält Fan Est: ER sit?) für 21,68. Roollen f | 
 H37 (1,68 644) = 741° 91, (1,68) — (1 + 009) 2 I, (1,68) positives Argument 2 + 1,68 = 1,6840 1 
: - — j-213 7 -F ra. = 1,68 ee Reelles negatives Argument 
- : +2 
“ Der Pfeil H% ist rechts herum, der Pfeil 2 J links herum — 1,68 - 1,68 = 1,6847. 
um 120° zu drehen. Das ergibt . 


0,61,571:121 41,16 s+1:33 — — 0,69 + 0,40; 

das ist ein Pfeil von der Länge 0,80 und dem Richtungswinkel + 150° (Bild 18). 

Zu einem ganz anderen Wert gelangt man nach einem vollen Umlauf von & rechts herum: 

H\h (1,68 i-*) = i+28 A), (1,68) + 521° 2 9,7, (1,68). 
Jetzt ist H% um 120° links herum und 2 J um 60° rechts herum zu drehen. Das re! 
0,61 v» 545 4 1,16 59:667 — +0,12 — 0,60. 
Das ist ein Pfeil von der Länge 0,61, der unter 79° nach unten geneigt ist. (Bild 18.) 
24* 


7 Ic 1 e 
Hls(zit) = — Hi), (e®). 


in der Tat lehrt die allgemeine Umlaufsformel für m halbe Umläufe? nr 


Heide — 


sin mp2” 


RER, 
sin pP % 2 Js (2) D 


daß die rechte Seite immer dann = 4 HW) ist, wenn das Produkt mp aus dem Parameter p 


und der Zahl m der halben Umläufe, aber nicht der Parameter p selbst eine ganze Zahl ist. Hier, 


wo p= 1/3 ist, findet das zum ersten Malnach m = 3 Umläufen statt. Dadurch gelangt das 
Argument aus einem Blatt der Riemannschen Fläche in ein anderes. Ein ähnliches Ver- 
halten zeigt schon die einfache Funktion w = z!/®, nur daß dabei die Beträge der einzelnen 


Funktionszweige nicht verschieden sind. 


Wenn man also beim Übergang von der positiven reellen Halbachse zur negativen drei 


halbe Umläufe um den Nullpunkt einschaltet, so verhält sich Z us als ungerade Funktion. 
Das gilt nicht nur für die reelle Achse, sondern auch für jede andere (im allgemeinen schiefe) Achse. 


Eingegangen 13. 5.48. 


KLEINE MITTEILUNGEN 


Eine einfache Methode zur Ermittlung der 

charakteristischen Daten eines gedämpft 

schwingenden Systems zweiter Ordnung mit 

an einer neuen Auftragung der Resonanz- 
urven. 


Es ist im allgemeinen bei Resonanzversuchen üblich, 
die Amplitude und den Phasenverschiebungswinkel 
graphisch in Abhängigkeit von der erregenden Frequenz 
aufzutragen, um daraus die interessierenden Daten des 
Schwingungssystems zu entnehmen. Dies geschieht 
mit Hilfe bekannter, aber meist nicht einfacher Me- 
thoden. Zu. weit einfacher auszuwertenden Resonanz- 
kurven gelangt man indessen, wenn man, anstatt von 
den unmittelbar auffallenden Merkmalen der Ampli- 
tude, Phasenwinkel und Frequenz, von Größen aus- 


geht, die sich aus diesen ableiten lassen. Diese ver-' 


wendet man zweckmäßigerweise sowohl zur theoreti- 
schen Darstellung als auch zur praktischen Auswertu 
von Resonanzversuchen. Dies aufzuzeigen ist das Zie 
der nachfolgenden Darstellung. 

Im Interesse einer möglichst allgemeinen Anwend- 
barkeit der Ergebnisse möge von dem Energieinhalt 
ausgegangen werden. Es liege ein dynamisches System 
von einem Freiheitsgrad (x) vor, dessen Energieinhalt 
durch 


1 t 1 t z 
Zt Spritze) + 
0 0) 


geßeben sei. Durch Differenziation nach der Zeit t 
erhält man daraus 


-" + 
ie “N 4 6 
Führt man die Kreisfrequenz &, = \ P_ der unge- 
” * . x 
dämpften Schwingungen ein, so wird daraus 


22 eine oin— lt): 


°p 
Y%p °k K 


Die Eigenbewegung des Systems wird bekanntlich 
durch die Gleichung 


ee 


2yep Fu E ( cp ) un 
%—%ep'e * co8[ @&l/ L-| F——) Do 
las 
(3) 
beschrieben, wobei x;, und 9, die Integrationskon- 


stanten sind. Nach Lehr [1] führt man die von ihm 
als „Dämpfungsmaß“ oder „Dämpfung“ bezeichnete 
dimensionslose Größe 


ED DR 


ein. Es wird dafür der Ausdruck Dämpfungsver- 
hältnis vorgeschlagen, da D physikalisch identisch 
ist mit dem Verhältnis der tatsächlich vorhandenen, 
Dämpfungskonstanten c,, zur sog. kritischen Dämp- 


% Pa > 
fungskonstanten ey y-j, = 2 Cx * Cp, welchen Werte,,, 


wie aus (3) hervorgeht, im Falle des aperiodischen 
Grenzzustandes annehmen muß. Es ist also 


DD °D 


Cpkrit 


Mit dem durch (4a) definierten Dämpfungsverhältnis 
vereinfacht sich die Differentialgleichung (2) zu 


. 1 
i+2Do,:+o!z= fit)» 2...» (2a) 
K 


und deren Integral im Falle f(t)=0 zu 


” —Do 
Ko Rore 


"008 (m Dit 9) (Ba), 


woraus man die Kreisfrequenz &, der gedämpften 
Schwingung zu 


erhält. 


Br im allgemeinen Nbakenähiei, man nun ne Voir EN 
 rungsverhältnis is v =: 40 und die Phasenverschie- 
Br bung a in Abhängigkeit vom. ‚„Abstimmverhältnis‘ © ax 

Zweckmä äßiger ist es, das Quadrat des Reziprokwertes. , 


an)» u ‚Funktion von en JE und cotar als Funktion 


we (9) © 
5 von ie au auftragen, wie esin Bild 1 "und 2 Baer 


Wie’ man nach der Elimination von D aus (11a, p) 
Aka genügen &, und 77, der ODER UNE 


m—1l-—#. : (12). 


Gl. (9a) stellt eine Schar von Parabeln — Resonanz- 
parabeln — dar, welche alle unter sich kongruent, nur 
gegeneinander parallel verschoben sind. Der Scheitel 
gleitet dabei an der Parabel (12) entlang, welche ihrer- 
seits selbst ebenfalls kongruent zur Parabelschar (9a) 

ist. Der Schnittpunkt von (12) mit der &-Achse hat 

die Koordinate &,= +1. Der Resonanzpunkt ist ge 
kennzeichnet durch die Koordinaten &,,n,[Gl. (11a, b)], 

d.h. die Kreisfrequenz ®, im Resoanpunk Sn 

sich aus 


I 2 
| (2) -1-—2D% 
wo. | 


und die Amplitude aus 


A(0)\? 
(Oman, 


05 


Bent 
(ee fe)" 


Bild 1. 


Pe) 


ist. Zur analytischen Darstellung dieser Kurven geht Da die Eigenschwingung, wie man aus (5) entnimmt, 
man von (7a, b) und (8) aus und erhält durch die Beziehung 


A(0)\? a2, 12 @oe\? 3 - 
Ko) =) + er Be 


in Abhängigkeit von 


a | | u 
Für die Phasenabhängigkeit gelangt man eben- 

Be Kurven zweiten ae wenn man statt 8 
ort, | ea 


vol... (13a) 


R 9) i : { "ih 
Ur PR N N Mn 13b 
STR \ ) 
 aufträgt. Aus (7b) erhält man mit (13a,b) 


Ua BT ae 
Hole ). . 28 .. (14), 
uw=—D-'v+y1+DWw . (14a) 
bzw. 


EIDWETENL 25 Mahl. 


Die Diskriminante A dieser Kurve zweiten Grades ist 
= —D?2. Es ist also stets A<0 [2]. Es he 

sich daher für jeden Wert von D eine Hyperbel — Pha- 

senhyperbel. Ihre Asymptoten sind durch 


und 


gegeben. Der Mittelpunkt aller Hyperbeln der Schar 


(14) fällt‘ daher mit dem Koordinatenursprung zu- 


sammen. Weiterhin gehen alle Hyperbeln durch den 
Punkt w= +1, v%y=0. Bemerkenswert ist, daß die 
zu diesem Punkt gehörige Ordinate der Asymptote 


1 
(15b) den Wert u = + SD besitzt. Den Reziprok- 


wert dieser, Größe, d.h. 2 D unmittelbar erhält man, 
von der bekannten Tatsache [3] ausgehend, daß 


u =2D 


= — 1 T%b= +1 
ist. Man kann dieses Ergebnis verallgemeinern, indem 
man die Differenz derjenigen Kreisfrequenzen », und ®, 
bildet, für welche die dazugehörigen Phasenverschie- 
bungen a, und &, symmetrisch zu /2 liegen, d.h. für 
welche cot &, =—cot a,, oder nach (13a) y =—ı, = ö 
ist. Dann gilt, wie man aus (14a) entnimmt, 


2. (16). 


Aus (14) entnimmt man, daß alle Hyperbeln aus einer 
einzigen (z.B. derjenigen für D=1) durch affine 
Verzerrungen hervorgehen, indem man nämlich alle 
v-Koordinaten proportional D verkleinert und die Ab- 
szissen u unverändert läßt, 

Die Darstellung in dimensionslosen Zahlen gemäß 
Bild 1 und 2 beansprucht nur theoretisches Interesse. 
Praktisch liegt die Aufgabe so, daß D und &, meistens 
nicht bekannt sind und daß sie erst aus Resonanzver- 
suchen ermittelt werden müssen. Man kann daher 


er =2D% 


0) 
nicht &E = (2) als Variable wählen, sondern man 
0 


muß für die Frequenz irgendeine passende, aber sonst 
willkürlich wählbare Maßstabeinheit &y einführen. 
Mit den neuen Variablen 


Ya o 2 
pn (2) et rt (17a) 
und 
'- “(171 
WRITER 2) sr ‚) 
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BB, Are 
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Unter Benutzung dieser Beziehungen 
die Aus von Resonanzversuchen 


zu entnehmen ist, ’ } 
gilt für die Phasenhyperbel (14) (Bild 4). Von den 
Asymptoten ist allerdings nur die eine (v =0) bekannt, 
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Bild 4. 


Die zweite ermittelt man mit Hilfe der bekannten 
Tatsache [5], daß die zwischen den Asymptoten und 
dem einen Hyperbelast auf einer Geraden liegenden 
Strecken einander gleich sind. Man legt daher in dem 
Phasendiagramm (Bild 4) eine beliebige Gerade B’C’, 
welche die aus dem Versuch sich ergebende Hyperbel 
in, den Punkten B’ und C’ ‚und die eine Asymptote 
(w = 0) im Punkte A’ schneidet und trägt über C” 
hinaus eine Strecke von der Länge A’B’ ab. Auf 


4 


es Tl 2 YE 


NA Die Vorteile dieser Methode sind i in n folgenden Punkten 
Be rd? erblicken: 

1. Die Resonanzkurve ist, eine zum Resonanzpunkt 
„Bymmebrische, analytisch einfach darstellbare Kurve. 
2. Es ist nicht notwendig, Amplitude und Phasen- 


gt 


it erhält die Bigenfroquenz ©. aus @? 


j ‚Wie. man sich durch eine einfache eelinube AZ 
zeugt, führt das Integral der Differentielgleichung 


Ban Bade nd Ad u m > m an 
D ach 


und bi 
"Ordinaten. an ist, ale: wie ‚man. aus 
VRR N { 


Baer he Be 
EREET we “ x A . (9) ; 
‚mit. denen des ersten Problems und kam dann darauf, Br 


pe Resonanzpunkt oder in dessen un- 
'mittelbarer Nähe zu messen; diese ändern sich dort 
besonders bei schwacher Dämpfung mit & sehz stark, 
was meßtechnisch sich oft unangenehm auswirkt. 
3. Man entnimmt @; unmittelbar aus dem m Dagpacn 
Ab os a% o 


&+2 Dove + 052 = bw? cos nt RR (20) 


- zu einer ‚stationären a von.der Form 


2 = Bio) ost — Moll... > (20a), 
wenn 
Bio) B(o) 1 
Re IR LEE Ta) 
b B(o , "2 
yYız@jf+b»® 
und 
20 
18 Bo)= — er (2lh) 
-@) 


ist. Durch Vergleich stellt man fest, daß (7b) und (21b\ - 


identisch sind und daß (21a) aus (7 a) hervorgeht, wenn 


2 
man € = (2) durch dessen Reziprokwert ersetzt; 


d. h. den durch die Gl. (20) beschriebenen Fall kann 
man mit Hilfe der Phasenhyperbel ohne Anderung und 
mit Hilfe der Resonanzparabel (9a) nach der Substi- 


-tution € = — auswerten. Bekanntlich liegt dieser Fall 


bei einer durch periodisch wirkende Massenkräfte her- 
vorgerufenen Störung [3a] und bei Schwingungsmeß- 
geräten mit relativer Federung und Dämpfung vor [4]. 


“Sehrifttum. 

Ba E. Lehr: sehr T. Bd.. Berlin 1930. 
S. 240 

[2] Siche; z. B.: Hütte, I. Bd, 27. Aufl., S. 135—136. 

[3] E. Lehr: Schwingungstechnik, II. Bd. Berlin 1934. 
8. 149. 

[3a] E. Lehr: Schwingungstechnik, II. Bd. Berlin 
1934. 8. 176—182. 

[4] F. Rixman: Wirkung der Dämpfung bei Schwin- 
gungsentstörung und Schwingungsmessung. 
VDI 85 312—315. 

[5] Siehe z. wien T. Bd., 27. Aufl., S. 138. 


München. H.St. Stefaniak. 


ch langeı 
wird die Te im i 
telle des Wärmeeintritts. 'eiter 
e.die Wärmezufuhr gi muß, da- 
immten Zeit; unkt 


die ee en ist. Me 
Für das erste Problem ist die Lö 

bekannt; ich bringe sie hier nur 

zweite Problem. Bei dessen Untersuchung fand 

erst eine gewisse Übereinstimmung der Gleichungen 


ua 


daß beide Probleme identisch sind. Die, Lösung des 
zweiten Problems ist schwer zu übersehen, während 
das erste (wenigstens für mich) viel anschaulicher ist. 
Aus diesem Grunde will ich hier Beide a, 80 Mi 
meinsam behandeln. ehe 


I.Problem. ‚Der. Feasaslelende, Stab. re 
sich von x = 0 bis x = ©. Zur Zeit i wird durch die 
' ‚Stelle x in der Zeiteinheit die Wärmemenge q.(t) durch- 
geleitet. Die im Anfangsquerschnitte je Zeiteinheit 
zugeführte Wärme ist q9 (?), die Zeit ist hierfür mit z 
bezeichnet. Die ganze, bis zum Be i ee Ib; 
Wärme wird; u 


n=fni Bay 
oder Br DER ; s ) IR 
| ala 
‘ Die Temperatur im Stabe ist 9(z,t). Setzen wir 
alle physikalischen Größen gleich eins, so erhalten wir 
die Gleichungen: 


0.0 | 
Er =— 357, (Wärmebilanz) 
0» | 


—3%:  (Wärmeleitung) 
Hieraus folgt die allgemeine Differentialgleichung: 
99. 09 
o2: 9 


Zuerst behandeln wir den Fall, daß zur Zeit nl). 
plötzlich die Wärmemenge Q eingeleitet wird, 
die sich dann allmählich im Stabe verteilt. Die Lösung 
erfolgt mit Hilfe einer Ähnlichkeitsbe- 
trachtung, die auch für den zwei- und drei- 
dimensionalen Fall durchgeführt werden kann. Bild la, 
b und e zeigt: den Temperaturverlauf im Stabe zur 
Zeitt=0, t und t+.dt. Ändert man die Maßstäbe 
der Bilder Ib und Ic in geeigneter Weise, so müssen 
sie ineinander übergehen. Die Stelle ® im Bilde 1b 
entspricht der Stelle x (l+.cdt) im Bilde le, während 
die Temperaturen an diesen Stellen und®# : (1 -+ cdt) 
=d9(1— cedt) sind. 

.Im Bilde 1b betrachten wir die Dosnurkn Stellen . AR 
x und 2 + de, wobei wir da = xcdt machen. Der Sc mpeg 
‘Wert c ist für den ganzen Stab einheitlich, wird aber 5 
eine Funktion ‚von { sein. 

Vergleichen wir die Bilder 1b und lc, so sehen wir, 
daß wir an der Stelle x + d« in der Zeit di eine Tem- 
peraturänderung erhalten von der. Größe 


00 
ia, c.din 


Dieses ist aber die Temperaturänderung FE Hieraus 


folgt 
0%) 09 


0% 
9 1l—ced)— +25, zei) = 


Durch Ita ih Bieruun: Fa E nr 5 * 


l 


ı X(Ircdr) ı ° u; 
RNENEE ra 
Bild 1. 


Die von t abhängige Integrationskonstante wird gleich 
Null, da für x = 0 auchz, = 0 ist. Aus Gl. (4) folgt 


d) 


Hierin ist k eine Funktion von t, Setzen wir jetzt ® 
nach Gl. (5) in Gl. (1) ein, so erhalten wir ce = 1/21 


und k = kıfli. Berücksichtigt man, daß die gesamte 
Wärme im Stabe ! dd. gleich der zugeführten Wärme 


daher e"i2, 


Q sein muß, so läßt sich hiermit k, durch Q ausdrücken 
und wir erhalten: 


ee ENDEN SET REEL 


Wird die Wärme allmählich zugeführt, so daß 
go = 9% (7) ist, so erhalten wir für die kurze Zeit dt die 
Wärmezufuhr 9 dr. Hieraus folgt die Gleichung 


t—T 


Sy) = Hr AT dar m. 


Zur Zeit A 20, id de 


ER EN nach Le 


Vom Zitpunkte 1, an il ine weitere Wärmenahr 
derart erlolgen, daB. Ai Tentperatur 95 = 


wer 


halten bleibt. Die Bee rn 


77 DEE [we Bi 


Bild 2 zeigt, wie ich mir den Vorgang 


Wir wollen die Zeiten vom Zeitpunkiag; an rechnen 
und bezeichnen hierzu 


hg 1-u=t. 


ıi-E 


Dann folgt aus Gl. (9) 
t* 
EIN. 
Yyer+t + Ye—# 


x 10 (ı*) dt — 


Se 


denke. 


1 


Yi. 


et. DR 


(T 
. Mel ran | te 
 ) GR 


ö Führen wir im RER die Veränderliche 2ER ein | 8 
_ und machen wir den Ansatz q” ig >=.on DAR: so 
erhalten x wir A 


d er = > en. 


@ ya 
ja — (7) 


De ER EL eine reine Ar, ee folgt = SE. 
3 &=n; jetzt kann cn für jeden Wert n berechnet 
fe ER - werden; hiermit ist qn (7*) bestimmt und g, (r*) als 


= 5 Reihenentwicklung gefunden. EN FRRENEE Bildd. SER 
: U. Problem: Bei einem schwingenden in ET DEREN UNE DEREN I TERREE WEN 
massensystem sei das Potential V der rückführenden _ RENTE ERS Y 
Kraft eine beliebige Funktion des Ausschlages x. Wir RAN RR: or » 
j wollen nur annehmen, daß V mit x zunimmt und EIN GT Pi zw) - a 
4 I (+2)=V (—) ist. Ist der BD ENIBE a en de dV =x (v) 4. 
En ‚gegeben, so ed die Geschwindigkeit und 1.10) gibt: av 
57, | ER x SEE 
RE BR je 3m [/le) — V(@)] 1 
Be: T (a) = “Em | re re 
u und die a - ’ VEIT RNE 
Er- = 3 5 L 4 
2 “ Diese Gleichung stimmt vollkommen mit Gl.(8 
= Fe? T —4/2m Va V.v(a) — Ve) - - - (10)... überein. Es ist nur an Stelle von x’ (V) die Wärme)- 
4 = zufuhr g, (r) zu setzen. Als Beispiel nehmen wir 


Wir verlangen, dab für a = & = 8; die Schwingunge- anwächst, daß die Schwingungszeit praktisch unend- 


I 
2 Er an, daß von a—=0 bisa =, das Potential so langsam 
\ zeit T (a) konstant ist. Gibt es hierfür nur die triviale 


ij ! lich groß wird, Im Punkte a, habe das Potential einen. 

Br Lösung V(e)=—-c2?? \ plötzlichen Anstieg, sodaß man für «a = (a, + &).—o 

eine endliche Schwingungszeit erhält. Von dieser an 

EN en nach kann man für den ersten Soll Y(a) so mit « wachsen, daß die Schwingungszeit 

Ausschlagsintervall bis a, eine beliebige Funktion für konstant bleibt. Wir tragen © als Funktion von V 

L das Potential annehmen. Man muß nur für das zweite auf, Bild3. Von «— 0 bis a, wächst V nur gering- 

f Intervall a, <a <a, das Potential entsprechend er- fügig. Für x =a, vergrößert sich V plötzlich auf 
mitteln. Es fällt eine gewisse Ähnlichkeit zwischen einen Wert V (a,). Der weitere Verlauf ist, noch unbe- 
der G1.(10) des II. Problems und Gl. (8) des stimmt. Bild 3 zeigt den Verlauf von 2’ —dx/dYV 
I. Problems auf. Um eine Übereinstimmung als Funktion von V. Das ist aber derselbe Verlauf wie 
zu erzielen, betrachten wir V nicht als Funktion beim Beispiele des 1. Problems, ‚siehe Bild 2. 

von x, sondern x als Funktion von V; das bedeutet, Ich überlasse es den Lesern, sich zurechtzufinden. 
daß nicht jedem Ausschlage ein Potential, sondern Manches ist ohne Worte schon klar. 

5 aeden Potentiale ein Ausschlag zugeordnet ist. Schlewecke über Bockenem. G.-Weber, 


Ei Me \ er, 


Er _ BUCHBESPRECHUNGEN RE, { 


- Dr.W.W eizel (0. Prof. a.d. Universität Bonn), Ein- zum Regpetieren von Vorlesungen, Übungen und an- 
führung in die Physik, 1. Band: Mechanik, derem Physikunterricht dienen. In dem ersten der drei 
2. Band: Elektrizität und Magnetismus (2. verb. Aufl.). Bände des Gesamtwerkes werden die Grundlagen der 
Meyers kleine Handbücher Bd. 42. 145 S. mit 119 Abb. Mechanik nicht als Extrakt zum Auswendiglernen,son- > 
bzw. Bd.43. 159 S. mit 115 Abb. Leipzig 1947. Ver- dern mit der Anleitung des erfahrenen Lehrers zum tie- 
lag: Bibl. Inst. Preis je 2,50 DM. feren Verständnis geboten, wobei 119 einfache Skizzen 
Diese handlichen Bändchen, fast unveränderte Neu- und Photobilder sowie wenige elementare mathema- 
drucke des ersten Auflage, sind aus Vorlesungen für tische Ableitungen helfen, auf historische Hinweise aber 
Mediziner hervorgegangen. Wie im Vorworte betont weitgehend verzichtet wird. 
wird, sollen sie im wesentlichen nur zur Ergänzung und Abgesehen von kleinen Verbesserungsvorschlägen, 


arWuT /- 
Er 
VE 


Zar Ar Fr 
ER TEN 


a nt 
wW " 


= 


Dh Aue a 


Ba 


Zweckbestimmung Z 
nitionen des Differentialquotienten und des bestimmten 
Integrals zu verzichten und scharf zu unterscheiden 
zwischen den Zeichen da, Ax und dx, anstatt nur de 


anzuwenden, ] Zahlenbeispiele für die Benutzung 
der Formeln ana Maßeinheien würden diese Einfüh- ]i 
‚ung erheblich wertvoller machen... . 


. Das fesselnd geschriebene Buch kann den Teil „‚Me- 
chanik“ der schlechten Repetitionsbücher vorteilhaft 
ersetzen. Es ist ihm a em eine seinen begrenzten 
Zielen entsprechende Verbreitung, z. B. bei Physikvor- 
lesungen an Volkshochschulen, zu wünschen. _ 

Der zweite Band beginnt mit dem Coulombschen Ge- 
setz und schließt mit der Rundfunktechnik. Ein Fünftel 
des Raumes ist dem Elektrizitätstransport durch Gase 
mit zahlreichen modernen Anwendungen met. 
Wie die meisten deutschen Lehrbücher der Elektrizi- 
tätslehre leidet dieser kleine Leitfaden trotz der klären- 
den Arbeit des Ausschusses für Einheiten und Formel- 
größen (AEF) unter dem Wirrwarr der noch fehlen- 
den Einigung über die Schreibweise der Gleichungen 
und die Definitionen der Größen. Wenn z. B, das Cou- 
lombsche Gesetz mit dem Faktor 4e, im. Nenner als 
Größengleichung geschrieben wird, so gilt diese für be- 
liebige Einheiten, also nicht nur, ‚‚wenn die Ladungen 
in Cb, die Entfernung in Meter und die Kraft in Groß- 
dyn gemessen wird“. Trotz meiner Befangenheit in 
dieser Richtung muß ich doch bemängeln, daß der als 
„ungeheuer wichtig‘ bezeichnete Ferromagnetismus 
zusammen mit Para- und Diamagnetismus bei nicht 
einwandfreier Darstellung noch nicht 2% des Buches 
beansprucht. Grundbegriffe wie Hysterese und Rema- 
nenz kommen überhaupt nicht vor. Eine Hysterese- 
schleife ist unter den 153 Abbildungen nicht zu finden. 

Weitere, kleinere Mängel schließen nicht aus, daß 
auch dieser zweite Band, besonders dank einiger vor- 
züglicher Abschnitte mit lehrreichen Abbildungen, in 
seinen vom Verfasser ausdrücklich gesetzten Grenzen 
nützliche Dienste leisten wird, M. Kersten. 


Dr. Friedrich Sehilling (em. o. Prof. d. Techn. Hoch- 
schule Danzig), Die Bewegungslehre im 
nichteuklidischen hyperbolischen 
Raum. Bd.I: VI+1388.,, Bd. I: I’+149S8, 
m. 118 Abb. München 1948. Leibniz-Verlag. Preis 
brosch. 35.— DM. | 

Hat man die Bewegungen eines Raumes, so hat man 
den Raum selbst. Seit Klein beschreibt man geradezu 
methodisch Räume durch die darin möglichen Bewe- 
gungsgruppen. So gibt Verf. nicht nur irgend etwas 
Spezielles im hyperbolischen Raum, sondern eine 
Lehre dieses Raumes selber. 

Das Bauwerk der nichteuklidischen Geometrie hat 
4 Eingänge: 

1. Die Entwicklung der Sätze aus den Axiomen, 
Dies ist historisch am treuesten und philosophisch am 
aufschlußreichsten. 

2. Die Geometrie auf Flächen konstanter Krüm- 
mung im euklidischen. Raum. Das ist im 2-dimensio- 
nalen sehr anschaulich, im 3-dimensionalen nicht mehr. 

3. Als Sonderproblem der allgemeinen projektiven 

seometrie. Das ist die eleganteste Methode und gibt 
eine Überschau von höchster Warte. 

4. Verf. bettet den nicht-euklidischen Raum in einen 
euklidischen mit rechtwinkligem Koordinatensystem 
ein, Dadurch erreicht er, daß er so gut wie keine Vor- 
kenntnisse vorauszusetzen braucht, immer auf dem 
Boden des voll überschaubaren, leicht faßlichen, nach- 
rechenbaren bleibt, und für alles, was er behandelt, 
am Zeichenbrett auszuführende Konstruktionen an- 
geben kann. Die übersichtliche Gliederung und An- 
ordnung tut ein.übriges, das Verständnis und das 
Durcharbeiten zu erleichtern. Er beginnt mit Dre- 
hungen und Schraubungen um die x-Achse, um erst 
von da aus zu allgemeinen Bewegungen aufzusteigen. 


umes 
. plexen Ebene, 


- speziell über die Bewegungslehre zu vermitteln. 


die vornehmlich den Verfasser interessieren, möchte ich ' Diese kann 
‚hier empfehlen, auch in einem Buche dieser Art und _ 
nicht auf einwandfreie kurze Defi- 


ge h Achsen der 
je zweier von Ahneh. — Die B indung der Trigono- 


erhält, wenn man auf krummen Flächen, insbesondere 
Kugeln und Abstandsflächen von er die kürze- E 
an ng; „Geraden“ nennt. Sie sind euklidisch "jee4 
elliptisch. 2 IN 


oder 
as Buch ist in gleicher Weise geeignet, dn An- 
fänger in die nichteuklidische hyperbolische Geo- 
metrie einzuführen, wie dem Kenner neue 

Dresden, Ott-HeinrichKeller. 

Archiv für Mathematik, Bd.1, Heft1. Herausge- 
geben vom Math. Forschungsinstitut Oberwolfach, 
Schwarzwald unter Leitung von W. Süss, Freiburg 
(Breisgau). 888. m. 6Abb. Karlsruhe. Verlag 
G. Braun GmbH. Preis 8 DM. ö 

Das Archiv für Mathematik, von dem das erste 
Heft vorliegt, soll zweimonatlich erscheinen. Es soll 
kürzere Originalarbeiten und Selbstreferate, Kleine 
Mitteilungen, Zusammenfassende Berichte, Nach- 
richten aus dem mathematischen Leben und Buch- 
besprechungen bringen, ist also ähnlich gegliedert wie ‘7 
die Z. angew. Math. Mech. Sein Hauptzweck ist die 
rasche Veröffentlichung kurzer Berichte über neue 
wissenschaftliche Ergebnisse, aber so, daß der Kenner 
die wesentlichen Schritte der Beweisführung versteht. 
Von den im ersten Heft erschienenen Arbeiten inter- 
essieren die Leser der ZaMM insbesondere die fol- 
genden: Blaschke, Über Integrale der Kinematik, 
Grammel, Eine Verallgemeine der Kreis- und 
Hyperbelfunktionen (Selbstreferat eines Aufsatzesim 
Ingenieurarchiv 16 (1948)) und Strubecker, Zur 
graphischen Integration der linearen Differentialglei- 
ehung n-ter Ordnung (u, D+ 1): (-,D+1)-... 
(,D+Dy=s(k). 

i dem vielfach sehr unregelmäßigen und seltenen 

Erscheinen mancher Fachzeitschriften ist eine Zeit- 
schrift, die sich das obengenannte Ziel gesetzt hat, ins- 
besondere, wenn erreicht wird, daß die Arbeiten trotz 
Kürze verständlich sind — es ist dasin entsprechenden 
Zeitschriften des Auslandes oft nicht der Fall — sehr 
zu begrüßen, 


Dresden. 


Willers. 


Dr. F. Hund (o. Prof. a.d, Univ. Jena), Einführung 
in die Theoretische Physik, 2.Band: 
Theorie der Elektrizität und des Magnetismus, 3. Band: 
Optik. Leipzig 1947. Meyers kleine Handbücher Bd. 
47/48. 297 8. mit 75 Abb. bzw. Bd.49/50. 2088. mit 
100 Abb. Verlag: Bibl. Inst. Preis brosch. je 5.80 DM. 

Dem unterzeichneten Referenten kommt es nicht zu, 
ein Werk dieses Verfassers zu ‚‚kritisieren‘‘, Wohl aber 
darf er vom Standpunkt des Nichttheoretikers darauf 
hinweisen, daß dieses kleine Buch nicht nur einen erst- 
klassigen Begleiter für Vorlesungen darstellt, sondern 
auch einen ungewöhnlich verbauen Ratgeber für 
Mathematiker, Experimentalphysiker und andere ‚‚ne- “ 
benamtliche Theoretiker‘, die sich um einwandfreie 
Vorstellungen und Definitionen aus der klassischen 
Theorie des elektromagnetischen Feldes und der Optik 
bemühen und bereit sind, sich in die konzentrierte 
Sprache dieser inhaltreichen Einführung zu vertiefen. 


Standpunkte klar zur Geltung kommt. Auf die ele- En 
'gante Behandlung der Lorentzinvarianz der Maxwell-. 
‚ schen Gleichungen, den historischen Anhang und die 


EINGEGANGENE BÜCHER 


L 


dritten Bandes tragen di 


a 


Bei der’ Schriftleitung sind folgende Bücher eingegangen (ausführliche Besprechung bleibt vorbehalten). 


Dr. Gerhard Kropp, Beiträge zur Philoso- Dr.Ing.habil. Alf Pflüger (Dozent ander Technischen. 


phie, Pädagogik und Geschichte der 
Mathematik. Mit,einem Anhang: Die mathema- 
tikgeschichtliche - Forschung. Geometrische Integra- 
tionsmethoden bei Lalouve£re.. 103 S. mit 7 Abb. Berlin 
1948. Fr. K. Koetschau-Verlag. { 


Prof. Dr. Wilhelm Blaschke, Analytische Geo- 
metrie. (Bücher der Mathematik und Naturwissen- 


. „.sehaften. Herausgegeben von Dr. Henry Poltz). 152 8. 


mit 67 Bildern. Wolfenbüttel-Hannover 1948. Wol- 


‚fenbütteler Verlagsanstalt G.m.b.H. (Notdruck). 
' Preis brosch. 10,50 DM. 


Dr. Heinrich Blasius (Oberstud.-Rat an der Ing.- 
Schule Hamburg), Mechanik, Physikalische 
Grundlagen vom technischen Standpunkt. I, Teil: 
Statik. -Vierte Auflage. VIII + 247 8. mit 260 
Fig. und 113 Aufg. Hamburg 1948. Verlag Eckardt 


- & Messtorff. Preis: brosch. 9,— DM. 


- Dr. Ing. habil. Hermann Athen, Ebene und 
sphärische Trigonometrie (Bücher der 
Mathematik und Naturwissenschaften. Herausgeber: 
Dr. Henry Poltz.) 112 S. mit 51 Abb. Wolfenbüttel- 
Hannover 1948. Wolfenbütteler Verlagsanstalt G. m. 
b. H. (Notdruck.) Preis brosch. 7,50.DM. 


Dr. Ing. habil. Hermann Athen, Vektorrech- 
nung. (Bücher der Mathematik und Naturwissen- 


‚ schaften. Herausgeber: Dr. Henry Poltz.) 90 S. mit 


60 Abb. Wolfenbüttel-Hannover 1948. Woltenbüt- 
teler Verlagsanstalt G. m. b. H. (Notdruck.) Preis 
brosch. 6,50 DM. 


. brosch, 4,80 DM. 


Hochschule Hannover), Einführung in ’die 
Schalenstatik. (Bücher der Technik, Heraus- 


geber; Dr. Ing. Alfred Kuhlenkamp.) ‘918. mit 
54 Abb. Wolfenbüttel-Hannover 1948. Wolfenbütte- 


ler Verlagsanstalt G. m. b. H. (Notdruck.) Preis 
a { | es ‚ 
Dr. Wolfgang Haack (vorm. o. Prof. an der Techni- 
schen Hochschule Berlin-Charlottenburg), -Diffe- 
rential-Geometrie Teil ’II (Bücher der 
Mathematik und Naturwissenschaften. Herausgeber 
Dr. Henry Poltz.) 131 S. mit 3 Abb. Wolfenbüttel- 
Hannover 1948. Wolfenbütteler Verlagsanstalt G. m. 
b.H. (Notdruck.) Preis brosch. 7,— DM, AM 


Prof. Dr. C. Schaefer (Univ. Köln), Prof.Dr.L.Brg- 


mann, Grundaufgaben des physikali- 
schen Praktikums, bearbeitet von Prof. Dr. 


‚W. Kliefoth. VI+ 2378. mit 180 Abb. und: 


18 Tabellen. Leipzig 1948.: B.G. Teubner Verlags- 


„gesellschaft. Preis brosch. 5,80 DM. 


Dr. Herbert Rost, Die Valorisation der 
Materie, VIII + 2998. mit 8 Abb. 'S. Hirzel- 
Verlag. Leipzig 1948. Preis geb. 20,— DM. 


Dr. Ing. Hayo Föppl, Messung der Span- 
nungen und der plastischen Verfor- 
mungen an oberflächengedrückten 
Probekörpern. (Mitteilungen des Wöhler-In- 
stitutes Braunschweig, Heft 41.) 67.8. mit 40 Abb. 
Braunschweig 1948. Friedrich Vieweg u, Sohn, Preis 
brosch. 4,— DM. 


NACHRICHTEN 


Tagung der Gesellschaft für angewandte Mathematik 
und Mechanik in der britischen Besatzungszone in 
Göttingen vom 22. bis 24. 9. 1948. In der Zeit vom 20, 
bis 22. September fand in Göttingen die erste Tagung 


der Gesellschaft für angewandte Mathematik und Me- | 
. chanik in der britischen Besatzungszone statt, an. der 


etwa 70 Mathematiker teilnahmen. Unter ihnen. be- 
fanden sich auch eine Reihe von Gästen aus der ameri- 
kanischen und der russischen Besatzungszone. Nach 
einem zwanglosen Treffen am Dienstag Abend begann 
am Mittwoch Vormittag die eigentliche Arbeitstagung, 
die durch Begrüßungsworte von Prof. Prandtl 
eingeleitet wurde. Der erste Tag brachte zunächts 
anter dem Vorsitz von Prof. Tollmien, Göttingen, 


dann von Prof. E Schmidt, Braunschweig, Tol- 
gende Vorträge aus dem Gebiet der Strömungslehre:; 


F.Vandrey, Göttingen: Die Reflexion schwacher 
Störungen an Unstetigkeitsflächen einer ebenen 
Unterschallströmung. 

E.Krahn, Göttingen: Näherungsverfahren zur Be- 
rechnung. kompressibler Unterschallströmungen. 

K. Oswatitsch, Emmendingen: Gesetzmäßig- 
keiten der schallnahen Strömung. 

M. Schäfer, Göttingen: Über die Berechnung der 
Ausbreitung von Störungen in einer Überschall- 
strömung nach der direkten Hodographenmethode. 

L. Prandtl, Göttingen: . Erzeugung von Zirku- 
lation beim Schütteln von Gefäßen, 


‘A. Betz, Göttingen: Neuere Geräte für Strömungs- 


LT 


\ ‚die Mitgliederversammlung der 


1. Gi inze ei RER 
schichten. 
messungen. 


H.Fromm, Alzey: Zur Theorie der Gefügebildung 
beim Erstarren von Metallen 


' M.Kobler, Horb/N.: Der Begriff der statistischen EB 


GEL drodynamik. 

F. Ri egels, Göttingen: Dreidimensionale - Strö-. 
mung um schlanke, nicht rotationssymm 
Körper mit ‚gerader Achse. 92 

H. Ludwieg und W. Tillmann, Göttingen: 
Untersuchungen über die Wandechuhppanmuig 
turbulenten Reibungsschichten. 

H.Reichardt, Göttingen: Zur Frag 
spannungsmessung in turbulenter «Strömung, - 

J. Rotta, Göttingen: Neue Rechnungen zur Sta- 
tistik der isotropen Turbulenz. s 


Am Vörmittag des folgenden Tages fand zunächst 
Ilschaft statt. 
Prof. Prandtlerinnerte daran, daß die Gesellschaft 
vor 25 Jahren in der schwersten Inflationszeit in Mar- 
burg gegründet wurde. Es erfolgten sodann: die Wah- 
len. Gewählt würden: 


zum Vorsitzenden: Prof. Prandtl, Göttingen, 


zum stellv.}Vorsitzenden: Prof. E. Schmidt, 
Braunschweig, 

zum Schriftführer: Prof. G. Voge 1 pohl, Göt- 
tingen, 


zu wissenschaftlichen Beiräten: Prof. C. Weber, 
Schlewecke,‘ Prof. Betz, Göttingen, Prof. 
Tollmien, Göttingen, 

zu Kassenprüfern: Prof, C. Weber, Bhkrere: 
Dr.-M. Schäfer, Göttingen. 


Nach längerer Diskussion wurde der Mitgliederbei- 
trag an 3,— DM. jährlich festgesetzt. 

Dann übernahm Prof. Walther, Darmstadt, für 
die folgenden Vorträge aus dem Gebiete der Mathe- 
matik den Vorsitz: 


F. Emde, Stuttgart: Über komplexe Nullstellen 
von Zylinderfunktionen. 

F. Tölke, Karlsruhe: Quasistrenge geschlossene 
Formeln für elliptische Normalintegrale, Theta- 
funktionen und elliptische Funktionen. 

H.Schäfer, Braunschweig: Transformationen der 
Variationsrechnung und ihre Anwendung auf tech- 
nische Eigenwertpröbleme. 

H.Bückner, Göttingen: Abschätzung von Lösun- 
gen inhomogener linearer Differentialgleichungen. 


*H, Gebelein, Braunschweig: Einige Bemerkun- 


gen über die Orthogonalpolynome zu einem end- 
lichen Intervall mit einer vorgeschriebenen Bele- 
gung. 

An diese Vorträge schloß sich eine Besichtigung des 
Max-Planck-Institutes für Strömungs- und Reibungs- 
forschung an. Der Nachmittag vereinigte die Teil- 
nehmer zu einer Wanderung nach Nikolausberg. 

Der Freitagmorgen war im wesentlichen Fragen der 
Festigkeit und Elastizität gewidmet, Unter dem Vor- 
sitz von Prof. Pöschl, Karlsruhe, sprachen: 
K.Ludwig, Hannover: Die Biegung der Rechteck- 

platte ohne die Bernoullische oder andere An- 
nahmen. 

E.Pestel, Hannover: Träger mit wandernder Last 
(strenge Lösung). 

O0. Föppl, Braunschweig: Die : “igenspannungen, 
dieim Kern und in der äußeren Schicht von ober- 
flächengedrückten zylindrischen Stäben erzeugt 
werden. 


zur Berechnung laminarer use ren he 


der Schub- N 


Air 
A. Walz, Teningen (Baden): Ein 
Gerät für harmonische Analyse und Ye 
An fast alle Vo hlossen sich aı ‚oft, 
sehr a D = Bene leider N 
lich wegen. Mangels an Zeit 
ten. Ein Bierabend, rege Mittag- und Abend- 
essen sowie die Wanderung nach Nikolausberg een. ni 
lichten eine Weiterführung der Me a 
nen Gruppen und eine em en 
Fragen im kleineren Kreise. Berichte über die einzel- 
nen Vorträge werden in einem der nächsten Hefte der. j 
Zeitschrift erscheinen. Die Tagung we dank der sorg- 2 
fältigen Vorbereitung durch die Herren Porleepagae RT 
Prandtl und Vogelpohl ganz ausgezei 
gelungen, und ich glaubeim Namen aller Teilnehmer 
zu sprechen, wenn ich diesen Herren für ihre mühe- \ 


d Li 


vollen Vorarbeiten auch hier den herzlichsten Dank SS“ 
ausspreche. 1 Bei 
Dresden. ; 


Willers. 


Braunschweig: Herr Pröf. Dr. Hermann 
Blenkerhielt Lehraufträge für Praktische Mathema- 
tik und Differentialgeometrie. OR 
Herr Prof. Dr. Rudolf Weyrich erhielt einen ER 
Lehrauftrag für Mathematik, B 
Herr Prof. Dr. Hermann Schmidt erhielt die 
venia legendi und wurde als Privatdozent für Mathe- 
matik zugelassen, 


Druckfehlerberichtigung. 


W. Wuest : Zur Theorie des gegabelten Verdich- 
tungsstoßes, Z. ang. Math, Mech. 28 (1948) S. 73. 


1 1 
In Formel (2) ist Ba statt I zu lesen, in. 
Formel (8) ist 7, und 7, zu vertauschen, in den Glei- 


chungen (13) und (14) ist auf der linken Seite =. i 
statt zT mu setzen. In der anschließenden Formel 
für P ist im zweiten Zählerglied 1/(1 + 27) statt 
1(M?} — au — 1) zu setzen. Außerdem sind in den 
Formeln für P und Q@ die vor dem Bruchstrich stehen- 
den Vorzeichen umzukehren, Formel (15) liefert 
M, = 1,2447 statt 1,2456. 
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